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Measurements of the Thermal Conductivity and Electrical 
Resistivity of Platinum from 100 to 900 °C 


D. R. Flynn and M. E. O’Hagan 


Institute for Applied Technology, National Bureau of Standards, Washington, D.C. 20234 


(August 17, 1967) 


Measurements have been made of the thermal conductivity and the electrical resistivity of com- 
mercial grade platinum (99.98% pure) in the temperature range 100 to 900 °C. The measurements have 
been made with a view to providing accurate data on the thermal conductivity of platinum to serve 
as a basis for establishing platinum as a thermal conductivity standard reference material. Two meth- 
ods of measuring the thermal conductivity have been employed, one an electrical method and the 
other a nonelectrical method. In the electrical method, a direct current passed through a necked-down 
portion of the specimen and the thermal conductivity was determined in terms of the temperature 
and electrical potential distributions in the necked-down region. The second method was of the abso- 
lute guarded longitudinal heat flow type. The experiment was designed to permit measurements by 
both methods in the same apparatus and on the same specimen thereby providing as direct a compari- 
son as possible between the methods. The data given by the two methods agree within experimental 
error and show the thermal conductivity of platinum to be a smoothly increasing function of tempera- 
ture in the measured range. Additional measurements on samples of differing purities are necessary 
before platinum could be adopted as a thermal conductivity reference material. 


Key Words: Conductivity, electrical conductivity, electrical resistivity, heat conductivity, Lorenz 
function, platinum, reference material, resistivity, standard, thermal conductivity. 
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1. Introduction 


Standards and standard reference materials are the 
basis of a consistent and accurate measuring system. 
The need for standard reference materials in thermal 
conductivity measurements is two-fold. In the first 
place, such materials are required for comparative 


measurements in which the thermal conductivity of 
the material under test is determined in terms of that of 


the standard reference material. Secondly, such ma- 


terials are required in evaluating the accuracy of 


apparatus designed for thermal conductivity measure- 
ments. The degree to which the measured value of the 
thermal conductivity of the standard reference ma- 
terial agrees with the accepted value is a check on the 
accuracy of the apparatus in which the measurements 
were made. 

The basic requirements for any standard reference 
material are that it be stable, reproducible and appro- 
priate for the measurements at hand, and that the 
property in question be uniform throughout the ma- 
terial. In the case of standard reference materials for 
thermal conuctivity other desirable requirements are 
that the standard be usable over a wide range of tem- 
perature, that it be chemically inert so as not to be 


*This research was the subject of a dissertation submitted by M. E. O'Hagan to the 
Faculty of The School of Engineering and Applied Science, The George Washington 
University, in partial satisfaction of the requirements for the degree of Doctor of Science 
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who supplied the platinum and financial assistance in the form of a grant to the National 
Bureau of Standards. M. E. O'Hagan was supported in part by Grant AFOSR-1025-66 
from the Air Force Office of Scientific Research and in part by the Institute for Industrial 
Research and Standards, Dublin, Ireland. 


**Present address: Institute for Industrial Research and Standards. Dublin. Ireland. 


affected by or affect other materials in the system and 
that the thermal conductivity of the reference material 
be close in value to that of the materials which are to 
be measured in terms of it. 

The advantages of using platinum as a thermal con- 
ductivity reference material have been pointed out by 
Powell and Tye [1]! and by Slack [2]. Platinum is 
available in high purity in pieces of substantial size. 
It has a fairly high melting point (1769 °C on the 1948 
International Practical Temperature Scale), has no 
known transition points, and is relatively stable chem- 
ically in air and other atmospheres, with the excep- 
tion of hydrogen, even at high temperatures [3, 4]. 
Its thermal conductivity, although relatively high for 
use as a reference material with nonmetals, is about 
the geometric mean for metals and alloys. 

Since the thermal conductivity of a pure metal is 
strongly correlated with the electrical conductivity of 
the metal, it is highly desirable that the electrical con- 
ductivity of a metal which is intended for use as a 
thermal conductivity reference material be very 
stable under varying heat treatments. The stability of 
the electrical conductivity of platinum is evidenced by 
the fact that the International Practical Temperature 
Scale is defined by a platinum resistance thermometer 
in the temperature range — 182.97 to + 630.5 °C [5, 6]. 
Studies are currently underway at NBS [7] and other 
laboratories to investigate the possibility of extending 
to the gold-point (1063 °C) the range over which a 
platinum resistance thermometer is used to define the 
temperature scale. 


'Figures in brackets indicate the literature references at the end of this paper. 
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Platinum appears to be, in every way save one, an 
ideal material to use as a thermal conductivity refer- 
ence standard. The exception is that the spread among 
the literature values for the thermal conductivity of 
platinum is considerable, to say the least. 

Powell, Ho, and Liléy [8] show a plot of essentially 
all of the published thermal conductivity data for 
platinum through the year 1965. The spread in the 
data increases from about 10 percent at room tempera- 
ture to over 30 percent at 1000 °C. Even if many of 
the older data are discounted, the picture is not par- 
ticularly improved. O’Hagan [9] has summarized the 
methods used for previous measurements of the 
thermal conductivity of platinum and also has sum- 
marized the characterizations of the various samples. 

Prior to the measurements of Powell and Tye [1]. all 
thermal conductivity values reported for platinum at 
temperatures above 100 °C were obtained by methods 
in which the temperature gradient in the specimen was 
produced by Joule heating due to passage of an electric 
current directly through the specimen. While the re- 
sults of all but one [10] of these higher temperature 
investigations employing electrical methods essen- 


tially agree with one another, they disagree with those | 
of later investigations [1], 11, 12], which employed - 


nonelectrical methods. This raised the question as to 
whether or not electrical methods yield results that are 
intrinsically different from those of nonelectrical 
methods. This could mean that the theory which has 


been .used in analyzing electrical methods is in error, : 


or it could mean that heat conduction is significantly 
dependent on electric current density, at least in the 
case of platinum. 

The considerations discussed above pointed to the 
need for a comprehensive investigation of the thermal 
conductivity of platinum. It was felt that both an 
absolute steady-state method without an electric cur- 
rent flowing in the specimen, and also an absolute 
steady-state method with a current flowing in the 
specimen should be employed. 

As regards the nonelectrical method, experience at 
NBS with guarded longitudinal heat flow methods in- 
dicated that such a method could be made to yield 
accurate results on a material having as high a thermal 
conductivity as platinum, provided a specimen of 
sufficient cross-sectional area was used. With a fairly 
conductive metal, there were no particular advantages 
in going to a radial heat flow method: furthermore, to 
do so would have required a much larger sample. 

As regards the electrical method, an arrangement 
utilizing quite large current densities would be more 
likely to reveal deviations due to a dependence of 
thermal conductivity on current density. It was also 
desirable to use the same method as that used by most 
of the previous investigators. Fortunately, these two 
desiderata both pointed to the necked-down sample 
configuration utilized for measurements on platinum 
by Holm and Stormer [13], by Hopkins [14], and by 
Cutler, et al. [15]. 

To give a direct and accurate comparison between 
the two methods it was considered desirable to com- 
bine both sets of measurements in one apparatus and 


on the same specimen thereby eliminating a number of 
uncertainties which would arise irg comparing data 
derived from measurements in different apparatus 
and on different specimens. The above considerations 
led to an apparatus, described in section 3, in which 
thermal conductivity measurements can be made by 
both the usual longitudinal heat flow method and by an 
electrical method. 

It was also felt that thermal conductivity measure- 
ments should be made on platinum samples of at least 
two purities. Samples were obtained of a high purity 
platinum (resistance thermometer grade) and of a 
somewhat lower purity platinum (commercial grade). 
As of this writing, only the measurements on the lower 
purity platinum sample have been completed. The 
results obtained on that sample are presented in this 
paper and are compared with the results of other 
investigators. 


2. Description of Sample 


In order for a valid comparison to be made between 
the results of different investigators who measure on a 
particular kind of material, it is necessary that their 
specimens be characterized as extensively as possible 
so that differences in specimens may be accounted for. 
For this reason a number of pertinent measurements 
were made in an attempt to characterize the speci- 
men used in the present investigation. These meas- 
urements and the results thereof are described below. 

The platinum was provided by Englehard Industries, 
Inc., in the form of a solid bar 2.04 cm in diameter by 
31 ecm long, and was classified as being of commercial 
purity. The fabrication and cleaning procedures used 
in preparing the platinum have been described by 
O’Hagan [9]. 

The as-received bar was annealed in air for 5.5 hr at 
770 °C in a horizontal tubular furnace and furnace 
cooled at a rate of approximately 120 deg/hr. Shortly 
thereafter the bar was accidentally dropped causing 
it to deform slightly at one end. After correcting the 
damage the bar was reannealed for 1.5 hr at 680 °C 
and furnace cooled at a rate of approximately 90 deg/hr. 
Its thermal conductivity was then measured in the 
NBS Metals Apparatus [34, 35] over the temperature 
range — 160 to +810 °C [32]. 

Following these measurements the bar was ma- 
chined and ground to 25.4 em long by 2.000 cm 
diam. The density of this bar was measured (see 
below) and the electrical resistance was measured 
at ice and liquid helium temperatures (see below). 
The thermal conductivity specimen (18.4 cm long by 
2.000 cm diam) was then fabricated from one end of 
this bar. A length of 6.4 cm was cut from the other 
end for separate low temperature thermal conductivity 
measurements [33]. The remaining disk, approxi- 
mately 1 cm long, was reserved for metallographic and 
spectrographic analyses. 

A thin neck, approximately 0.1 cm in diameter and 
0.3 cm long, was machined in the thermal conductivity 
specimen at approximately 4 cm from one end. To 
cleanse this necked-down region of oil and any con- 
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tamination contracted from the cutting tools, the fol- 
lowing cleaning procedure was followed. Degreasing 
was effected by inimersion for half an hour in tri- 
chlorethylene vapor. The neck was then pickled for 
10 min in hot 50 percent nitric acid. Following this 
the neck was washed in distilled water, pickled for 
10 min in 50 percent hot hydrochloric acid, and again 
washed in distilled water. There was equal likelihood 
of the rest of the specimen having slight surface con- 
tamination as a result of machining but its effect on 
the bulk properties of the specimen would not be nearly 
as grave as the corresponding effect in the neck region. 
It was considered sufficient to clean the surface of the 
specimen with toluene and carbon tetrachloride. 

A general qualitative spectrographic analysis, per- 
formed by the NBS Spectrochemical Analysis Section 
on the 2-cm-diam by 1-cm-long platinum disk men- 
tioned above, detected Ag (10-100 ppm), Pd (10-100 
ppm), Fe (< 10 ppm), and Mg (< 10 ppm). O’Hagan [9] 
reported the results of a quantitative spectrographic 
analysis made on a portion of a 0.05 cm wire drawn 
from the same platinum ingot as the thermal con- 
ductivity specimen. 

Photomicrographs were also made of the platinum 
disk: cut from the bar sample. They showed grain size 
to be of the order of 0.02 cm. No evidence appeared 
of inclusions of foreign matter or of any irregularities 
in microstructure. Hardness measurements were made 
on the original platinum bar, after the second anneal, 
using a Vickers Pyramidal Diamond Tester with a 10 kg 
load. Values ranged from 36.5 to 38.0 Vickers hardness 
number 

The density of the bar, when it was 25.4 cm long, 
was determined by mass and dimensional measure- 
ments to be 21.384 g/em® at 21 °C, accurate to within 
+ 0.002 g/em*. This density is an average value for the 
whole bar and there is no guarantee that the density 
was uniform to that degree throughout the bar. 

The ratio of the resistance at the ice-point tempera- 
ture to that at the boiling point of helium at atmospheric 
pressure is a measure of the extent and condition of 
impurities in a material and of the crystallographic 
state of the material. This ratio was determined on the 
2-cm-diam bar. 

The ice-point resistances were measured both be- 
fore and after the helium point measurements. The 
current was supplied from a regulated d-c power 
supply and was measured using a calibrated resistor 
and a precision potentiometer. The voltage drops in 
the specimens were measured on a high precision 
6-dial potentiometer. In each case the resistance was 
measured at three or four different current levels and 
the value corresponding to zero current obtained by 
extrapolation. The ratio of the resistance of the sample 
at the ice-point to that at the helium-point was found to 
be 393. This value is believed to be accurate to within 
1 percent. However, it corresponds to an average 
value over a considerable length of the sample and the 
sample may not have been uniform in purity throughout. 

A set of knife edges of known separation was fas- 
tened to the 2-cm bar during the first set of ice-point 
resistance measurements. The knife edges acted as po- 
tential taps. Using the known separation of the knife 


edges and the cross-sectional area of the bar, the ice- 
point resistivity, corrected to 0 °C dimensions, was 
determined to be 9.84740 ecm, accurate to within 
+ 0.010402 cm. 

A length of 0.05 cm platinum wire, drawn from the 
same ingot as the thermal conductivity sample, was 
electrically annealed in air for 1 hr at about 1450 °C. 
The electromotive force of this wire versus the plati- 
rum standard Pt 27 [16] was measured by the NBS 
Temperature Section with the reference junctions at 
0 °C. The values obtained (at 100 deg intervals) in 
creased in an essentially linear manner from 0 pV at 
0 °C to +15 pV at 1100 °C. This indicates that the 
sample used in the present investigation was less pure 
than Pt 27. 

The temperature coefficient of resistance, 
a@=(Rioo — Ro)/100Ro, between the ice-point and the 
steam-point is often used as an indication of the purity 
of resistance thermometer grade platinum. The limiting 
value of q@ for extremely pure platinum is given by 
Berry [17] to be 0.0039285. The size and low resistance 
of the thermal conductivity specimen precluded a 
highly accurate direct determination of @ on the speci- 
men using existing equipment. On the basis of the 
electrical resistivity measurements (described later) 
on the necked-down portion of the specimen, a had a 
value in the range 0.003876<a<0.003916. The rela- 
tively large uncertainty in @ arises from the use of 
platinum versus platinum—1l0 percent rhodium ther- 
mocouples to measure the temperature near 100 °C. 
Berry [17] gives a plot which correlates a with the ratio 
of the resistance of a sample at absolute zero to that 
at the ice-point. He also correlates the resistance at 
absolute zero with that at the helium-point. On the 
basis of Berry’s correlations, the ice-point to helium- 
point resistance ratio for the thermal conductivity 
specimen used in the present investigation corre- 
sponds to 0.003907<a<0.003916. This is not incon- 
sistent with the range of values obtained from the 
electrical resistivity measurements. 

Although @ could have been measured directly with 
high accuracy on the 0.05 cm platinum wire, this was 
not done since the results would not necessarily be 
valid for the 2 em bar, owing to possible differences in 
purity and annealing. Corruccini [18] gives an em- 
pirical expression, due to Wm. F. Roeser, correlating 
a with the electromotive force versus Pt 27 with the 
junction at 1200 °C and the reference junctions at 
0 °C. Extrapolation of the emf measurements men- 
tioned above, indicates an emf of +16 wV at 1200 
C, corresponding on the basis of Roeser’s expres- 
sion, to a ~ 0.003914 for the 0.05 cm wire drawn from 
the same material as was used to fabricate the thermal 
conductivity specimen. 

\ cooperative project, involving NBS and several 
producers of thermometric grade platinum, is cur- 
rently underway to study the properties of pure plat- 
inum. If,as a result of this project, it appears that 
further characterization is indicated for the platinum 
used in the present investigation, such characteriza- 
tion will be performed on material which is being 
reserved for that purpose. 
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3. Method and Apparatus 


3.1. Method 


As discussed in the introduction, it was decided to 
build a single apparatus in which thermal conductivity 
could be measured by both an absolute guarded longi- 
tudinal heat flow method (nonelectrical method) and 
by a method in which a necked-down portion of the 
sample was heated directly by passage of an electric 
current (electrical method). The specimen configura- 
tion selected for these measurements is shown in fig- 
ure 1. The specimen (A) was raised to the desired 
temperature level by means of the heaters Q; and 
Q:. In the longitudinal heat flow method, the heater, 
Q., located slightly above the center of the bar pro- 
duced a temperature gradient along the portion of the 
bar below Q».. Heat from this heater was prevented 
from flowing up the bar by adjusting Q3 so that there 
was negligible temperature difference across the 
necked-down region of the specimen. Lateral heat 
losses from the bar were minimized by matching the 
temperature distribution along the guard to that along 
the specimen. Thermal conductivity was calculated 
from the measured temperature distribution along the 
lower portion of the bar, the power input to the central 
heater, and the geometry. 
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SCALE, cm 


x — THERMOCOUPLE 

<— CONTROL COUPLE 

@— HEATER WINDING 
Schematic 


FIGURE 1. diagram showing locations of heaters and 


thermocouples in the apparatus. 


In the electrical method, where the sample was 
directly heated by passage of an electric current, 
the voltage drop across the necked-down region of the 
specimen was measured as a function of current, while 
the maximum temperature rise in the neck was com- 
puted from the change in electrical resistance (due to 
a given change in current) and the temperature coef- 
ficient of resistance of the material. The thermal 
conductivity was determined from the voltage drop 
across the neck, the computed maximum temperature 
rise in it, and the electrical resistivity of the material. 

The apparatus is described in detail in this section. 
The experimental test procedures and calculation pro- 
cedures for the nonelectrical method are described in 
sections 4.1 and 4.2: those for the electrical method 
are described in sections 5.1 and 5.2. 


3.2. Mechanical Configuration 


The mechanical configuration of the apparatus is 
illustrated diagrammatically in figure 2 and described 
in detail below. 

a. Specimen 

The specimen (A) was a bar 2 cm in diameter by 
18.4 cm long with a 0.11 em diam by 0.33 cm long neck 
machined in it 4.1 em from the upper end. A special 
technique, described by O’Hagan [9], had to be de- 
veloped for machining the neck due to its structural 

















SCALE, cm 


FIGURE 2. The apparatus (components are identified in the text). 
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weakness. Hollow molybdenum extensions (B) of the 
same diameter as the specimen were screwed to the 
specimen (A) at both ends. The open ends were brazed 
to copper blocks (C) that served as heat sinks. The 
molybdenum extensions were filled with high purity 
“coral” alumina. The lower end of the specimen as- 
sembly (B-A-B) was bolted to, but electrically in- 
sulated from, a brass flange which was welded to a 
water-cooled brass column (D). This column was 
firmly bolted to plate (E) which served as a base for 
the apparatus. 

One of the major problems with the necked-down 
specimen was that of protecting the neck from mechan- 
ical strain due to tension, compression, torsion, or 
bending. Any clamp supporting the neck would have 
had to be electrically insulated from the specimen and 
differential thermal expansion between the clamp and 
specimen could have introduced strain in the neck. 
As an alternative to a clamp it was decided to counter- 
balance the load on the neck due to the weight above 
it so that there would be only a small net force on the 
neck. The weight of platinum above the center of the 
neck was computed from dimensional measurements 
and from the measured density of the specimen, and 
the components extending from the upper end of the 
specimen were weighed before assembly. 


The counterweight (W) was suspended from a string 
which passed over two pulleys and attached to an 
aluminum hanger (H). The pulley wheels were mounted 
on low-friction bearings having a starting force’of less 
than 1 g each. A molybdenum well (F), which was 
brazed to the upper copper block (C), passed through 
a linear bearing (L) which served to maintain the up- 


per part of the specimen in precise alinement with the 
lower part, and presented negligible resistance to the 
free vertical motion of the specimen resulting from 
thermal expansion. This bearing was mounted on the 
upper plate (E) but electrically insulated from it. The 
two aluminum plates (E, E) were connected by three 
tie bars (K) to form a rigid framework for maintaining 
proper specimen alinement. An auxiliary device, de- 
scribed by O'Hagan [9], prevented the upper part of the 
specimen from rotating but still allowed free vertical 
motion. Current was introduced to the specimen via a 
copper rod (N) at the lower end and through a hollow 
molybdenum electrode (QO) at the upper end. The upper 
electrode was brazed to a copper support (S) which was 
fastened to, but electrically insulated fram, the upper 
plate (E). The molybdenum well (F) contained a liquid 
metal alloy into which the electrode dipped thereby 
affording a flexible current connection. The buoyant 
force of the liquid metal on the electrode contributed 
to the load on the neck and was compensated for in the 
counterweight. The electrode was fixed but the molyb- 
denum well moved upwards with the specimen due to 
thermal expansion during test runs. This changed the 
buoyant force and consequently put a load on the neck. 
The maximum change in buoyant force was only 3 g, 
however, which would not strain the neck signifi- 
cantly. The liquid metal used was a gallium-indium 
eutectic alloy chosen primarily for its low vapor pres- 


sure and its comparatively low freezing temperature of 


15.7 °C. The requirement for low vapor pressure was 
dictated by a need to evacuate the system. Preliminary 
tests were run to evaluate the uncertainty in buoyant 
force due to surface tension and sticking of the gallium- 
indium to the molybdenum surfaces. A very definite 
hysteresis effect was observed as the electrode was 
moved relative to the well and then returned to its 
initial position. The largest uncertainty in buoyant force 
was determined to be about 4 g. The choice of molyb- 
denum as the electrode and well material stemmed 
from its compatability with gallium, which reacts with 
most other metals, and from the fact that molybdenum 
is wetted by gallium. The current feed-in system also 
served as a heat sink for the upper part of the speci- 
men assembly. The hollow molybdenum electrode was 
internally cooled by circulating water at a temperature 
higher than the freezing point of the gallium-indium 
eutectic alloy. 


b. Furnace and Guard 

The inner core (G), or the guard as it is called, was 
a molybdenum tube of 5.7 cm inside diameter. Since 
the inner core acted as a thermal guard to prevent 
heat losses from the specimen, it was considered more 
desirable to make it from metal rather than from ce- 
ramic so that the temperature distribution along it 
could be more easily controlled and more accurately 
measured. The bottom of the guard was attached to a 
water-cooled brass plate (P). A water-cooled brass ring 
(Q) was attached to the upper end of the guard. The 
outer furnace core (V) was an aluminum oxide tube 
supported top and bottom by three l-cm diam alumi- 
num oxide rods (U). 

The exterior portion of the furnace consisted of a 
water-cooled shell (X) supported between two water- 
cooled plates (P). Attached to the upper plate of the 
furnace was a split nut. This nut engaged a lead screw 
mounted between the plates (E, E) and by turning the 
screw the furnace could be moved up or down. Ready 
access to the specimen was thereby afforded. The 
three rods (K) acted as guide rods for the furnace. 
Six linear bearings (Y) attached to the furnace plates 
ensured alinement and permitted the furnace to move 
up and down freely. The correct vertical location of the 
guard relative to the specimen was determined when 
a probe attached to the upper end of the guard made 
electrical contact with a plate attached to the molyb- 
denum well at the ypper end of the specimen assem- 
bly. When the furnace had been positioned, the indi- 
cating probe attached to the guard was removed. 

The space between the specimen and the molyb- 
denum guard and that between the guard and the 
water-cooled shell were filled with fine high-purity 
aluminum oxide powder of low thermal conductivity 
and low bulk density (0.16 g/cm’). 


c. Environmental System 
The entire apparatus was mounted inside a 24-in 
diam metal bell jar to enable operation in an inert 
atmosphere. A 4-in oil diffusion pump and a 5 cfm 
mechanical pump were used to evacuate the system 
prior .o refilling with argon or helium. Initial evacua- 
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specimen heater. 


Schematic diagram construction of the 


tion was controlled to avoid disturbance of the fine 
powder insulation. 


3.3. Thermal Configuration 


The heater and thermocouple locations on the 
specimen and guard are shown in figure 1. 


a. Specimen 


The heaters, Q; and Q;. used to raise the mean tem- 
perature of the sample were located in the molybde- 
num extensions. These heaters consisted of six 
series-connected helical coils of 0.02 em diam plat- 
inum—thirty percent rhodium wire insulated from the 


surrounding metal by thin-wall aluminum oxide tubing. 
Platinel ? thermocouples were attached adjacent to 
heaters Q, and Q, for use in controlling the tempera- 
tures at these locations. 


The central heater, Q.. was contained in six holes 
drilled through the platinum bar and was constructed 
as shown in figure 3. The inner four holes were 0.1 cm 
in diameter and accommodated helical elements con- 
tained in thin-wall aluminum oxide tubing. The ele- 
ments were made from 0.013 cm diam platinum—10 
percent rhodium wire, the outside diameter of the helix 
being 0.05 cm and the pitch 0.025 cm. The outer two 
holes were 0.16 cm in diameter and accommodated 
“swaged elements” having platinum—l0_ percent 
rhodium sheaths insulated from platinum — 10 percent 
rhodium heater wires by compacted MgO powder in- 
sulation. The swaged elements were a snug fit in the 
holes so that there was good thermal contact between 
the sheath and the bar, and consequently good thermal 
coupling between the heater and the bar. The six ele- 
ments were connected in series as in heaters Q; and 
Q;. Platinum heater leads, 0.05 cm in diameter, 
were welded to the ends of the swaged elements. 
The good thermal contact in the swaged elements en- 
sured that the temperature at the ends of the heater 


Platine! has a high thermal emf, approximately that of Chromel P versus Alumel 


The negative leg of the thermocouple is 65 percent Au, 35 percent Pd alloy (Platinel 5355) 
and the positive leg is 55 percent Pd, 31 percent Pt. and 14 percent Au (Platinel 7674) 


closely approximated that of the specimen. Moreover, 
the current leads extended radially from the heater in 
an isothermal plane. The combined result was to mini- 
mize heat losses via the leads. Two 0.02 cm platinum— 
10 percent rhodium potential leads were welded to each 
of the current leads, one at the junction of the heater 
and the current lead, and the other about 1 cm back 
along the current lead. The two platinum—10 percent 
rhodium potential leads together with the intervening 
section of platinum current lead served as a differential 
thermocouple to determine the temperature gradient 
in the current lead. By taking potential readings with 
the current flowing in the forward and reverse direc- 
tions, the IR drop in the current leads could be ac- 
counted for and the temperature gradients therein 
determined. These data were used in computing heat 
flows along the leads. The potential drop across the 
inner taps was used in computing the power generated 
in the heater. The distances from the heaters to the 
nearest thermocouples and potential taps were such 
that perturbations in heat flow and electric current 
flow generated by the presence of the heaters decayed 
to an insignificant level at the position of the thermo- 
couples or potential taps (see appendix B of O’Hagan 
[9}). 

Five thermocouples, spaced 2 cm apart, were lo- 
cated in the gradient zone of the specimen, with the 
lowermost one (designated 4 in fig. 1) being 2 em from 
the end of the specimen. The thermocouples were 
fabricated from 0.020 cm diam platinum and platinum — 
10 percent rhodium wires which were annealed 
in air at about 1450 °C for 2 hr and then butt-welded 
together. They were pressed into 0.018 em wide by 
0.023 cm deep horizontal slits in the surface of the 
specimen thereby replacing the metal removed in 
machining the slits. By virtue of the fact that the speci- 
men was fairly pure platinum with essentially the same 
absolute thermoelectric power as the platinum leg of 
the thermocouple the junction of each thermocouple 
was effectively at the point where the platinum—10 
percent rhodium leg first made contact with the speci- 
men, and the temperature measured was the tempera- 
ture at that point. The platinum—10 percent rhodium 
wire emerging from its groove extended a short way 
around the specimen in the same isothermal plane 
insulated from the specimen in broken ceramic tub- 
ing—so as to minimize the amount of heat conducted 
away from the junction. Similar thermocouples were 
located in the molybdenum extensions, three in each, 
to measure the temperature distribution along them. 
This information was essential to the mathematical 
analysis of the system. 

Additional thermocouples were located on either 
side of the neck (locations 9, 10, 11, and 12 in fig. 1). 
These were fabricated from annealed 0.038 cm diam 
platinum and platinum—10 percent rhodium wire and 
pressed into slits. In addition to measuring tempera- 
ture, these thermocouples were wired at the selector 
switches so that the platinum legs could be used to 
measure voltage drops across the neck when an electric 
current was flowing through the neck. With no current 
flowing, the platinum—1l0 percent rhodium legs, in 
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conjunction with the platinum neck, could be used as a 
differential thermocouple to control the differential 
temperature across the neck in the longitudinal heat 
flow method of measurement. 


b. Guard 


The guard had three heaters, Q;, Q;, and Qg, at 
positions corresponding to those on the specimen 
assembly. All three were swaged heaters with plat- 
inum—l0 percent rhodium sheaths and heating ele- 
ments, and MgO insulation. They were pressed into 
grooves machined in the guard, thus giving good ther- 
mal contact. The end heaters (Qs and Qs) were used to 
keep the guard at the desired temperature. The cen- 
tral heater (Q;) was used to produce a temperature 
gradient in the guard matching that in the specimen. 

The guard was electrically grounded but the heaters 
were isolated. Platinel control thermocouples were 
peened into the guard adjacent to each of the heaters. 

Twelve thermocouples were located on the guard, 
three in the gradient zone, three in the isothermal zone, 
and three in each of the end zones. All the thermo- 
couples were 0.038 cm platinum versus platinum — 10 
percent rhodium with the junctions pressed into slits 
machined in the guard. The thermocouple wires were 
taken one turn around the guard in broken ceramic tub- 
ing in an isothermal plane and cemented to the guard 
with high purity alumina cement. This helped to temper 
the thermocouple leads and reduce the amount of heat 
conducted away from the junction by the leads. Within 
the furnace all the thermocouples were insulated in 
single-bore ceramic tubes. For the remainder of their 
lengths the wires were insulated in flexible fiber-glass 
sleeving. All the thermocouples, both from the guard 
and the specimen assembly, went to a junction box 
mounted on the inside of the feedthrough ring. There 
they were torch welded to identical wires which were 
taken through wax vacuum seals to an ice bath. All 
the Platinel control couples went to terminal strips en 
the upper plate of the furnace. There they were spot- 
welded to Chromel P and Alumel wires coming from 
the temperature controllers. 

The aluminum oxide outer core (V) was provided 
with a heater winding (0.1 cm diam molybdenum) to 
bring the furnace as a whole to temperature and to 
reduce heat losses from the molybdenum guard and 
the power load on its heaters. A Platinel control 
thermocouple was mounted on the outer core (V). 


3.4. Instrumentation 
a. Temperature Control 


In the longitudinal heat flow method of measuring 
thermal conductivity, the platinum—1l0 percent rho- 
dium legs of the outer pair of thermocouples in the 
neck region were used in conjunction with the necked- 
down portion of the specimen as a differential thermo- 
couple to control the power to heater Qs, and thus 
maintain essentially a zero temperature differential 
across the neck. The signal from this thermocouple 
was amplified by a chopper-stabilized d-c amplifier 


and fed into a:current-adjusting-type proportional con- 
troller incorporating automatic reset control and rate 
control. The output of the proportional controller regu- 
lated the power to the heater by means of a transistor- 
ized current amplifier fed by a regulated d-c power 
supply. In the electrical method of measuring thermal 
conductivity, an electric current flowed through the 
specimen and the above system of control could not 
be employed. In this case the Platinel control couple 
adjacent to heater Qs; was put in series opposition with 
a signal from an adjustable constant voltage source 
and the resultant signal fed to the proportional con- 
troller which regulated the power to Q3. The external 
signal was manually adjusted to give zero temperature 
differential across the neck. 

The specimen heater (Q2) was fed constant voltage 
(+0.01%) from a regulated d-c power supply. Power 
to heater Q; and to the three guard heaters (Q,, 
Q;. and Qs) was supplied by variable-voltage trans- 
formers, which in turn were fed by voltage-regulated 
isolation transformers. Power to each heater was 
regulated by individual thermocouple-actuated con- 
trollers. Power to the heater winding (Q;) on the 
alumina core was supplied by a variable-voltage trans- 
former fed by a voltage-regulated isolation transformer. 
The current was manually ratioed among the three 
heater sections. The total power to this heater was 
regulated by a single thermocouple-actuated con- 
troller. All heaters were supplied by separate isola- 
tion transformers or power supplies to minimize 
current leakage effects. 


b. Temperature Measurement 


The noble metal leads of the thermocouples were 
brought to an ice bath, where they were individually 
joined to copper leads. The copper leads. went in 
shielded cables to a bank of double-pole selector 
switches of the type used in precision potentiometers. 
The selector switches were housed in a thermally 
insulated aluminum box with 1 cm thick walls. The 
copper leads were thermally grounded to (but elec- 
trically insulated from) the switch box to minimize 
heat transfer directly to the switches. The emfs of the 
specimen thermocouples were read on a calibrated 
six-dial high-precision potentiometer to 0.01 wV, using 
a photocell galvanometer amplifier and a secondary 
galvanometer as a null detector. (Due to thermal emfs 
in the potentiometer and circuitry, these emfs were 
probably not meaningful to better than +0.05 pV.) 
The emfs of all other thermocouples were read on a 
second precision potentiometer to 0.1 wV using an 
electronic null detector. 


c. Power Measurement 


Power input to the specimen heater was measured 
using a potentiometer in conjunction with a high- 
resistance volt box to measure the drop across the 
inner set of potential taps, and a standard resistor in 
series with the heater to measure the current. 
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d. Resistance Measurement 


All electrical resistance measurements were made 
by measuring the current (from a 0-100 A regulated 
d-c power supply) flowing through the specimen, 
utilizing a calibrated 0.001 © standard resistor, and 
by measuring the appropriate voltage drop in the sneci- 
men, using a potentiometer. 


4. Longitudinal Heat Flow Method 
4.1. Experimental Procedure 


a. Preliminaries 

The furnace was heated to 150 °C and the system 
evacuated to 3X10-* torr. Initial pump-down was 
through a needle valve to give a sufficiently low rate so 
as not to disturb the very light powder insulation. 
When the pressure had fallen below 10-' torr, the dif- 
fusion pump was turned on, pumping initially through 
the needle valve, and, when the pressure was below 
10-2 torr, through the 4-in port. After pumping for 
24 hr, the pumps were turned off and the system back- 
filled with high purity (99.99%) argon. The argon was 
bled in slowly through a needle valve to avoid dis- 
turbing the powder. The pressure was allowed to build 
up to almost 1 atm before the valves were shut off 
and the cycle of evacuating and backfilling repeated. 
The final argon pressure after the second backfilling 
was about three-quarters of an atmosphere. 

The cooling-water flows to the system were adjusted 
to the desired levels as indicated on flow-meters. The 
water was pressure-regulated to maintain constant flow 
rate. A control thermocouple on the specimen was 
wired to shut all the heaters off if its temperature ex- 
ceeded a predetermined level. This was a safety pre- 
caution in the event of an interruption in the cooling- 
water flow. 

The test procedures are described in detail below. 
To facilitate the discussion let us refer to figure 1. 
The region of the specimen below the specimen heater, 
Q»., will be referred to as the lower part of the speci- 
men, and the region above Q» as the upper part of the 
specimen. 


b. Description of Tests 


In the longitudinal heat flow method of measuring 
the thermal conductivity each datum point was com- 
puted by simultaneous solution of three tests: 

1. an “isothermal” test with no power input to the 
specimen heater (Q2) and with the temperature distri- 
bution on the guard adjusted to closely match that on 
the specimen. 

2. a “matched” gradient test with sufficient power 
input to the specimen heater to maintain the desired 
longitudinal temperature gradient in the specimen and 
with a matched temperature distribution on the 
guard. 

3. an “unmatched” gradient test with the power 
input to the specimen heater and the temperature at 
the center of the measuring span the same as in the 
“matched” gradient test, and with the temperature 
distribution on the guard parallel to that on the speci- 
men but 10 deg cooler. 


Matched Gradient Test 


The furnace temperature was raised by means of 
heater Q;. The power to the specimen heater, Qo», 
was adjusted to give a temperature gradient of 5 
deg/cm in the lower part of the specimen. Power to 
Q:; was automatically controlled, using a_propor- 
tional controller, to maintain a minimum tempera- 
ture drop across the neck. The Pt-10 percent Rh 
legs of thermocouples 9 and 12 were used-in con- 
junction with the necked-down portion of the speci- 
men as a differential thermocouple activating the 
proportional controller. The temperature drop across 
the neck never exceeded 0.1 deg. The specimen was 
maintained at the required mean temperature by 
thermostatting the power to the lower heater Q). 

The temperature distribution along the guard was 
forced to match that along the specimen by adjusting 
the controllers for heaters Qs, Q;. and Qs. Tempera- 
tures at corresponding locations on the specimen and 
the guard generally agreed to within 1 deg. 

After allowing time for the system to come to equi- 
librium, readings were taken of the thermocouple 
emfs and the voltage and current to the specimen 
heater, Q2. Normally these data were taken three 
times over a period of about 2 hr. The temperatures 
never drifted more than a few hundredths of a degree 
from one set of readings to the next and the three 
sets of data were averaged. When the drift between 
the first and second sets of readings was less than 
0.01 deg, the third set of readings ,was not taken. 
On completion of the last set of readings, the voltage 
drops between the inner and outer taps of each cur- 
rent lead were measured with the heater current 
flowing normally and then reversed. 


Unmatched Gradient Test 


Upon completion of the “matched” gradient test 
the controllers for the guard heaters Q,, Q;, and 
Q« were adjusted to lower the temperatures on the 
guard by 10 deg while maintaining the temperature 
distribution parallel to that on the specimen. With 
the guard at the lower temperature heat losses from 
the specimen to the surrounding insulation were sig- 
nificantly increased and the heat flow in the specimen 
reduced. As a result, the temperature gradient in the 
specimen and its mean temperature were decreased. 
The power to heater Q; was adjusted to restore the 
specimen to the initial mean temperature. The system 
was then allowed to equilibrate and the same data 
were taken as for the “‘matched” gradient test. 

isothermal Test 

For the “isothermal” test heater Q. was shut off 
and heater Q, adjusted until the specimen was approxi- 
mately isothermal. No adjustments had to be made to 
Q; as the controller automatically adjusted to maintain 
71: — T19=0. The guard heaters were likewise adjusted 
until the temperature distribution on the guard once 
again matched that on the specimen. When the sys- 
tem was in equilibrium, the data mentioned above were 


taken, with the exception of the power to the specimen 
heater which was shut off. 
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At 1100 °C an additional ‘“‘matched”’ gradient test was 
run after the three regular tests were completed. The 
extra data obtained in this test served as a check on 
drifts in thermocouple calibrations during the testing 
period. 

c. Testing Sequence 

All of the measurements described above were 
made at a number of temperatures. Tests were first 
run in air at 100 °C, then in argon at 100, 300, 500, 
700, 600, 400, and 200 °C in that order. One of the 
guard heaters would short out at about 750 °C and 
consequently the upper limit on the first run was 
700 °C. After completing that run in argon the system 
was evacuated and backfilled with helium (99.99% 
pure) to a pressure of about three-quarters of an atmos- 
phere. The helium, having a much higher thermal 
conductivity than argon, changed the effective thermal 
conductivity of the insulation surrounding the speci- 
men. Tests were run in helium at 200 and 400 °C to 
experimentally evaluate heat losses from the specimen 
to the insulation. The system was then opened up and 
the trouble with the guard heater corrected. The fur- 
nace was filled with fresh powder insulation, the pow- 
der being packed lightly around the neck to ensure 
that the necked-down region was uniformly filled with 
insulation. The system was refilled with argon as de- 
scribed at the beginning of this section and tests were 
run at 300, 700, 900, and 1100 °C. Heater Q; burned 
out while tests were being conducted at 1100 °C by 
the electrical method and those tests were incomplete. 
The thermocouples started drifting at 1100 °C and 
losing their calibration due to contamination. Conse- 
quently, it was decided to terminate the tests at that 
point. 


4.2. Calculation Procedures and Uncertainties 


For steady-state heat flow, the 
total heat flow, Q, through the specimen is given by 


one-dimensional 


(1) 


where A is the thermal conductivity, A is the cross- 
sectional area of the specimen, 7 is the temperature 
and z is the longitudinal coordinate. For moderate 
temperature ranges, the thermal conductivity of the 
specimen can be assumed to vary linearly with 
temperature; then (1) becomes 

Q=—doA {1+ BAT —T)} FZ, (2) 


where Avy is the thermal conductivity of the specimen 
at a reference temperature, 7, and Bp is its correspond- 
ing temperature coefficient. The difference between 
the heat flows in two tests is given by 


Q-Q’=- rel (S)- = | 


— ho T—Ty)(2)-(7' =o (4 |. (3) 


dz 


where quantities of one test are distinquished from 
those of the other by use of primes. If we define the 
reference temperature as 


_1T(dT/dz) —T' (dT /dz)' 
(dT/dz) — (dT/dz)' 





To (4) 
the second term on the right-hand side of (3) vanishes, 


and the thermal conductivity at the reference tempera- 
ture, 7, is given by 





A|(dT/dz) — (dT/dz)'J) ; 
The purpose of computing the thermal conductivity 
from data corresponding to two different powers was 
to correct for errors that did not depend on the power 
transmitted through the specimen. The most obvious 
errors of this type are thermocouple errors. In a 
simultaneous solution each thermocouple in effect 
measures a temperature difference so that errors in 
calibration of the thermocouples cancel out to first 
order. Further possible sources of error will become 
evident below. Determination of A» involved measure- 
ment of the cross-sectional area of the specimen, the 
total heat flow in the specimen and the longitudinal 
temperature gradient in the specimen for each of 
the two tests. These quantities were evaluated at the 
position of the middle thermocouple in the measuring 
span. 


a. Cross-Sectional Area 


The effective cross-sectional area of the specimen, 
after correction of the measured diameter for surface 
roughness, was determined to be 3.1331 cm? at 21 °C. 
The uncertainty * in this area was estimated to be less 
than 0.02 percent. The diameter at temperature ¢ 
(°C) was computed from that at 25 °C using the equation 


D,= D2;(0.99978 + 8.876 X 10-&t+ 1.311 x 10-%?). (6) 


This equation was derived from smoothed thermal ex- 
pansion data for platinum [19]. The cross-sectional 
area was then computed from the diameter. 


b. Heat Flow 


The total heat flow through the specimen at the 
position of the center thermocouple was calculated 
using the expression 


QO=P—qa— qo—Gn— Gi— Ge; (7) 


where P is the measured electrical power input to the 
specimen heater; 
da and q» are the heat losses along the two leads 
carrying current to the specimen heater; 
Gn is the heat loss across the necked-down portion 
of the specimen; 


* Uncertainties stated in this paper represent either (a) statistical uncertainties based on 
results of calibrations or (b) limits to errors conservatively estimated by the authors. 
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FIGURE 4. Circuit diagram for the specimen heater. 


gi is the heat loss into the insulation surrounding 

the specimen: and 

qe is the heat loss down the thermocouple wires and 

ceramic tubes next to the specimen. 
Each of the quantities in (7) is considered separately 
below. 

Power Input to Specimen Heater (P): The electrical 
circuit for the specimen heater is illustrated diagram- 
matically in figure 4. The power input to the heater 
was computed using 


P=nE,.14+ R,/R,\Es/Rs—nE,/R,). (8) 


where E,. is the voltage drop across the output of the 
volt box as measured with a potentiometer, n is the 
resistance ratio of the volt box, F, is the voltage drop 
across the standard resistor, R,, as measured with a 
potentiometer, X, is the total resistance of the potential 
leads and R,. is the total resistance of the volt box. 
The potentiometer, voltbox, and shunt box were 
each calibrated to 0.01 percent or better. The emf of 
the standard cell was known to 0.01 percent or better. 
The correction terms in (8) for the voltage drop in 
the potential leads and for the current through the 
voltbox were small (a few tenths of a percent) and 
uncertainties in these corrections could not have 
introduced more than 0.01 percent additional error 
in P. Thus the percentage uncertainty in the measured 
electrical power input was less than 0.05 percent. 
Heat Flow in Current Leads (qa. qu): The circuit 
diagram in figure 4 shows only one set of potential 
leads coming from the heater. In fact there were two 
sets as shown in figure 5, but only the inner set, i.e., 
the potential taps closer to the heater, was used in 
measuring the power input to the heater. The two sets 
of potential leads were required in measuring heat 
conduction along the current leads. The current leads 
were platinum and the potential leads platinum-10 
percent rhodium. Consequently. each current lead 
could be used along with its two potential leads as a 
differential thermocouple to measure the temperature 
drop, AT, between the potential taps. With current 
flowing to the heater the voltage drop measured be- 
tween 1-2 or 3-4 was the algebraic sum of the /R drop 
between adjacent potential taps (where / is the current 
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FIGURE 5 {rrangement of potential leads for the specimen heater. 


and K is the resistance of the current lead between 
the potential taps) and the Seebeck emf due to the 
temperature drop, AT, between the potential taps. 
Assuming no_ heat from the current leads. 
it can be shown [9, 20-26] that the heat conducted 
along the leads is given by 


losses 


Ap on 
da. b 3 AT - (9) 


I*R 
2 . 


where p is the electrical resistivity of the current 
lead and J its thermal conductivity. By taking measure- 
ments with the current flowing forward and reversed. 
R and AT were determined for each lead, and con- 
sequently ga and qp. 

The heat conducted along the leads, as determined 
using (9), was less than 0.05 percent of the heat flowing 
in the specimen. These corrections were sufficiently 
accurate that no errors larger than 0.02 percent are 
believed to have been introduced into the measured 
thermal conductivity values by uncertainties in qa 
and q». We will discuss below the possible magnitude 
of heat losses from the current leads into the surround- 
ing powder insulation. 


Heat Flow across Neck (q,): Referring to figure 1, 
it is seen that heat could be conducted across the 
necked-down region of the specimen by the neck itself 
and by the powder insulation surrounding it. The con- 
ductance of the powder was kz(b? — a”)/2/ where « is 
the thermal conductivity of the powder, 6 the radius 
of the specimen, a the radius of the neck, and 2/ the 
length of the neck. The conductance of the specimen 
between thermocouple positions 10 and 11 was AF 
where A is the thermal conductivity of the specimen 
and F is a geometric factor which was determined from 
ihe corresponding equation for electrical conductance, 
!=oaFV.where | is the current, o the electrical con- 
ductivity and V is the voltage drop between the thermo- 
couples. These data were available from measurements 
by the electrical method. The heat flow across the 
necked-down region was given by 


w(b2— a?) 


dn AF +k ry) AT, (10) 


where it is assumed that the temperature differential 
across the insulation was the same as that measured 
between the thermocouples. 

The correction for heat flow the necked- 
down region of the specimen was always less than 
0.05 percent of P. Errors in the temperature drop 
across the neck due to possible inhomogeneities in 
the thermocouple leads or to stray thermal emfs would, 
for the most part. cancel since they were common to 
Gn and q,, (i.e., corresponding to Q and Q’). At 300 
°C a series of tests were run in which the tempera- 
ture drop across the neck was held in turn at about 
—5. 0, and +5 deg, a range 20 times larger than that 
which occurred during normal measurements. The 
corresponding values for (g,—g}) were about +1 
percent, 0 percent, and —1 percent, respectively, of 


across 
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P. The three thermal conductivity values obtained, 
using (10) to effect the correction for heat flow across 
the neck, fell within a range of less than + 0.02 per- 
cent. For normal tests, in which the temperature across 
the necked-down region was maintained quite small, 
it is felt that any uncertainty in the measured thermal 
conductivity values due to this source was less than 
0.02 percent. 

Heat into the Insulation (qi): In order to 
determine gi, the heat exchange between the speci- 
men and the surrounding insulation, it was neces- 
sary to perform an extensive mathematical analysis. 
if the temperature distribution along the guard exactly 
matched that along the specimen there would have 
been no radial heat exchange between the specimen 
and the guard. However, there would still have been 
an exchange of heat between the specimen and the 
surrounding insulation in order to provide the longi- 
tudinal heat flow in, the insulation adjacent to the 
specimen. 


Loss 


The heat flow from an elemental length of the surface 
of the specimen was 


(ED) 


where a is the radius of the specimen, «x the thermal 
conductivity of the insulation, @ the temperature in 
the insulation relative to an arbitrary fixed tempera- 
ture, r the radial coordinate, and z the longitudinal 
coordinate. The net heat flowing across the surface 
r=a between z; and z» was 


gi(z1. 22) = 27a [ ms ( 2 dz, 


J2 Or /r=a 


(12) 


where k is, in general, temperature dependent. Let 
us define a new potential, €, that satisfies the relation 


ko VE= k(8) VO, (] 


3) 


where Ky = k(0). Integrating (13), 


"#4 


K(0)dé, (14) 


Koy 
where we have selected the integration constant so 
€=0 when 0=0. Writing (12) in terms of € we 

vet 
gilz1. 22)=KoD(z,. 22), (15) 


where 


(16) 


a 


z» a 
D(z, z22)= 27a [ (5) dz. 
0 2) gr r 


The factor D(z,, z2) was determined by analyzing 
the heat flow in the hollow cylinder of powder insula- 
tion between the specimen and the guard, using the 
measured temperature distributions along the speci- 


men assembly and along the guard cylinder as bound- 
ary conditions. Polynomial expressions relating tem- 
perature to longitudinal position were used to describe 
these temperature distributions in the regions between 
the heaters. In the intervening regions near the heaters, 
smoothing cubics [27] were used which provided conti- 
nuity of temperature and longitudinal temperature 
gradients. The evaluation of D(z,, z) is described more 
fully in appendix A. 

Since the heat flow in the specimen was to be evalu- 
ated at the position of the center thermocouple in 
the gradient region, all the heat loss to the insulation 
from the location of the heater down to the position of 
the center thermocouple had to be considered. In 
addition, heat losses in the region between the speci- 
men heater and the neck had to be considered since 
these had to be provided by the specimen heater. 
The neck, in effect, could be considered as the upper 
end of the specimen for purposes of this analysis, 
since any heat exchanges between the powder and the 
specimen above the neck did not affect that part of 
the specimen below the neck as long as zero tem- 
perature differential was maintained across the neck. 
Therefore, in evaluating qi, the limits of integration 
for D(z,. 2) were the position of the center thermo- 
couple (z;) and the position of the center of the neck 
(Z2). 

The correction (gi—qi') for heat exchange with the 
powder insulation was potentially a large source of 
error and considerable effort was expended to, first, 
keep this correction small and, second, evaluate it 
accurately. Evaluation of this correction, as seen from 
(15), required a knowledge of the integral, D(z:. zz), 
and of the thermal conductivity, Ko, of the insulation. 
_ Numerous factors could have adversely affected the 
determination of D(z, z). The use of logarithmic 
functions to define the radial temperature distribu- 
tion across the ends of the hollow cylinder of insula- 
tion was an approximation. However, it is easily 
shown that the potential distribution near the speci- 
men was not significantly affected by the boundary 
conditions at the remote ends of the extensions. 

In the mathematical analysis it was assumed that 
the temperatures on the inner surface of the guard 
were the same as the temperatures measured on the 
outer surface. The molybdenum guard had high 
thermal conductivity so that any radial temperature 
sradients in the guard would be small and the asso- 
ciated errors would tend to cancel on simultaneous 
solution of the gradient and isothermal tests. Angular 
variations in the temperature distribution on the guard 
could have arisen if the specimen and guard were not 
concentric or if the insulation between the specimen 
and the guard was not packed uniformly. Great care 
was taken to avoid both of these conditions. Any 
angular variations would have been approximately 
the same for two tests at the same mean temperature 
and so the associated errors would in large part cancel 
under simultaneous solution. Such would not be the 
case for errors arising from uncertainties in the longi- 
tudinal positions of the thermocouples since the 
temperature distribution along the guard cylinder in 
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the gradient test differed from that in the isothermal 
test. The necessary steps were taken to ensure that 
the longitudinal position of the guard was accurately 
known relative to that of the specimen, and the loca- 
tion of the thermocouple slits both on the guard and on 
the specimen assembly were measured accurately 
prior to installation of the thermocouples. 

The details of the actual temperature distributions 
in the heater regions where smoothing cubics [27] 
were used could conceivably have influenced D(z,. z2). 
Except for the regions of Q. and Q;, such effects 
should have been about the same in the gradient and 
in the isothermal tests, and hence would have canceled 
on simultaneous solution of these tests. Since the 
guard heaters were on the outside of a considerable 
thickness of high-conductivity metal, we feel that the 
temperature on the inside surface of the guard cylinder 
varied smoothly with position in the regions of the 
guard heaters. Thus, no significant errors were 
believed to be introduced by the use of a smoothing 
function in the region of Qs. 

R: W. Powell [28] has pointed out that a possible 
additional source of error, not specifically discussed 
by O’Hagan [9], is heat loss into the powder insulation 
from the external platinum jumpers which connected 
the elements of heater Q» (see sec. 3.3a and fig. 3). 
There were five such jumpers, each about 0.3 cm long, 
contained in aluminum oxide tubing to electrically 
insulate them from the specimen. Although the heat 
generated in the jumpers was only a small fraction of 
the total heat generation in the heater, these jumpers 
were heated, by conduction from the helical heater 
coils inside the specimen, to a temperature above that 
of the specimen. This would have resulted in a heat 
flow into the powder insulation surrounding the jump- 
ers. A portion of this heat would have flowed back 
into the specimen but, at least in principle, a net por- 
tion of the power input to the specimen heater could 
have been lost from the heater jumpers with a corre- 
sponding error in the measured thermal conductivity 
values. 

Since the temperature rise of the jumpers was de- 
pendent upon the power input to the specimen heater, 
the heat loss discussed in the previous paragraph was 
not eliminated, or even reduced, by the simultaneous 
solution of a matched gradient test and an isothermal 
test. The use of an unmatched gradient test also did 
not help in evaluating this source of heat loss. The 
smoothing functions used were quite adequate for the 
isothermal tests but did not truly represent the temper- 
ature distributions in the region of Q» for the gradient 
tests. In view of the importance of this potential source 
of error, heat loss from the jumpers is considered 
further in appendix B, where a mathematical analysis 
is used to approximately evaluate this source of error. 
This analysis indicates that errors in the measured 
thermal conductivity values due to heat loss from the 
jumpers on heater Q». were less than 0.2 percent at 
100 °C and less than 0.5 percent at 900 °C. 

A calculation (see appendix I of O’Hagan [9]), based 
on the degree of tempering provided by the swaged 


elements to which the heater current leads were 


attached, indicated that the current leads were only 
1 to 2 deg C hotter than the adjacent specimen ma- 
terial. By comparison to the discussion of heat loss 
from the jumpers (see appendix B), which may have 
been over 100 deg C hotter than the specimen, we see 
that any errors due to heat loss from the heater leads 
into the powder insulation were negligible. 

In analyzing the data the “isothermal” test was 
combined with the “matched gradient” test and with 
the “unmatched gradient” test to give two equations 
of the form (5). Inspection of eqs (5), (7), and (15) 
shows that the value obtained for Ay depends in a 
linear manner on the value assumed for ko. The ef- 
fective thermal conductivity of the insulation § sur- 
rounding the specimen depends on the density of the 
powder and the pressure and type of gas present, and 
is best determined under experimental conditions. 
This was done by simultaneous solution of two equa- 
tions of the form (5) which yielded values both for the 
thermal conductivity of the specimen, Ao, and the 
thermal conductivity of the insulation, ko. The thermal 
conductivity values obtained for the aluminum oxide 
insulation, in argon and in helium, were given by 
O’ Hagan [9]. 

If there was a significant heat exchange between the 
specimen and the insulation that was not being 
adequately corrected for, one would expect a system- 
atic difference between the values obtained for the 
thermal conductivity of the specimen in helium and 
those obtained in argon, due to the large difference 
between the thermal conductivity of the powder in 
the different atmospheres. In fact, however, the values 
measured in helium fell within the scatter band of 
those measured in argon indicating that any uncor- 
rected heat exchanges were certainly less than 0.2 
percent, the width of the scatter band (see sec. 4.3). 

For all of the tests taken the correction |qj—q; | 
was less than 0.1 percent of P. It is felt that D(z, z) 
and ky were each known to better than 10 percent and 
hence the uncertainty in (g;—q;) less than 0.02 per- 
cent of (P—P’). However, in view of all the factors 
which conceivably could have influenced this correc- 
tion, an uncertainty of 0.1 percent is assigned to the 
measured thermal conductivity values due to possible 
errors in (qi—qi). To this must be added the un- 
certainty due to heat loss from the jumpers on heater 
Qs». 


Heat Loss along Thermocouple Wires and Insulators 
(qe): The heat loss, q-, along the thermocouples and 
ceramic insulators next to the specimen was computed 
from the expression 


~ 


(17) 


- dT 


ail 2 Ci ( dz 


where C; is the longitudinal thermal conductance of 
the ith wire and its insulator n is the total number of 
wires crossing the plane where the thermal conduc- 
tivity was evaluated, and d7/dz is the temperature 
gradient at that plane. Each C; was computed from 
the thermal conductivities and dimensions of the wire 
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and insulator. This correction was rather large, falling 
from about 0.8 percent at 100 °C to 0.3 percent at 
1100 °C, the falloff being due to the rapidly decreasing 
thermal conductivity of the ceramic tubing. Although 
the thermal conductivity of the thermocouple wires 
was known fairly accurately (5%) the thermal conduc- 
tivity of the ceramic tubing was known only approxi- 
mately (15%). Furthermore, the cross-sectional areas 
of the wires and tubing were not accurately known 
(5%). Thus the total conductance of the wires and the 
ceramic tubing was only known to about 25 percent. 
The corresponding uncertainty in the measured 
thermal conductivity values was 0.2 percent at 100 °C 
and 0.1 percent at 900 °C. 

Departure from Steady-State: The ratio of heat 
absorbed (or released) to that conducted in the speci- 
men is given approximately by 


wcAL(dT/dt) L (dT/dt) 13 

A(dT]dz)  D(dT/dz)’ ™ 
where w is density, c is specific heat, A is area, L is 
total length of specimen below the necked-down 
region, d7/dt is time rate of temperature change, J is 
thermal conductivity, d7/dz is temperature gradient, 
and D=)/wc is thermal diffusivity. Temperatures in 
the system did not drift at a rate greater than 0.03 
deg/hr (i.e., 10~° deg/sec): the length, L, was about 
14 cm: the temperature gradient was 5 deg/cm:; and 
the thermal diffusivity of platinum in the temperature 
range 0-1100 °C is always greater than 0.2 cm?/sec 
[12]. Hence the ratio of heat absorbed (or released) 
to that conducted was less than + 0.02 percent. No 
correction was made for departure from steady-state 
conditions. 

Total Uncertainty in Heat Flow: Using standard 
propagation of error formulas, the estimated uncer- 
tainty in (Q—Q') was 0.4 percent at 100 °C and 0.7 
percent at 1100 °C. 


c. Temperature Gradient 


The temperature gradient in the specimen was 
computed from the measured temperatures at the 
five thermocouple positions in the gradient region. 
The separations between thermocouple grooves at 
room temperature were accurately measured before 
the thermocouples were installed: the separation at 
elevated temperature was computed using (6). Since 
temperature gradients in the specimen were rather 
small (less than 5 deg/cm) it was essential that the 
conversion of thermocouple emfs to temperature not 
introduce any additional uncertainties. The equation 


E=15.83952 eo) — 9.18328 (sean) 


i \ t 4 
Ma 572 (T5595) — 9275: (—) 
+ 7.30572 | 7999) — 1.92793 | F999 


— 2.50480 (1.0—exp | — 4.18312(t/1000)]), 


(19) 


TABLE 1. 


where ¢ is temperature (°C) and E is emf (mV), was 
found to fit to within 1 «V the calibration data for 
the platinum versus platinum—10 percent rhodium 
thermocouple wire from which the specimen thermo- 
couples were fabricated. This equation was used for 
converting the thermocouple voltages to temperatures. 

The temperature gradient was computed by evaluat- 
ing the slope, at the center thermocouple location, of 
the quadratic equation of least-squares fit to the five 
temperatures and thermocouple positions. The uncer- 
tainty in the temperature gradient due to uncertainties 
in effective thermocouple positions is estimated to have 
been less than 0.2 percent. Errors in reading thermo- 
couple emfs did not introduce an uncertainty of more 
than 0.05 percent for the temperature gradients used. 
Errors due to heat conduction along thermocouple 
leads should have been negligible and were certainly 
less than 0.05 percent (see Appendix G of O'Hagan 
[9]). Due to the use of a simultaneous solution of a 
gradient and an isothermal test, errors in the measured 
temperature gradient due to variations between in- 
dividual thermocouples are estimated to have been less 
than 0.05 percent. It is estimated that the conversion 
of thermocouple emfs to temperatures introduced 
errors of less than 0.2 percent in the temperature 
gradients. 

The estimated overall uncertainty in the temperature 
gradient is estimated to have been 0.3 percent. 


d. Mean Temperature 


In addition to the uncertainties discussed above in 
the area, heat flow, and temperature gradient, there 
is an uncertainty in the temperatures to which the 
thermal conductivity values correspond. For a 0.5 deg 
uncertainty in temperature, the associated uncertainty 
in thermal conductivity is less than 0.001 percent at 
100 °C and less than 0.02 percent at 900 °C. 


4.3. Results 


The experimental values obtained for the thermal 
conductivity of our platinum specimen by the longi- 
tudinal heat flow method are given in table 1. The 


Experimental values for the thermal conductivity of 
platinum as measured using the longitudinal heat flow method 


The values given are corrected for thermal expansion. 


Mean Thermal 
temperature | conductivity 
( 


W/cm deg 


Run Atmosphere 


Air 99.6 
Argon 99.8 
Argon 301.0 
Argon 501.5 
Argon 701.2 
Argon 601.3 
Argon 100.2 
Argon 201.7 
Helium 199.6 
Helium 400.6 
Argon 300.0 
Argon 701.3 
Argon 900.0 
Argon 1100.5 
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TABLE 2. Typical set of data from measurements of the thermal 
conductivity of platinum by the longitudinal heat flow method 


The data correspond to Test No. 6 in table 1. 


Matched 


gradient 


Unmatched | Isothermal 


gradient 


Temperature distribution along the 

specimen... ( P | 80.4 579.0 
990.9 588.9 
601.3 598.8 
611.7 608.9 
622.1 619.0 


600.7 
600.6 
600.3 
600.1 
600.0 
Temperature gradient at the location of 


thermocouple 7; deg/cm ~ 5.18 4.96 0.09 





Power generated in specimen 
heater W 
Heat flow across the necked-down 
region W 0.003 0.008 
Heat flow along the current leads W | 004 O17 
003 006 
Heat flow to the insulation W ] 006 542 


12.853 


12.818 


Heat flow along the thermocouple 
wires and insulators W ‘ 051 | 
i 


, 








~ 


thermal conductivity values given there have been 
corrected for thermal expansion of the specimen. A 
typical set of data, including the various correction 
terms, is given in table 2. 

A cubic equation of least-squares fit was found to 
fit the data (corrected for expansion) from the longi- 
tudinal heat flow method with the residuals having a 
standard deviation of 0.08 percent. This was signifi- 
cantly less than the standard deviation of the residuals 
from a parabola. The equation, valid over the range 
100 to 900 °C, is 


A=0.713 4+ 0.683 X 10-°t + 0.173 X 10- ®t? 


= 0.313 10788. 


(20) 
where ¢ is temperature in °C and J is in W/em deg. 
Departures of the data points from (20) are plotted in 
figure 6. With the exception of the point at 1100 °C 
all the data points, including the two obtained in 
helium, fall within + 0.1 percent of the curve. There 
were no significant differences between the values 
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ductivity 
method. 


Percentage departures from eq (20) of the thermal con- 
data points obtained by the longitudinal heat flou 


TABLE 3. Summary of individual uncertainties contributing to 
the overall uncertainty in the thermal conductivity results obtained 
using the longitudinal heat flow method 


Uncertainty, % 


Source of uncertainty 


1 100 ¢ 900 °¢ 


Cross-sectional area 
Heat flow 
Power input to specimen 
heater 
Heat flow in current leads 
Heat flow across neck 
Heat loss into the insulation | 
Heat loss along thermocot 
ple wires and insulators 
Departure from steady-state 
Temperature gradient 
Mean temperature 





Combined 


Each combined uncertainty was obtained by t the sum of the 


squares of the individual uncertainties 


obtained on heating and cooling in the first run and 
1100 °C the thermo- 
drifting giving rise to significant 
uncertainties in temperature measurement. After 
completing the “matched” gradient, “unmatched” 
gradient, and “isothermal” tests at 1100 °C a second 
“matched” gradient test was run. Two values for the 
thermal conductivity were obtained by simultaneous 
solution of each of the two ““matched” gradient tests 
with the “isothermal” test. The first value (0.820) 
was below the value predicted by (20) and the second 
value (0.905) was above it. Assuming that the tempera- 
ture drifts were linear with time, interpolation to the 
time of the “isothermal” test gave a value of 0.866 
for the thermal conductivity of the specimen at 1100 
°C. This value is 0.6 percent above the value given by 
extrapolation of (20) to 1100 °C. The 1100 °C point, 
due to the larger uncertainty associated with its value, 
was not used in deriving (20). 

The thermal conductivity. values corresponding to 
(20) are believed to be uncertain by not more than | 
percent over the temperature range 100 to 900 °C. 
The estimated uncertainties arising from the various 
known sources of error are summarized in table 3: 
these uncertainies were discussed in section 4.2. 


couples started 


5. Electrical Method 


5.1. Experimenta! Procedure 


Measurements of thermal conductivity by the 
method in which the sample was heated directly by 
passage of an electric current were made following 
each measurement by the longitudinal heat flow 
method. Thus the preliminary procedures and testing 
sequence described in section 4.1 also apply to these 
measurements. Power to heaters Q; and Qs was con- 
trolled to maintain 7,9 and 7); constant with |7\; — Tio! 
less than 0.2 deg. Temperatures along the guard in 
the region opposite the neck were maintained at the 
same value as 79, Data were taken with currents (both 
normally and reversed) of approximately 10, 58, 82, 
and 100 A de flowing through the specimen. These 
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current levels correspond approximately to even 
increments of power generation. Readings were taken 
of thermocouples 9, 10, 11, and 12 (see fig. 1) on the 
specimen and of the guard thermocouples between 
heaters Q; and Qs. The voltage drops between the 
platinum legs of thermocouples 10 and 11 and also 
between 9 and 12 were measured with the current 
being measured simultaneously on another poten- 
tiometer. 

The temperatures never drifted more than a few 
hundredths of a degree during a test (approximately 
three-quarters of an hour), and although the current 
showed drifts of 0.05 percent the measured resistances 
remained constant to 0.01 percent with current and 
voltage being read simultaneously. 

5.2. Calculation Procedures and Uncertainties 

Consider a conductor of arbitrary geometry which is 
perfectly insulated, both thermally and electrically, 
except for the ends where electric current enters and 
leaves the conductor. Joule heat will be generated in 
the conductor by the electric current. Let 7, be 
the maximum temperature in the conductor. Let V be 
the voltage drop between two surfaces, one on either 
side of the surface of maximum temperature, which 
are both at some lower temperature, 7». It is shown 
in appendix C, where earlier work is referenced, 
that the following relation holds: 


~ \pdT | " Apdo, (21) 


0 


where 0=T—T . 0n=Tm—To, 
‘ductivity of the conductor, 
resistivity of the conductor. 

Equation (21) shows that the maximum temperature 
rise, @,. in an electrically heated conductor with 
negligible lateral heat losses and with both its — ds 
held at the same temperature is a function only oi the 
voltage drop | and of the thermal and 
electrical conductivities of the material, and is in- 
dependent of the geometry of the conductor provided 
the geometry is suc h as to make lateral heat losses 
negligible. If the maximum temperature rise can be 
measured as a function of the applied potential, Ap 
is readily determined. 

The maximum temperature rise, @,,, can be meas- 
ured indirectly using the specimen as its own resist- 
ance thermometer. This method was employed by 
Holm and Stérmer [13] and Cutler et al. [15]. and 
requires that the temperature coefficient of resistance 
of the material be known over the temperature range 
of interest and that it be large enough to yield sufficient 
sensitivity. Platinum satisfies both requirements very 
well and the resistance method of measuring the 
maximum temperature rise was adopted in the present 
experiments. 

It is shown in appendix C that the following relation 
holds for a perfectly insulated conductor: 


\ is the thermal con- 
and p is the electrical 


across it, 


oe: 99 
R FJo F = 


Re 1 dol nO 


where 


(23) 


4 "2 
{2{ Apae 


and R is the electrical resistance between the two 
surfaces, one on either side of the surface of maximum 
temperature, which are both at the same temperature, 
T. The quantities Ro and po are the values R and p 
would have if the conductor were isothermal at some 
reference temperature, 7). The quantity yo is the 
coefficient of linear thermal expansion at 7». 
Assuming that A and p can both be represented by 
linear functions of temperature—a valid assumption 
over small temperature intervals—we can write: 


P= po(1+ a8). (24) 


nN ho( 1 t Bo), (25) 


Ap = Aopo( 1 + no), (26) 
where Q@ is the temperature coefficient of resistance 
evaluated at @=0,. By is the temperature coefficient of 
thermal conductivity evaluated at 0=0 and no=ao 
+ Bo. The term in a@oBo0? in (26) has been neglected. 
Substituting the above expression for Ap into (23) 
and integrating we get 


F (0, Om) 
») 


| 2aups 4m (1 b= A» | -O(1 


If the integration is evaluated for T 
have, from (27) and (21), 


1/2 


2a (1+) 


/ 
© J 


and 


6). 


SAoPo0Am | | (29) 


Substituting F from (27) into (22) and performing the 


indicated integration we have for @=0, 


Ro 


»] 
— arctan G, 
( 


where 


No ; 
= a): (31) 


G2 =2noOm (1 


The coefficients a» and By are the “true” coefficients 
in that they correspond to electrical resistivity and 
thermal conductivity values which have been corrected 
for thermal expansion. If one considers “‘apparent”™ 
coefhicients a) and Bo which correspond to “apparent” 
thermal conductivity values which not 


have been 
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corrected for thermal expansion, then it is apparent 
that ap=ao+ yo. Bo=Bo— Yo. No= No. and (30) can 
, be written as 


. > ; 
Ro _ Bo , a : arctan G, (32) 
R Tho To G 
where, reiterating, a>, Bo. and no=aot+ By are ap- 
parent coefficients corresponding to electrical resis- 
tivity and thermal conductivity values which have 
not been corrected for thermal expansion. Thus 
the maximum temperature rise in the conductor 
can be computed without knowing the coefficient of 
expansion. 

Equation (30) gives @,, implicitly in terms of known 
measured parameters, and its value can be determined 
_by iteration. In practice it would generally be neces- 
sary to measure K at at least two different levels of 
heating current and then extrapolate to zero current 
to obtain Ry. Once 0,» is obtained using (30), its value 
can be used in conjunction with (29) to obtain the 
thermal conductivity. In the present investigation, 
R and V were measured at four or five different current 
levels. The data obtained were analyzed as described 
below. 


An adjusted voltage, V*, was defined as 


In the limit of small ao9». (30) reduces to 


es * Tee > Me (oO ) 
aoRo o 


An adjusted resistance, R*, was defined for the case 
where G0, is not sufficiently small for (34) to be valid: 


(35) 


R*=Ro L-4 = 200m ). 
re) j 


In the limiting case of small a9. (35) reduces to (34) 
and R* to R. 


Substituting (33) and (35) into (29), we obtain 


: : R ae ra 
y*= 12depo— : . 
aoRo 


Differentiating (36) and rearranging, we obtain 


= aoKo dV *2 


a (37) 

12po dR 
The thermal conductivity values in the electrical 
method were obtained from the slope of the V*? 
versus R* curve, using (37). From (33) and (35) we see 
that 


(38) 


In the analysis given in appendix C and also that 
given just above, we explicitly assumed that there 
was no loss or gain of either electric current or heat 
across the surface of the conductor. In the present 
investigation a necked-down sample was employed, 
the neck and the specimen as a whole being sur- 
rounded by an insulating powder. Since this insulating 
powder had an electrical conductivity many orders 
of magnitude lower than that of the specimen material, 
the assumption of no flow of electric current across 
the boundaries was completely valid. The powder 
surrounding the neck prevented heat loss by convec- 
tion, and radiation through the powder was negligible. 
However, the powder conducted heat away from the 
neck and it was necessary to analyze this heat loss 
and develop an appropriate correction for it. Since 
most of the temperature rise between thermocouple 
positions 10 and 11 was in the neck itself, only heat 
losses from the neck were considered and other heat 
exchanges were neglected. O'Hagan [9] has shown that 
in the presence of small heat losses from the neck, 
(37) is replaced by 

1 ao db 


ho=— s 
| T & 12po dR c 


(39) 


where 


(40) 


Ky being the thermal conductivity of the powder 
insulation at the reference temperature, To, and Ao the 
thermal conductivity of the specimen at Ty. The geo- 
metrical factor, Q. is given by 


Q) 


=|) 2!) x " 
Ko(aal2l) R2" ot) 


Ta 


where 2a is the diameter of the neck and 2/ is the 
length of the neck, and Ky and KA, are the modified 
Bessel functions of second kind and order zero and 
one, respectively. 

The calculation procedure used was as follows: 
At any given nominal temperature, 
were made at n current levels. An approximate 
value for Ro computed from the resistance 
corresponding to the lowest current level, for which 
Joule heating was minimal and the neck was nearly 
isothermal. For each experimentally determined re- 
sistance a corresponding approximate value of Om 
was computed using (34). Values of aj) were computed 
from the data obtained in the resistivity measure- 
ments (see sec. 5.3a). Using this value of 6» as the 
first trial value, Newton-Raphson iteration was used 
to compute the value of 0» which satisfied (30).4 
Values for R* and V** were then computed from (35) 
and (33) respectively. Since it was impractical to 
hold Ty (i.e. the average temperature of thermo- 


measurements 


was 


‘Values for B, computed from the thermal conductivity values obtained by Watson and 
Flynn [32] on the same specimen were used 
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couple locations 10 and 11) at exactly the same value 
for different current settings, it was necessary to 
adjust slightly the V*" and R* values to correspond to 


a common reference temperature. A straight line of 


least-squares fit through these values of R* and V*? 
gave Ry as the ordinate intercept and dV*2/dR* as 
the slope. This value for Ry was a more accurate value 
than the initial value used, as it was based on an 
extrapolation to /=0. The calculation was then re- 
peated using this value of Ry. This in turn led to 
improved values for Ry and dV*?/dR*. The iteration 
converged in a few passes and the final values for Ro 
and dlV*?/dR* were used in (37) to compute the ap- 
parent thermal conductivity. Using this apparent 
value for Ao, and eqs (40) and (41), the correction 
factor for heat losses from the neck was computed 
and used in (39) to give the final value for the thermal 
conductivity. 


a. Electrical Resistivity 
The 


from 


electrical resistivity values were computed 


Po PiceRo ( Ro hices ( 42) 


where pice is the independently measured ice-point 
resistivity, Ro is the resistance of the necked-down 
region at zero current, and (Ro) je is the value of Ro 
extrapolated to the ice point. These values were then 
corrected for thermal expansion. The uncertainty 
iN Pice Was not more than 0.1 percent. The uncer- 
tainties in the measured values of Ro did not exceed 
0.05 percent. However, an additional uncertainty 
must be assigned to the values of Ry due to the uncer- 
tainty in the temperature to which these values cor- 
respond. Quite conservatively, temperatures were 
known to within 0.5 deg, corresponding to 0.14 percent 
of Ro at 100 °C and 0.04 percent at 900 °C. The un- 
certainty in (Ro)iee. Which was based on an extrapola- 
tion, was estimated to be less than 0.25 percent. At 
lower temperatures the uncertainties in Ro were 
correlated with the uncertainty in (Ro). so that the 
uncertainties in Ro/(Ro)ico were lower than the un- 
certainties in either Ro or (Ro)ice at lower tempera- 
tures but not necessarily at higher temperatures. The 
overall uncertainty in electrical resistivity was esti- 
mated to be less than 0.1 percent at 0 °C and less than 
0.4 percent at 900 °C. 

The thermal conductivity was computed, using (39), 
from the ratio Ro/po. Since po was taken from a 
smoothed curve, this quantity was further uncertain 
by the scatter in the individual Ro data points around 
a smooth curve, or about an additional 0.05 percent 
(see sec. 5.3a). The overall uncertainty in Ro/po was 
probably less than 0.4 percent. 


b. Temperature Coefficient of Resistance 
The uncertainty in the sensitivity, ao, of the speci- 
men as a resistance thermometer was essentially the 
uncertainty in the sensitivity of the thermocouples 
used plus a small uncertainty in the resistance meas- 
urements and the thermocouple emf measurements. 


That this is so is seen by considering that Ro was 
measured as a function, say Ro=/f(E), of the thermo- 
couple emf, E. Then 
1 dRo 1 df(E) dE _ 

a-D . 


ate _ | GOEL = 
Ro dT f(E) dE dT 


2(E) a 


(43) 


The uncertainty in the sensitivity, dE/dT, of the 
thermocouples is estimated to have been less than 0.2 
percent while the uncertainty in the resistance and 
voltage measurements was less than 0.05 percent. 
Thus ap is believed to have been known with an un- 
certainty of less than 0.25 percent. 


c. Slope of |’** Versus R* Curve 


The voltage drop across the neck and the resistance 
of the neck were each measured with an uncertainty 
of about 0.02 percent. The change in the resistance of 
the neck for the different current levels was small 
compared to the resistance itself, so that any error in 
measuring resistance would be greatly magnified in 
computing the rate of change of resistance as a func- 
tion of the voltage drop. For all of the data taken, the 
departures of the R* values from the least-squares 
straight line fitted to the R* and V*? values were less 
than 0.01 percent and for a majority of the tests they 
were less than 0.005 percent. However, since (R — Ro) 
<R, the small scatter in R was highly magnified in cal- 
culating dV*?/dR*. What departures did exist tended 
to be systematic rather than random and tended to 
indicate that the plot of R* versus V** was very 
slightly concave downward rather than linear, as 
assumed. Such an effect could possibly have been due 
to neglecting higher order terms -in the temperature 
dependence of the thermal conductivity and the elec- 
trical resistivity. The uncertainty in dV*?/dR* is 
estimated to have been less than 1.5 percent. 


d. Heat Loss Correction 

The uncertainty in the heat loss correction, C, may 
have been as large as 50 percent of C. For the tests in 
argon this corresponds to an uncertainty in the meas- 
ured values for the thermal conductivity of the speci- 
men of 0.1 percent at 100 °C and 0.6 percent at 900 °C. 
For the tests in helium the corresponding uncertainties 
were 0.5 percent at 200 °C and 0.7 percent at 400 °C. 


e. Departures From Theory 


In the derivation of the mathematical expressions 
used to compute thermal conductivity in the electrical 
method, there were a number of explicit or implicit 
assumptions made which could lead to erroneous re- 
sults if these assumptions were not valid. 

It is shown in appendix C that the Thomson effect 
cancels out to first order provided the temperatures 
at the two potential taps (used to measure V) are 
approximately the same. The requirement that the 
potential taps be at the same temperature is also 
necessary for the Seebeck emf and Fermi energy 
terms in (C-50) to drop out. In all of the measure- 
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ments, the temperatures at the inner potential taps 
(thermocouple positions 10 and 11) agreed within 0.2 
deg or better, as compared with values of 6,, (for 100 
A current) ranging from 35 deg at 100 °C to 100 deg 
at the highest temperatures. Neglecting the Seebeck 
term in (C-50) for 
uncertainty in 
values of 


T, — T,|<0.2 deg corresponds to an 
the measured thermal conductivity 
than 0.05 percent. Neglecting the 
Fermi energy term in (C-50) also involves an uncer- 
tainty of less than 0.05 percent (see O'Hagan [9, 

127)). 

It is explicitly assumed in appendix C that the con- 
ductor is homogeneous and isotropic: this should be a 
valid assumption for platinum of the purity used in the 
present investigation. 

As stated in appendix C, the electron current is con- 
sidered to be the only mass current. In principle it 
would be possible for platinum, or any impurities 
which might be present, to migrate under the in- 
fluence of the electric potential gradient (electromigra- 
tion) or under the influence of the temperature gradient 
(Soret effect). If such mass motion existed there would 
be an associated entropy flow and equation (C—4) and 
all following equations in appendix C would have to be 
modified to include a term involving the mass current 
density and _ its entropy transport. A 
rigorous analysis including mass migration would be 
quite,complex and was felt to be beyond the scope of 
the present investigation. However, a qualitative dis- 
cussion of the probable effect of such mass motion if it 
were to occur is given below. 


less 


associated 


In electromigration, the energy flow is proportional 
to the electric field. Similarly, the energy flow due to 
the Thomson effect is essentially proportional to the 
electric field. Thus if one were to introduce an electro- 
migration term into (C-4) and rigorously go through 
the analysis, the equivalent expression to (C-31) 
would include an electromigration term of a form 
analogous to (i.e., proportional to the electric field) 
the term in (C-31) involving the Thomson coefficient. 
Provided this term were small compared with the first 
two termis in (C—35) it would, to first order, cancel in a 
similar manner to the the Thomson term can- 
celed, provided the temperature at the two potential 
taps were the same and the medium was homogeneous 
and isotropic. Physically what happens is that the 
energy transport due to electromigration adds to the 
energy transport by conduction on one side of the 
surface of maximum temperature and subtracts on 
the other side, with no net effect on the maximum 
temperature rise in the conductor or on the voltage 
drop between the potential taps, at least to first order. 
To reemphasize, the temperatures at the taps must be 
the same for this cancellation to occur and also the 
medium must be homogeneous and isotropic. 


Way 


In the Soret effect, the energy transport is propor- 
tional to the temperature gradient. Thus the ratio of 
energy transport due to the Soret effect to that due to 
heat conduction by other mechanisms is independent 
of the temperature gradient. The Soret effect is a dif- 
fusion process, similar to heat conduction. For a homo- 
geneous medium, therefore, the Soret effect simply 


behaves as an additional mechanism for “heat conduc- 
tion” and thus is a legitimate augmentation (or deple- 
tion, depending on the sign of the heat of transport) 
of the thermal conductivity and should be included in 
the total thermal conductivity value. 

For the conditions of the present experiments, there 
would be no errors involved due to electromigration or 
Soret effect provided the medium remained homo- 
geneous and isotropic. Both of these effects can change 
the distribution of impurities and vacancies in a solid. 
During the experiments the current was in one direc- 
tion through the specimen about one-half of the time 
and in the other direction about one-half of the time. 
Thus, it is doubtful if there was significant inhomo- 
genity introduced due to electromigration effects, 
even if they were occurring. The Soret effect, if large 
enough, could cause a redistribution of impurities in 
the specimen. 

For platinum of the purity used, the thermal con- 
ductivity at high temperatures would not be expected 
to be significantly affected by changes in the impurity 
and vacancy distribution. All of the data at low tem- 
peratures, where impurity concentration could affect 
the thermal conductivity, were taken before the speci- 
men was heated to temperatures where significant 
mass migration was likely to occur. In view of the 
above discussion it is felt that neither electromigra- 
tion nor Soret effect had any adverse effect on the 
results. 


5.3. Results 


a. Electrical Resistivity 
The ice-point resistivity determined on the 
platinum bar before the neck was machined in it (see 
section 2). The value obtained was 


was 


(Po ice = 9.847 + 0.010420 em, (44) 
C dimensions. 

In conjunction with the electrical method of meas- 
uring thermal conductivity, the resistance of the 
necked-down region of the specimen at each tem- 
perature was determined at a number of current levels. 
The resistance, Ro. at temperature 7. was evaluated 
by extrapolation to zero current as described in section 
5.2. In the first series of tests, or the first run as it 
will be referred to, measurements were made in the 
following order: 100, 300, 500, 700, 600, 400, and 
200 °C in argon: and 200 and 400 °C in helium. The 
apparatus was then opened up for repair of the guard 
heater as described in section 4.lc. In the second 
series of tests, or the second run as it will be referred 
to, measurements were made at 300, 700, 900, and 
1100 °C in argon, in that order. 

The values for the resistance of the neck obtained 
during the second run were about 1.5 percent higher 
than the values obtained at corresponding tempera- 
tures in the first run. This increase in resistance was 
observed both for measurements at the inner potential 
taps (10, 11) and for measurements at the outer poten- 
tial taps (9, 12), indicating that the increase in values 


corrected to 0 
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were due either to contamination or to a change in the 
geometrical factor of the necked-down region rather 
than to, for example, a short circuit which changed 
the effective position of one of the potential taps. 
Since data taken in the first run at lower temperatures 
after the sample had been heated to 700 °C agreed 
quite closely with results obtained before the sample 
was heated to 700 °C, it is felt that there was no 
evidence indicating that any chemical contamination 
of the necked-down region occurred during the first 
run. Furthermore, 700 °C is rather too low a tempera- 
ture to expect any significant contamination of plat- 
inum in a relatively clean environment. Exactly what 
happened is not known but it is believed that the 
necked-down region suffered some slight geometrical 
change between the first and second runs, perhaps 
while the guard heater was being repaired. In deter- 
mining the electrical resistivity, data from the first 
run had to be treated separately in view of this change 
in the geometric factor. 

A quadratic equation was fitted, by the method of 
least squares, to the measured values of Ro from the 
tests in air, argon, and helium from the first run. 
The equation, extrapolated to 0 °C and normalized 
to the ice-point resistivity, (44), gave 


p= 9.847(1 + 0.3963 x 10>-7t— 0.5389 x 10> &t?) (45) 
as the electrical resistivity (corrected for thermal 
expansion) of the neck from 0 to 700 °C, t being in 
°C and p in wQ cm. The standard deviation of the 
residuals from this equation was 0.03 percent. Depar- 
tures of the data points from this equation are plotted 
in figure 7 where the overall scatter is seen to be 
+ 0.05 percent. All but three of the points fall within 
0.025 percent of the curve and there is no significant 
difference between the values obtained on heating 
and cooling. 

A second quadratic equation was fitted to the data 
points (corrected for thermal expansion) at 300, 700, 
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FIGURE Percentage departures from eq (45) of the electrical 
resistivity data points. 
The values obtained during the second run h 
in the text 


ive been normalized at 700 °C as discussed 


and 900 °C obtained in the second run. The equation, 
normalized to the value given by (45) at 700 °C, gave . 


49 & 


p= 9.767(1 + 0.4033 x 10-?t — 0.5802 x 10> §t?), (46) 
t being in °C and p in pO cm. While the geometric 
factor of the neck changed between the first and 
second runs there is no reason to believe the resistivity 
would have changed. Any strains introduced when the 
geometric factor changed would have annealed out 
at 700 °C. Thus normalizing the data from the second 
run to that from the first run at 700 °C is felt to be 


justified. The departures of the normalized data points 


at 300, 700, 900, and 1100 °C from the values given by 
eq (45) are indicated by triangles in figure 7. The 1100 
°C. data point was not used in deriving (46). The 
thermocouples started to drift at 1100 °C due to 
contamination so that greater uncertainty had to be 
assigned to data obtained at that temperature. The 
1100 °C data point agrees with the extrapolated value 
given by equation (46) to within 0.02 percent. The 
deviations from the quadratic equation (45) at the 
higher temperatures as shown in figure 7 are in line 
with what one would expect. Recent measurements 
at NBS [29] in conjunction with work toward extending 
platinum resistance thermometry to the gold point 
showed that the measured resistance of a certain high 
purity platinum resistance thermometer at the gold 
point (1063 °C) was about 0.07 percent below the value 
obtained by extrapolating the Callendar equation [5] 
for that thermometer from 630.5 °C to the gold point 
(1063 °C). Laubitz and van der Meer [30], in measuring 
the electrical resistivity of high purity (99.999% pure) 
platinum, found that between 800 and 1200 °C their 
experimental results fell consistently below the values 
extrapolated from an equation similar to (45) by an 
average amount of 0.10 percent. 

As stated in section 5.2a, the uncertainty in the 
electrical resistivity values is estimated to be less 
than 0.1 percent at 0 °C and less than 0.4 percent at 
900 °C. 


b. Thermal Conductivity 


The experimental values obtained for the thermal 
conductivity of the platinum sample by the electrical 
method are given in table 4. The values for po and ap 
used in computing these thermal conductivity values 
were computed using (45) in the first run and (46) in 
the second run. The thermal conductivity values given 
in table 4 have been corrected for thermal expansion. 
A typical set of data is given in table 5. 

The following equation, obtained by the method of 
least squares, was found to fit the thermal conductivity 
data (corrected for thermal expansion) obtained in air 
and in argon, using the electrical method, with a 
standard deviation of 0.30 percent: 


A\=0.716— 0.247 X 10-°t + 0.182 X 10- St? 


— ().783 x 10 


1073 


(47) 


where ¢ is temperature in °C and A is in W/cm deg. 
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TABLE 4. Experimental values for the thermal conductivity of 


platinum as measured using the electrical method 


The values given are corrected for thermal expansion. 
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Argon 

Argon 
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TABLE 5. 


Typical set of data from measurements of the thermal 


conductivity of platinum by the electrical method 


The data correspond to Test No. 6 in table 4 
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800 


Percentage departures from eq (47) of the thermal con- 


ductivity data points obtained by the electrical method. 
‘ 


This equation is valid from 100 to 900 °C. Deviations 
of the data points from eq (47) are plotted in figure 8. 
The scatter is random and there are no significant 
differences between values obtained on heating and 
cooling or between the first run and the second run. 
The value obtained in helium at 200 °C agrees within 
0.05 percent of the value given by (47) but the measure- 
ment in helium at 400 °C shows a deviation of + 0.9 
percent. The measurements in helium were made 
during the first run and there is a possibility that the 
powder insulation did not uniformly and completely 
fill the necked-down region. If the powder did not 
completely fill the space, the boundary conditions 
assumed in deriving the heat loss corrections would 
not have been met. Any errors in the heat loss correc- 
tion would be amplified in helium, particularly at the 
higher temperature, due to the high thermal conduc- 
tivity of this gas. In the second run the powder was 
carefully packed around the neck so that the necked- 
down region was completely and uniformly filled with 
insulation. The two data points obtained in helium were 
not used in deriving (47). 

The thermal conductivity values corresponding to 
(47) are believed to be uncertain by not more than 2 
percent over the temperature range 100 to 900 °C. 
The estimated uncertainties arising from the various 
known sources of error are summarized in table 6: 
these uncertainties were discussed in section 5.2. 


TABLE 6. Summary of individual uncertainties contributing to 
the overall uncertainty in the thermal conductivity results obtained 
using the electrical method 


Uncertainty, ‘ 
of uncertainty 


900 * 


Geometrical factor (Ro/p. 
Temperature 
Slope of }* 
Heat loss correction (tests in a 
Seebeck effect 


coefhcient of resistance 
versus R* curve 

re 
Temperature coefficient of Fermi energy 


Combined 


Each 


squares of the individual uncertainties 


combined uncertainty was obtained by taking the square root of the sum of the 


6. Comparison of Results With Other 
Investigations 


The smoothed experimental results given in sections 
1.3 and 5.3 are tabulated in table 7. The electrical 
resistivity values correspond to equation (45) from 
0 to 700 °C and to (46) at 800 and 900 °C. The thermal 
conductivity values for the longitudina: heat flow 
method and the electrical method correspond to eqs 
(20) and (47), respectively. The values for the Lorenz 
function, Ap/T, were computed from the smoothed 
values for the electrical resistivity, p. the two sets of 
values for thermal conductivity, A, and the absolute 
temperature, 7. 

In this section the results of other investigations on 
platinum are compared with the results of the present 
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investigation. Since many investigators do not correct 
their thermal conductivity and electrical resistivity 
results for thermal expansion, it was necessary to 
convert all data to a common basis in order to make 
valid comparisons. 


TABLE 7. Smoothed experimental values for thermal conductivity, 
electrical resistivity, and Lorenz function of platinum 

All values are corrected for thermal expansion 

¥ 
Lorenz function 
V*/deg? 


Thermal conductivity 
Electrical W/cm deg 
resistivity a = ae — 
BQ cm Longitu- | Electrical} Longitudinal 
dinal method method 
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Electrical 
method 
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Electrical resistivity of platinum: comparison of pre- 
vious data with the results of the present investigation. 
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FIGURE 9. 


@ Powell and Tye {1} 


ratio of the ice point to helium point ‘resistivity of their No 


Reported purity was 99.999 percent for both samples tested. The 
1 specimen, which was 
purer than their No. 2 specimen, was previously reported by Powell. Tye, and Woodman 
[11] to be about 760 The uncertainties shown represent possible round-off error due to 
results only being reported to nearest 0.1 w#Q cm. 

Laubitz and van der Meer {30|. The reported purity was 99.999 percent. The ratio of 
ice point to helium point resistivity was about 1890. 

Roeser [31]. His sample was probably thermometric grade platinum having a purity 
of about 99.999 percent 

Martin, Sidles, and Danielson [37]. The values shown correspond to their sample which 
had a reported purity of 99.999 percent and a ratio of ice point to helium point resistiv 
ity of 5000 


6.1. Electrical Resistivity 


The electrical resistivity values reported for plat- 
inum by Powell and Tye [1], Laubitz and van der Meer 
[30], Roeser [31], and Martin, Sidles, and Danielson 
[37] are compared with the results of the present in- 
vestigation in figure 9. The reported purities and 
resistivity ratios for the different samples are given 
in the figure caption. 


6.2. Thermal Conductivity 


The thermal conductivity values reported for plat- 
inum by Powell and Tye [1], Laubitz and van der 
Meer [30], and Martin, Sidles, and Danielson [37,.60| 
are compared with the results of the present investi- 
gation in figure 10. The base line in this figure is a 
weighted average of the two sets of data obtained 
in the present investigation, the data from the lon- 
gitudinal heat flow method being given twice the 
weight that was given to the data from the electrical 
method, which had a larger uncertainty. The derived 
curve of Slack [2] is also shown in figure 10; this 
represents essentially all of the thermal conduc- 
tivity values reported for platinum prior to 1962 with 
the exception of the data of Krishnan and Jain [10]. 

The previously unreported data of Watson and 
Flynn [32] shown in figure 10 were’ made at NBS on 
the bar from which the specimen for the present in- 
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FIGURE 10. Thermal conductivity of platinum: comparison of pre- 
vious data with the results of the present investigation. 

The curve in the upper drawing and the baseline in the lower drawing correspond to a 
weighted average of the two sets of data from the present investigation with the data from 
the longitudinal heat flow method being given twice the weight of the data from the elec 
trical method 

Present investigation — longitudinal heat flow method. 
Present investigation — electrical method. 
¥ Watson and Flynn [32|—same sample as that used in the present investigation. 
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vestigation was later machined. Watson and Flynn 
made their measurements in air in a longitudinal heat 
flow apparatus which has been described by Watson 
and Robinson [34] and by Ginnings [35]. 

Several recent sets of thermal conductivity values 
for platinum at lower temperatures are displayed in 
figure 11. The data of Powell and Tye [1]. curve D, 
were also shown in figure 10. The data of Bode [38]. 
curve E, were obtained using a longitudinal heat flow 
method: Bode stated that he considered his values 


¥ the longitudi 
f the sar 


f the sample 


of the s 


the samules-of f we Ive 


to be accurate to within 0.5 percent. Since the present 
investigation was begun, Powell. Tye. and Woodman 
[39], curve F, have used a longitudinal heat flow method 
to make low temperature thermal conductivity meas- 
urements on one of the samples of Powell and Tye. 
Also during this time, Moore and McElroy [40]. curve 
G, have made thermal conductivity measurements 
on the same sample which Powell, Tye. and Wood- 
man [39| used for their investigation. In the tempera- 
ture range shown in figure 11, the thermal conduc- 
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tivity of platinum could be affected by small differences 
in purity and this may account for some of the dif- 
ferences shown. However, the measurements of 
Powell, Tye, and Woodman [39] and those of Moore 
and McElroy [40] were carried out on the same sam- 
ple; hence the differences between the results of 


these investigators must be presumed to be due to, 


experimental errors. 

Additional references to previous measurements 
of the thermal conductivity of platinum are given by 
Powell,, Ho, and Liley [8]. Ciszek [61] has recently 
reported some high temperature thermal diffusivity 
measurements for platinum. 

Values obtained for the Lorenz function of plat- 
inum are compared with those of other investigators 
in figure 12. The present measurements indicate a 
Lorenz function increasing above the theoretical value 
as the temperature increases. Laubitz and van der 
Meer [30] point out that such behavior of the Lorenz 
function can be understood qualitatively if one as- 
sumes a low Fermi energy for platinum. 


7. Conclusions 


The good agreement among the thermal conduc- 
tivity values obtained for platinum by two different 
methods in the present investigation and by an in- 
dependent measurement by Watson and Flynn [32| 
lend considerable weight to these results. The values 
indicate that the thermal conductivity of platinum 
increases with temperature as found by Laubitz and 
van der Meer [30| rather than being essentially in- 
dependent of temperature as found by Powell and 
Tye [1]. However, there remains the possibility that 
there is a real difference between the samples of Powell 
and Tye [1] and those of other investigators. In order 
to explore this possibility, measurements should be 
made at high temperatures in another laboratory on 
one of the actual samples of Powell and Tye. Steps 
are being taken to see if this can be done. 

The results of the present investigation indicate 
that electrical methods of measuring thermal con- 
ductivity can yield equivalent results to those obtained 
by the more conventional nonelectrical methods, at 
least for the conditions of this investigation. This can 
be interpreted as indicating that the thermal conduc- 
tivity of platinum does not depend significantly on 
electric current densities in the range than 
10?A/em?. 

Before platinum can be established as a thermal 
conductivity reference standard, additional measure- 
ments by the same method on samples of differing 
purity are indicated. It is intended to repeat the 
measurements described in this paper on a sample 
of platinum of higher purity than that used in the 
present investigation. 
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8. Appendix A. Evaluation of D(z,, 22) 


In this appendix, the factor D(z,, z2) is evaluated. 
This factor was needed to correct for heat exchange 
between the specimen and the surrounding insulation 
in the longitudinal heat flow method. 

In order to evaluate D(z,;. z2) one must know the 
temperature dependence of the thermal conductivity 
of the powder insulation as well as the radial tempera- 
ture gradient at the surface of the specimen. The ther- 
mal conductivity of the insulation was assumed to be 
a linear function of temperature; 

K= Ky(1 + 609), (A-1) 
where Ko and 69 are evaluated at 6=0, 6 being the 
temperature coefficient of the thermal conductivity 
of the insulation. Then, from (14) 


é=60+ 


bo g2 (A-2) 


The radial temperature gradient at the surface of 
the specimen was determined by analysis of the 
hollow cylinder of powder insulation between the 
specimen and the guard with known temperature 
distributions at its boundaries. At the surface of the 
specimen, r=a, the potential distribution was rep- 
resented by 


NWTZ 


— go) + A, sin’ 
@ a @ (A 3) 


oO, + Oo 
BO Sw 


and at the inner surface of the guard, r=, the poten- 
tial distribution was represented by 


NTWTZ 


(=), b Ny (hy a hy) + >. By sin — 
oI o  (A-4) 


where €(r, z) has the values (a. 0) = go. E(a, @)= 8. 
E(b, 0) =ho. E(b, w) =h,, at the ends of the region, 
z=0 and z=c, which correspond to the outermost 
thermocouples on the guard. At these ends the radial 
temperature d_ -tribution in the powder insulation was 
assumed to be !ogarithmic: 


In r/a 
(ho— go) 7 7 
In b/a 


In r/a 


Bat (he (A-6) 


= Bee) 


In b/a 


Fourier coefficients 4, and B, in (A-3) and 
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(A-4) were determined from the measured tempera- 
ture distributions along the specimen assembly and 
the guard cylinder. In the gradient region of the 
specimen, and in the upper and lower specimen ex- 
tensions, the observed temperatures and thermo- 
couple locations (corrected for thermal expansion) 
were used to derive quadratic equations giving 6 
versus z, where @ was measured relative to the tem- 
perature at the center of the gradient region of the 
specimen. In the neck region, which was nearly iso- 
thermal at all times, a linear equation was used. The 
three quadratics in 6 transformed to quartics in € and 
the linear equation in 6 transformed to a quadratic in 
é with application of (A—2). In the intervening heater 
regions, smoothing cubics [27] were used which pro- 
vided continuity of temperature and longitudinal tem- 
perature gradients. Thus the &-distribution along the 





hes (hyo — go) 


In r/a ae )4 | ta (ha — Sw) 


In b/a \ @ 


In r/a 


In b/a 


{,F o(nar/w: 


specimen was described by seven smoothly joined 
polynomials: quartic (lower extension), smoothing 
cubic (lower extension heater), quartic (gradient region 
of specimen), smoothing cubic (specimen heater), quad- 
ratic (neck region of specimen), smoothing cubic (up- 
per extension heater), quadratic (upper extension). 
Equating this set of polynomials to the right-hand side 
of (A-3) and making use of orthogonality enabled calcu- 
lation of the 4,’s. A similar set of polynomials was used 
to represent the &-distribution along the guard, from 
which the B,’s were calculated. It would serve no 
useful purpose to explicitly display the expressions 
for A, or B, in this paper; they are quite lengthy and 
rather complex 

With the boundary conditions (A-3, (A-4), (A-5), and 
(A-6) the potential, €, at any point in 


the powder insu- 
lation is given by 


n7rb/w) — B,F o(nar/@: nza/w) 





t > 
ha 
n=1 

where 


F'o(x, .) ) = To(x )Ko(s ) — Ko( x Mos he 


— sn — 
Fy(nzra/@: n7b/w) 


I, = modified Bessel function of first kind and order m, 


K,, = modified Bessel function of 


Substitution of (A-7) into (16) yields D(z,. z»): 


D(z. z2) =27a (ho — g0) (22 —21) + (ho —ho— ga + 2 


{,F ,(n7a/w: nab/w) — (w/n7za)B,, 


second kind and order m. 


2m ja ln bla 


NW7Z2 





a 
z, 
n=1 


where 


F(x. vy) =11(x%)Ko(y) + Ki(x)Io(y). 


Fo(n7ra/w: n7zb/w) 
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9. Appendix B. Evaluation of Heat Loss 
from Heater Jumpers 


As discussed in section 4.2b, the “jumpers” shown 
in figure 3 were heated by conduction from the helical 
heater elements to a temperature somewhat above 
that of the adjacent portion of the specimen. Since the 
smoothing functions [27] used to represent the tem- 
perature distribution along the specimen did not allow 
for this local hot region, an uncorrected-for heat loss 
from the jumpers into the powder insulation must have 
occurred. The purpose of this appendix is to estimate 
the possible magnitude of this heat loss. 

We will assume that the five hot jumpers can be 
represented by equivalent hot regions at the convex 


surface, r=a, of the specimen. By integrating over the 
angular temperature distribution near the longitudinal 
position of the heater, we can obtain an effective aver- 
age longitudinal temperature distribution which can 
be used to estimate the heat loss from the jumpers. 
We assume that, for a gradient test, the effective longi- 
tudinal temperature distribution was as shown in the 
curve labeled g'(z), in figure 13. In our analysis, it 
was assumed that the temperature distribution was 
as shown in the curve labeled g(z) in figure 13. The 
difference between the effective longitudinal tempera- 
ture distribution, g’(z), and the assumed temperature 
distribution, g(z), is shown in the curve labeled g”(z) 
in figure 13. 

If we go through an analysis similar to that in ap- 
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yf 


Longitudinal temperature distribution in the vicinity 
of the specimen heater. 


FIGURE 13. 


The three different curves are identified in the text 


pendix A, but using g’(z) rather than g(z) to represent 
the longitudinal temperature distribution along the 
specimen, we obtain 


D’ (2, 22) = D(z, z2) + D"(z, 22), 


where D'(z;. 22) 


corresponds to eg (z), D (21, 22) is 
given by (A-8), 


and 


Fy (naalw: nizb/w) 


D" (2, Zz») 
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NT7zZ\ 
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where the Fourier coefficients A)’ are defined by 


: _ nwz 
e"(z) sin dz. (RB 
W@W * 


Let us assume that ¢”(z) can be represented by 


o"(z) GC. 


e'(z) Q. elsewhere 
where ¢ is the location of the center plane of the heater 
which has a thickness of 2e€. Substitution of (B—4) 
into (B—3) yields, after integration, 


1G . nit NTE 
ae” lh 


niT @ @W 


he a 
‘tn 


(B-5) 


Using (B—5) for A”, evaluation of (B—2) for the geom- 
etry of our apparatus yields D"(z,, z2) = 1.2G (deg cm). 
An upper limit for the temperature of the jumpers 
was provided by the average temperature of the heater 
winding. which was determined using the heater as 


a resistance thermometer. The jumpers shown in 
figure 3 occupied less than 25 percent of the perimeter 
of the specimen at the midplane of the specimen 
heater. Using this information an upper limit for G 
was determined for each of the tests in table 1. Using 
this in conjunction with the measured values of « and 
the value for D"’(z,. z2) obtained above, the maximum 
error due to the uncorrected-for heat loss from the 
jumpers was computed to increase from 0.2 percent 
at 100 °C to 0.5 percent at 900 °C. 


10. Appendix C. Theory of the Electrical 
Method 


Although the theory of electrical methods of measur- 
ing thermal conductivity had been developed in the 
literature, there was no single reference that ade- 
quately covered the subject. In the present investiga- 
tion, it was found necessary to digest a rather large 
number of papers in order to arrive at an understand- 
ing of the validity of the expressions used for comput- 
ing the thermal conductivity values. It is hoped that 
the fairly unified treatment given in this appendix will 
be useful to future workers in this area. 

In the electrical method the specimen supports 
an electrical potential gradient, as well as a tempera- 
ture gradient, and the corresponding thermoelectric 
effects must be analyzed. This is done most conven- 
iently by the methods of irreversible thermodynamics, 
applying the Onsager Reciprocal Relations [44—54}. 
The application of the relations to steady-state proc- 
esses involves an approximation; but as Callen [44] 
points out, it is an excellent approximation and is 
completely justifiable in this case. 

We define a set of current densities J;: 


Ji=¥ LX. 


} 
such that 


S Ji: Xi. (C 


i 


where X; are the “conjugate forces,” 
R(S) 
where R(S) is the rate of production of entropy in the 


system. Then the Reciprocal Theorem 
53, 54] states that 


Onsager 


Lij = L ji (C3) 
in the absence of a magnetic field. The rate of entropy 
production R(S) is uniquely defined by the system 
under consideration, but since R(S) can be = split 
into a sum of products in many ways, one is left with 
a choice of current densities and conjugate forces [45]. 

We define an electric current density J, an energy 
current density W. and an entropy current density S, 
so that the divergence of each of these current den- 
sities is the rate of change per unit volume of the 
corresponding thermodynamic variable. With these 
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definitions for the 


[44, eq 12]: 


current densities we can write 


Ts=w-yj. 


e 


where T is the absolute temperature. u is the electro- 
chemical potential, e is the electronic charge and 
J. the electric current density. is equal to the electron 
current density times the electronic charge. We are 
considering the electron current to be the only mass 
current. Mass transport and the associated entropy 
Hows due to electromigration (preferential migration 
of ions in a solid when a direct current is passed 
through the solid) and due to the Soret effect (mass 
transport in a solid due to a temperature gradient) 
are not considered. The reader is referred to | 11-43 ] 
for discussions of these two effects. 

We are concerned with steady-state conditions for 
whic h: 


Vw 
vs 


From (C—4) and (( 


Using the identity 


Vv (+) —VT, 
1 ‘Ke 


we can write the 


forces” 


current sities In 


“conjugate 


Lu - 
‘lev! 


where @ is the electrostatic potential, and @ is the 
chemical potential, which in the case of a metal is 
simply the Fermi energy [55]. Since the Fermi energy 
is a function of temperature only, we can write 


S WT + Vo. (C-13) 
e aT 
We now proceed to evaluate the coefficients re Lys: 
and Lo». 

The electrical conductivity is defined under iso- 
thermal conditions (i.e.. V7 =0) as 


Co 


fe. 
lve 
Q. Vip e) 


since 


for 
lows that 


V/ Vo. From (C 


10) it fol- 


(C—]5) 


The absolute thermoelectri: power is defined (see [56]) 


as 


Vile) ] 
Vl 


and from (C-10) it follows that 


he thermal conductivity is d 


W] 


| 


LVF 4s 


current 


[—Q | 
LWT is 


Q is the heat density. when 


since 


1. Lye. and Ls in (4 
1 W in 


terms of the 


10) and (C-11). 
defined thermo- 


w=|4+75 


AVT. 


21) 


Taking the divergence of W 
(C—6). 


and 


using (C—5) and 


V-W=J -V(ulet+ TS) — V-AVT=0. 


Therefore 


J+ [V(ule) + TVS + SVT] —V-AVT=0. 


(C—23) 
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In a homogeneous and isotropic medium, the param- 

eter S*"s is a function of temperature only, and we can 

write 

J: (Viule) + (TaS*s/dT) YT + SYT)|—-¥ - AVT=0. 
(C-24) 

The coefficient TaS*"’/dT is defined [44] as the Thom- 

son coefficient 7. From (C-20) and (C-24) 


J: (-J/o+7V7T\-—V -AVT=0. 


At this point, it is convenient to introduce a pseudo- 
potential & such that 


oV ws. (C—26) 


and 


Vb =Viule) + SVT. 


Writing (C—25) in terms of Y we get 
VeAVr +o V+ or¥u -Wr=0. 
(C-28) 


Equation (C-28) can be considered as the generalized 
equation relating the temperature distribution and 
the pseudo-potential distribution in an isotropic ho- 
mogeneous medium carrying an electric current. We 
now proceed to solve that equation. 

The second term in 


V-(WoVw) since 
Vv ‘a Uo Vu \ wV : (cVv) 4 ov ” Vu 
=o VU - Vv. since V-cVY =-V -J=0 (C—29) 


The third term in (C-28) can also be written in the 
form of a divergence, since 


(C-28) can be written as 


ie 
oVw V| 


J, td 


v-[o(| 


7 


dt) 98) 
[’ TdT ) Vv: (cV ws) 


ih | 
oVh VY | rdT=o79h VT. 


J! 


where 7™ is any arbitrary fixed temperature. 
Using these identities (C-—28) can be written as 


Vv: Awe oby+o (| 


1 


7 
rdT | ve | =0. (C 


Integrating, 


1 
VT +a Vi+o | | rd) Vv=Verte. 


a1 


where € is any potential satisfying the Laplacian 
V?é=0 and ¢ is a vector constant. Since Y -(oVus) =0 
we can write 


VéE=WoVy, 


where Wo is a constant which must satisfy the boundary 
conditions. Substituting for Yé in (C-32), we get 
C7 
AV T+ o0( b— Wo) Vb + (co Vu | TdT =e. (C—34) 
1 

If c=0, it can be seen from (C—34) that VT and Vv 
must be parallel at all points. This implies that if there 
is no electric current flow across a boundary, there 
can be no heat flow across that same boundary. More- 
over, equipotential surfaces are isothermal. 
More explicitly, the condition that e=0 requires that 
all the heat generated in an electrically insulated con- 
ductor must flow out at the ends of the conductor, 
there being no heat losses from the sides. This con- 
dition can never actually be met in practice since there 
is no perfect or even near-perfect thermal insulator. 
However, by appropriate choice of geometry, as 
employed in the present experiments, heat losses 
can be made very small so that we can take c= 0 as 
a valid statement of the boundary condition. Equation 
(C-34) then reduces to 


also 


1 
\pVT + (b—Wy) V+ (Vw | TdT =0, 
(C-35) 
where p is the electrical resistivity, p= 1/o. 

If there is a point in the medium where the tempera- 
ture has a’maximum value, 7). the gradient WT is zero 
at that point. Further, if we let 

T*=Tm» (C-36) 
the first and third terms in (C-35) are zero at T=T,,. 
Since there is a current flowing, cY Ww #4 0 and (C35) is 
satisfied only if &/— W)=0 at T=T,,. In other words 
w= Wy at T=T),. -37) 


so that ao is the value of Ww at the point of maximum 
temperature. 

Integrating (C-35) along a line from one isothermal 
surface S$,(7;. w,) to another isothermal surface 
So(T2. We.) and substituting 7,, for T 


S 


; ApVT- dr T | ; ( Ww — Wo) Vw f dr 


JS 


So rl 
+ | Vu “de | tdT=0, (C-38) 
+4 ' Tm» 


dT and 


where r is a position vector. Since Y7'- dr 


Vw -dr=dwu, we have 


rT, 2 CW. 7 
| ApdT + | (Ww— bo)db +| dus | rdT=0. 
1 hy hh Tm 


(C—39) 





FIGURE 14. A thermally and electrically 


arbitrary geometry. 


insulated conductor of 


An arbitrary isotropic and homogeneous conductor 
is represented in figure 14. 

There is no electric current flow across the bound- 
ary € and the geometry is selected to make heat 
losses from C negligible so that the condition for e=0 
discussed above is satisfied. The isothermal surfaces 
S, and Ss are taken to be at the same temperature, 
7». and the surface of maximum temperature is repre- 
sented by S,,. In the absence of the Thomson effect. 
the surface of maximum temperature would bisect the 
reduced resistance between S,; and S,. the reduced 
resistance—a geometrical quantity—being the re- 
sistance of the medium for electrical conductivity of 
unity. However, the Thomson effect moves the surface 
of maximum temperature away from the symmetrical 
location, this effect being referred to as the “Thomson 
shift” [57]. It is analyzed in detail by Davidson [57| 
and discussed by Holm [58] and Llewellyn Jones [59]. 
It is convenient to define W=0 and consider the inte- 
gral (C—39) from S,, where T=T,, and /=0 to any 
other surface S where J=wW. T=T. Rearranging we 


have 
Tm —_ 
TdT. 
7 


(C—40) 


el - 


wu ly» i 
| wd = | ApdT + | dv 
0 J7 J0 


It is convenient to write (C—40) in the following form: 


W —2ew — F?=0 (C-41) 


where F(T, Tm) =+ 


2 oe hod | 


and 


] y rT» 
e(T, Tn) =3 | db | rd. 
Ww 0 1 


Solving the quadratic equation (C—41), we g 
Ww=et (E+ F*)'?=e+F (1+ (€/F)’ 
(C—44) 


The Thomson effect is relatively small so that the 


second-order term (€/F)? in (C—44) can be neglected. 
Then 


w=exrF. (C—45) 
the plus and minus signs referring to the high and low 
potential sides of S,, respectively. 
Referring to (C—43) we see that although T is a fune- 
T 


. . ° . n et e . 
tion of W. the sign of TdT is independent of the 
Ji 

sign of & since 7), is always greater than 7. Conse- 
quently € is. to the first order, an even function of w. 
having the same value at corresponding points on the 
high and low potential sides of S,,. At S; we have 


Us E | To. ci ) F( To. fen ) (¢ : 16) 


and at S» 


Therefore. 


Us - Us, 


2F (To. Tm). 


From (C 


27) and (C-13) we have 


Vu=yo+! “WT 4 SY. 
ead 


Integrating from S,; to So. 


db» — db; 4 (Z2- 
e 


We — Wy 


—50) 


In an isotropic and homogeneous medium ¢ and S#"S 
are functions of temperature only. When the terminal 
temperatures are the same, as in the case we are con- 
sidering; i.e.. when 7)=72=T». 
(C—48) and (C—50) 


we can write from 


We — Ui = ho — od; 


2F(To.Tm), (C-51) 
where V is the voltage drop between the surfaces 
S, and S:. If the potential probes are not at the same 
temperature, a correction must be made as indicated 
in (C—50). Writing F explicitly we finally get from 
(C-51) and (C-42): 


V2 [Tm Am 
B | ApdT | Apdé, (C=32) 


J4o 


where 06=T—T , 0m =Tm—To. 

We now proceed to relate the maximum temperature 
rise in a conductor to the measured resistance of the 
conductor. This problem has been discussed by Meiss- 
ner [20] and by Holm [58]. 

Consider two geometrically identical conductors 
differing only in that one (real) conductor has a finite 
temperature-dependent thermal conductivity, \=A(@), 
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while the other (hypothetical) conductor has an in- 
finite thermal conductivity. Both conductors are 
assumed to have the same temperature-dependent 
electrical resistivity, p=p(@). We assume the lateral 
surface of each conductor to be insulated perfectly 
against the flow of both heat and electricity. Let each 
end of both conductors be in thermal equilibrium 
with a heat sink at a temperature 6=0. If an elec- 
trical current is passed through the real conductor, 
there will be an accompanying temperature rise due 
to the heat generated by the passage of the current. 
This in turn will cause a change in the electrical 
resistance of the conductor. Since the hypothetical 
conductor is assumed to have an infinite thermal 
conductivity, the temperature throughout this con- 
ductor will be @6=0, even in the presence of heat 
generation, and the electrical resistivity will have the 
constant value p= p(0)= pp at all points in the hypothe- 
tical conductor. We wish to consider the difference be- 
tween the electrical resistances of these two con- 
ductors due to the fact that one of them (real) is heated 
to a higher temperature due to having only a finite 
thermal conductivity. 

Let & and @ be the potential and the temperature 
respectively in the real conductor and let Wo and 4 
refer to the corresponding variables in the hypotheti- 
cal conductor. If the same electrical current is passed 
through these two conductors, the temperature and 
potential distributions in the two conductors will dif- 
fer due to the temperature-dependence of the proper- 
ties in the real conductor and also due to the thermal 
expansion of this conductor. If dws and dio are the 
potential differences between corresponding equipo- 
tential surfaces in the two conductors, we can, from 
(C—26), derive the relation 


Life 0 
<= 8 ( 


+ oO), 
dus p Yo 


(C—53) 


where yo is the coefficient of linear thermal expansion 
at 0=0. Neglecting terms arising from the Thomson 
effect, we have from (C—41), 


(C-54) 


) 


and from (C35), 


Combining the last three equations, 





dj =+ Apo (I = yo) dé 


Integrating, 





roi Apo(1 + yo) dé oe 
Wo = F . (C—-57) 


Since both conductors are assumed to carry the same 


current, 
Ro _ wo 


Rw’ (C—58) 


and we have 





Ro ] 9m = Apo(l ie 
=; | poll + yob) dd, (C-59) 
0 P 


R / 


where R is the electrical resistance between two sur- 
faces at temperature 6=0 in the real conductor and 
Ry is that between these surfaces in the hypothetical 
conductor having an infinite thermal conductivity. 
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1. Introduction 


Commercially pure iron, usually produced by the 
American Rolling Mill Company and known as Armco 
iron, has been used for many years as a thermal 
conductivity reference material, either as a “heat 
flow meter” in the case of comparative thermal con- 
ductivity measurements or as’a material to be used 
in checking apparatus in the case of absolute thermal 
conductivity measurements. Powell [1]! has reviewed 
the numerous published data available on the thermal 
conductivity of Armco iron and derived “most probable 
values,” which he estimated to be good within +2 
percent from 0 to 600 °C, and within about + 5 percent 
at 1000 °C. The range of experimental values in the 
literature, however, is about 7 percent between 0 
and 600 °C, and increases to almost 30 percent at 
1000 °C. In an attempt to resolve this discordance of 
experimental results, C. F. Lucks of the Battelle 
Memorial Institute, Columbus Laboratories, proposed 
a round-robin of thermal conductivity measurements 
on Armco iron. The Battelle Memorial Institute (BMI) 
subsequently obtained a quantity of Armco iron and 
distributed samples to: National Physical Labora- 
tory (NPL), Teddington, England; National Research 
Cuuncil (NRC), Ottawa, Canada; and National Bureau 
of Standards (NBS), Washington, D.C., and later to 
a number of other laboratories. Powell et al. |2], have 
reported the NPL measurements over the temperature 
range from — 200 to + 1000 °C; Laubitz has reported 
[3, 4] the NRC measurements from 30 to 1000 °C. 


The present paper reports thermal conductivity 
and electrical resistivity results obtained at the Na- 
tional Bureau of Standards on the sample of Armco 
iron supplied to NBS by BMI. Thermal conductivity 
measurements were made at NBS over the tempera- 
ture range —160 to +640 °C. Electrical resistivity 


' Figures in brackets indicate the literature references at the end of this paper. 


* See appendix for discussion of nonmetallic impurities 


275-286 O-67—3 


thermal conductivity. 


measurements were made over the temperature range 


— 195 to + 1380 °C. 
2. Sample Description 


The sample of Armco iron supplied to NBS by BMI 
was in the form of a round bar, nominally 1 in in diam- 
eter and 39 in in length. The ends of this sample were 
stamped No. 3 and No. 4, the No. 3 end being immedi- 
ately adjacent to the sample sent to NPL, and the 
No. 4 end adjacent to the sample sent to NRC. This 
Armco iron was stated to have been obtained from the 
Steel Sales Company, Chicago, Ill. The ladle analysis 
in weight percent, as given by the supplier, is as 
follows: C 0.02, Mn 0.030, P 0.006, S 0.023, Si 0.004, 
Cu 0.083, Fe (by difference) 99.834. Powell [2] reported 
that an NPL analysis showed no significant variation 
from the analysis given above, other than the presence 
of 0.083 percent Ni. A spectrochemical analysis by 
the NBS Spectrochemistry Section yielded the fol- 
lowing impurity content ? (in weight percent): Mn < 0.1, 
Si < 0.015, Cu 0.06, Ni 0.06, Cr 0.01, V < 0.01, Mo 
0.01, W < 0.02, Co < 0.01, Ti 0.006, Sr 0.02, Nb < 0.01, 
Zr < 0.003. 

Battelle Memorial Institute reported that they had 
annealed the sample at 870 °C for 1/2 hr. The Rockwell 
hardness of the sample, as received, was found to 
vary between B30 and B50. A photomicrograph of 
this specimen, as viewed at 100x, is shown in figure 1. 
The specimen material exhibits large grains with no 
particular orientation. 


3. Thermal Conductivity Measurements 


The thermal conductivity measurements on the 
BMI Armco iron sample were made in the NBS metals 
apparatus which has been described by Ginnings [5] 
and by Watson and Robinson [6]. In brief, the measure- 
ment was made by determining the electrical power 
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FIGURE 1. Microstructure of the BMI Armco iron sample (X 100). 


input to a heater in one end of a specimen 37.0 cm 
long and 2.386 cm in diameter, which was cooled at 
the other end by circulating water or liquid nitrogen. 
Temperatures were measured by means of thermo- 
couples (fabricated from wire which had been cali- 
brated by the NBS Temperature Physics Section) 
3.51 cm apart, along the central portion of the bar, 
thus permitting the calculation of six thermal conduc- 
tivity values, each at a different mean temperature, 
for each thermal equilibrium. A guard cylinder, con- 
centric with the specimen, was used to minimize heat 
exchanges between the specimen and the surrounding 
insulation, and corrections were made for such heat 
exchanges. 


TABLE 1. Test conditions for the thermal conductivity measurements 


Appa Heater 


ratus 


Atmos 
phere 


Number 
of data 
points 


winding Thermocouples Insulation 


Nichrome Chrome! P 


alumel 


Diatomaceous 12 
earth 

Nichrome Platinum — 10 Powdered 

rhodium 

platinum 


alumina 


Platinum — 10% Powdered Vacuum 18 


rhodium alumina 


platinum 


Although diatomaceous earth is a better thermal insulation, it was necessary to use 
powdered alumina in series 2 and 3 in order to avoid contamination of the noble metal 


thermocouples 
Note that none of the values of thermal conductivity and electrical resistivity given in 
this paper are corrected for the effect of thermal expansion of the specimen 
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FIGURE 2. The lower figure shows the data points obtained at NBS 
/ 


for the thermal conductivity of BMI Armco iron. 
The three different sy nditions presented in table 
1. The upper figure shows the percentage departures of the smoothed values of Powell 


et al. [2 


mbols correspond to the different test co 


and of Laubitz [4] from the smooth curve in the lower figure 


Data were taken on the same sample in two different, 
but similar, apparatus with varying experimental con- 
ditions as indicated in table 1. The smoothed NBS 
results.’ which have not been corrected for thermal 
expansion, at 50 °C intervals are presented in table 2 
(along with values for the electrical resistivity and the 
Lorenz function which will be discussed below). The 
curve in the lower drawing of figure 2 represents the 
thermal conductivity values given in table 2; the three 
different symbols represent data points correspond- 
ing to the three series of tests described in table 1. 
For the material which we tested, in the state in which 
we tested it, the thermal conductivity values in table 2 
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TABLE 2. Smoothed values for the thermal conductivity. electrical 
resistivity, and Lorenz function of BMI Armco iron 


Lorenz 
function 


Electrical 
resistivity 


Thermal 
conductivity 


Tempera- 
ture 


Cc it 
200 
160 : 2.é 
150 
100 
50 
0 
50 
100 
150 
200 
250 
300 
350 
L00 
150 
500 
550 
600 
640 
650 
700 
720 
740 
750 
760 
770 
780 
800 
850 


reye) 
6oU 


em deg 


BOD cm 


0 Ul of 








900 
910 
920 
950 

1000 

1050 

1100 

1150 

1200 

1250 

1300 

1350 

1380 


Pec CoOUNONOAUN 











*Extrapolated value. 


are estimated to be in error by not more than two per- 
cent. 

In the upper drawing of figure 2, Powell’s values 
[2| and Laubitz’ adjusted values [4| for the thermal con- 
ductivity of BMI Armco iron are shown as percentage 
departures from the thermal conductivity values given 
in table 2. 


4. Electrical Resistivity Measurements 


The electrical resistivity specimens were placed 
in series with a calibrated standard resistor and a 
regulated d-c power supply. The resistance of each 
specimen was determined by comparing the voltage 
drop across potential taps in the specimen with that 
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FIGURE 3. The upper sketch (a) shows the “‘squirrel-cage” specimen 
configuration used for the electrical resistivity measurements on 


BMI Armco iron. 


Phe lower sketch (b) indicates the path of current flow 


the standard resistor. In order to minimize 
thermoelectric effects, voltage drops in the speci- 
mens were measured with the current flowing normally 
and reversed. All voltage measurements were made 
using a precision d-c potentiometer. 

The electrical resistivity at the ice point, 0 °C, was 
determined by measuring the resistance between two 
knife edges 9.997 cm apart spanning the central por- 
tion of a 0.5003 cm diam bar machined from the ther- 
mal conductivity specimen described above. 

The temperature-dependence of electrical resistivity 
was determined by measuring the electrical resistance 
of two different “‘squirrel-cage”’ specimens (see fig. 3) 
as functions of temperature. Three thermocouples 
peened into these specimens served to measure the 
temperature of the specimen during testing and also 
served as potential taps to measure the voltage drops 
in the specimen. Comparison of the ice-point resistance 
of these two specimens with the ice-point resistivity 
value separately determined on the 0.5003 cm diam 
bar enabled calculation of resistivity values as func- 
tions of temperature. 

The first squirrel-cage specimen of BMI Armco 
iron was fabricated from the end of the NBS thermal 
conductivity specimen which had not been heated 
during testing. Data were taken in air in an isothermal 
cryostat from — 195 to +75 °C with Chromel P versus 


across 
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constantan thermocouples in the specimen and from 
0 to 1380 °C in helium in an isothermal furnace with 
platinum — 10 percent rhodium versus platinum and 
also platinum —30 percent rhodium versus platinum 
—6 percent rhodium thermocouples in the specimen. 
Data were first taken on heating to 850 °C. On cooling 
back down, the electrical resistance of the specimen 
at room temperature was found to be lower by about 
0.1 wQ cm, or 1 percent. On second heating. data 
were taken from room temperature to 1380 °C. Since 
Powell [2| had reported » change in the resistivity at 
910 °C, corresponding to the alpha-gamma transfor- 
mation in iron, data were taken at about 2-deg inter- 
vals from 900 to 920 °C, looking for the effect of the 
alpha-gamma transformation. No effect was found— 
the precision of the data was such that a change of 
0.1 percent could easily have been detected. Study of 
a large scale plot of the data indicated that the elec- 
trical resistivity of the iron had undergone a drop, 
estimated by extrapolation as about one-half percent, 
somewhere between 875 and 900 °C.4 Data taken on 
cooling from 1380 °C indicated a jump in electrical 
resistivity of perhaps two-tenths of 1 percent between 
895 and 890 °C. After cooling from 1380 °C, the elec- 
trical resistance of the specimen at the ice-point was 
about three percent lower than the original resistance 
before the specimen had been heated. 

The second squirrel-cage specimen was fabricated 
from an unused portion of the stock supplied to NBS 
by BMI, annealed in helium at 850 °C for about 1/2 
hr to relieve possible strain introduced during fabri- 
cation, and then cooled to room temperature. Data 
were taken from — 195-to +50 °C in air with copper 
versus constantan thermocouples in the specimen and 
from 0 to 850 °C in helium with platinum — 10 percent 
rhodium versus platinum thermocouples in the speci- 
men. These data were taken in the following order: 
on heating at 50 deg intervals to 700 °C, on heating 
at approximately 1 deg intervals from 745 to 770 °C, 
on cooling at approximately 0.5 deg intervals from 763 
to 749 °C, on heating at 800 and 850 °C, on cooling at 
100 deg intervals from 800 to 0 °C. A slight hysteresis 
(< 0.2%) was found between 754 and 761 °C, pre- 
sumably corresponding to the Curie transformation. 
The ice-point resistance of the sample was found to 
be 0.7 percent less than it had been at the beginning 
of testing. 

The electrical resistivity values given in table 2 
were obtained by combining and smoothing the data 
obtained on first heating of the specimens. No correc- 
tions were made for thermal expansion. For the mate- 
rial which we tested, in the state in which we tested 
it, the resistivity values tabulated are estimated to 
be in error by not more than 2 percent, or 0.1 wO cm, 
whichever is greater. Most of this estimated uncer- 
tainty arises from the vagaries of this material, rather 
than from experimental error. All electrical resistance 


*Since our measurements were made, Fulkerson, Moore. 
their electrical resistivity measurements on BMI 
gamma transformation was between 893 and 898 °¢ 


and McElroy [7| reported that 
Armco iron indicated that the alpha 
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FIGURE 4. The lower figure represents the data cbhiained at NBS 
for the electrical resistivity of BMI Armco iron on first heating of 
a specimen which previously had been annealed at 850 °C for 
Yo h 

2 hr. 


The upper figure shows the percentage dep 
et al 


artures of the smoothed valu f Powell 


2|. Laubitz [3]. and Fulkerson et al. [7|. from the curve in the lower figure 


measurements are believed to have been accurate to 
+0.1 percent or better. The mean temperature of the 
sample is believed to have been known within +3 
deg below 1100 °C and within +5 deg at 1380 °C. The 
uncertainty in resistivity directly attributable to the 
uncertainty in mean temperature would be about 0.8 
percent just below the Curie point and would be less 
at all other temperatures. 

The curve in the lower drawing of figure 4 represents 
the electrical resistivity values given in table 2. The 
location of the Curie temperature is indicated by the 
arrow labeled 7,.; the break in the curve corresponding 
to the alpha-gamma transformation is indicated by 
the arrow labeled a—y. In the upper drawing of figure 
4, the values of Powell [2], Laubitz [3], and Fulkerson, 
Moore, and McElroy [7] for the electrical resistivity 
of BMI Armco iron are shown as percentage depar- 
tures from our values as given in table 2. 
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5. Lorenz Function 


The values obtained for the Lorenz function, Ap/T 
(where A is thermal conductivity, p is electrical re- 
sistivity, and T is absolute temperature), are given in 
table 2. The curve in the lower drawing of figure 5 
represents these values. In the upper drawing of figure 
5, Powell’s values [2] and Laubitz” values, which we 
have computed from his electrical resistivity values 
as reported in [3] and his revised thermal conductivity 
values as reported in [4], for the Lorenz function of 
BMI Armco iron are shown as percentage departures 
from the values given in table 2. 


6. Comments 


C. F. Lucks of the Battelle Memorial Institute is 
preparing a paper in which he will discuss in detail 
the results of the round-robin thermal conductivity 
measurements on BMI Armco iron. In addition to the 
results reported in this paper, and those of Powell [2] 
and of Laubitz [3. 4] which have been published pre- 
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FIGURE 5. The lower figure represents the data obtained at NBS 


for the Lorenz function of BMI Armco iron. 
The upper figure shows the percentage departure of the smoothed values of Powell et al. 


[2|. and Laubitz |3, 4| from the curve in the lower figure 


viously, Lucks will discuss other unpublished results 
from other laboratories. In view of this, the authors 
feel there is no need to discuss the above data beyond 
what has been done. 

Powell [1] has summarized thermal conductivity 
measurements on Armco iron made through the year 
1960. Powell. Ho, and Liley [8] presented a compen- 
dium of essentially all of the known thermal conduc- 
tivity data on Armco iron and also on pure iron through 
the year 1965. After the present paper was written, 
Shanks, Klein, and Danielson [11] have published 
values for the thermal diffusivity, specific heat, and 
electrical resistivity of BMI Armco iron. 

In the appendix, some previously unpublished data 
on another sample of Armco iron, not from the BMI lot, 
are presented. 

7. Appendix 

Measurements over the temperature range from 
— 160 to +200 °C were made of the thermal conduc- 
tivity and electrical resistivity of a sample of Armco 
iron submitted by U.S. Army Missile Support Com- 
mand, U.S. Army Missile Command, Redstone 
Arsenal, Alabama. A spectrochemical analysis by 
the NBS Spectrochemistry Section yielded the fol- 
lowing impurity content (in weight percent): Mn < 0.1, 
Si < 0.015, Cu 0.03, Ni 0.04, Cr 0.01, V < 0.01, Mo 
0.01, W < 0.02. Co < 0.01, Ti 0.006, Sr 0.02, Nb < 0.01, 
Zr < 0.003. The Rockwell B hardness of this sample, 
as received, was found to be 72. A photograph of this 
specimen, as viewed at 100X, is shown in figure 6. 
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FIGURE 6. Microstructure of the cold-worked Redstone Arsenal 


sample of Armco iron (X 100). 
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This material exhibits an oriented microstructure 
typical of a cold-worked material, with the long direc- 
tions of the grains parallel to the long axis of the 
specimen. 

The thermal conductivity measurements were made 
using the equipment designated as Apparatus A in 
table 1. 

For this specimen, electrical resistivity measure- 
ments were made in the thermal conductivity appara- 
tus at the temperature conditions existing at the end 
of each pair of runs for determining the thermal 
conductivity, by passing a direct current along the 
bar. Observations were made of the potentials between 
the Chromel P and/or Alumel leads of the span 
thermocouples, with the current direction forward 
and reversed. Due to a slight warming of the bar 
during the period of current flow, the average re- 
sistivity for a span was assigned to correspond to the 
time-average of the span mean temperature over this 
period. A separate measurement on the same specimen 
was carried out in an ice bath using the thermocouple 
leads as potential taps. 


TABLE 3. Thermal conductivity, electrical resistivity, and Lorenz 


function of the Redstone Arsenal sample of Armco iron 


| 
| 
| 
| 





The values obtained for the thermal conductivity, 
electrical resistivity, and Lorenz function of the Red- 
stone Arsenal sample of Armco iron are given in table 3 
for the temperature range — 160 to +200 °C. These 
values are shown in figure 7 as percent departures 
from the values given in table 2 for the BMI iron. 

In general, cold-working raises the electrical resis- 
tivity and lowers the thermal conductivity of a given 
material, the effect on thermal conductivity increasing 
at lower temperatures (see, for example. the work of 
White [9] on gold, silver, and copper). Thus. on the 
basis of the cold-worked state versus the annealed 
state only, assuming identical chemical composition, 
the thermal conductivity of the Redstone Arsenal 
specimen would be expected to be lower than that of 
the BMI specimen rather than higher, as was found. 
The electrical resistivity of the Redstone Arsenal 
specimen was significantly lower than that of the BMI 
specimen (which had an ice-point resistivity of 9.88 
pO cm as compared to 9.36 wO cm for the Redstone 
Arsenal specimen), implying greater purity of the 
Redstone Arsenal specimen, and confirming the find- 
ing of a higher thermal conductivity for it than for the 
BMI specimen. 

Godfrey et al. [10], report (pp. 26-29) that a quanti- 
tative chemical analysis showed the presence of 0.086 
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ing values obtained for the BMI sample 
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conductivity. 


percent O.. 0.023 percent S, and 0.013 percent C (all 
weight percentages) in an Armco iron sample which 
they investigated. They report the presence, in this 
sample, of about 0.9 volume percent of a second phase, 
presumed to consist of oxides. sulfides, and phos- 
phides. On a sample of BMI Armco iron, they found 
1.3 volume percent of a second phase. The amount of 
nonmetallic impurities in a sample of Armco iron and, 
as pointed out by Godfrey et al. (p. 40), also the physical 
state of the impurities, may be of quite significant 
importance as regards the thermal and electrical con- 
ductivities. Impurities in solution would be expected 
to'have a much greater effect on thermal and electrical 
conductivity than would impurities present as a dis- 
persed second phase. 

\ significant difference in purity between the Red- 
stone Arsenal specimen and the BMI specimen could 
only be explicitly determined by a more complete 
chemical analysis of both specimens, including quan- 
titative analysis for nonmetallic impurities. A detailed 
microstructural analysis would also be required to 
determine the physical state of the impurities present. 
It is interesting whether the orientation of the micro- 
structure in the cold-worked specimen would result 
in anisotropy in the electrical and thermal conduc- 
tivities of the metal in its present state. 
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Expressions are developed which permit calculation of the temperature-dependent thermal con- 
ductivity of a cylindrical specimen in which heat is generated internally, e.g., by radioactive decay. 
The information needed consists of the experimentally determined temperature distributions on the 
surfaces of the cylinder, the heat flow through a central circular area at one end of the cylinder, and 
the rate of internal heat generation (which in general may be position-dependent). Numerical coefficients 
are tabulated for the case of uniform internal heat generation. The application of this calculation pro- 
cedure to published methods of thermal conductivity determination is shown and an example is given. 

Key Words: Heat conduction, heat generation, heat transfer. neutron absorption, radioactive 
decay, thermal conductivity. 


Introduction where temperature above an arbitrary datum is de- 


noted by the symbol. v. the temperature-dependent 
thermal conductivity by /(v), and the radial, longi- 
tudinal, and angular coordinates by r. z, and 0, respec- 
tively. Equation (1) can be reduced to a simpler form 
by introduction of a new variable, u. defined by the 
relation 


Subsequent to the publication of a paper [1]! in 
which a mathematical analysis of heat flow in a right 
circular cylinder with arbitrary tempezxature boundary 
conditions was made, we were requested by D. L. 
McElroy of the Oak Ridge National Laboratory to con- 
sider the problem of determining the thermal conduc- 
tivity of a cylinder in which internal heat generation ieee : 
is also present. Since most methods of experimentally * hiv’ )de’, (2) 
measuring thermal conductivity are not readily applica- 
ble to self-heating specimens, and since there is a real 
need for determining the thermal properties of radio- 
active materials or of materials which are being irradi- 
ated, we felt the present analysis to be of sufficient 
interest to warrant publication even in the absence of /7u= > ae = 
any experimental data. ar’ ror az’ 0 k* 


where /:* is the value of 4(v) at v=0. Making the sub- 
stitution (2), (1) reduces to 


ape 


Au lou, au, ldu_ Wr) 


2. Mathematical Development For Dirichlet boundary conditions (i.e., the tempera- 
ture distribution on the surface is specified), the solu- 
Consider a homogeneous, isotropic, opaque, solid tion of (3) can be written as 
right circular cylinder of radius, 6, and thickness, /, 
with position-dependent heat generation, W(r), per u=yrw, 
unit volume per unit time, where r is the position vec- 
tor. In general, the thermal conductivity of the cylinder | Where 
material may vary with temperature. The steady-state 
heat flow equation in cylindrical coordinates is é) —ig vite dG(r: Fr $ 
NY 


1S’ 


a On’ 
la / Ov 0 [ Ov 
V -(kVv) (A(oy )4 — ( k(v) ~ 
ror ) 


or) az dz 


l ra) 


r? 00 00 


2 a ee 
(kiv) ou )\=— Wir) w=wr, Z, N=} W (r')G(r; rv’) dr’. (6) 


In (5), n’ is the outward drawn normal to the surface, 
Ee ee oy en ene a ey dS' is an element of area on the surface, and the 
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integration is taken over the entire surface of the 
cylinder. In (6), dz’ is an element of volume and the 
integration is taken over the whole volume of the 
cylinder. The quantity G(r: r’) is a time-independent 
Green’s function which we choose to express in the 
two alternative forms [2, 3, 4, 5. 6]: 

ae > > 


Thk* a 
m 


Gir: r’) 


x 1 


cos m(é— 8’) J nh Qmnl DF lh Gtnt b) 


= 
Omnd & 4 





1(Qmn) sinh (Qyn»l/b) 


: z 2 l—z, a 
sinh | Qin sinh (cnn — : (7) 
: b ) 


where J, is the ordinary Bessel function of the first 
kind and order m, Qn» is the nth positive root of 
J m(Qmn)=0. ze is the lesser of z’ or z, and zs is the 
ereater of z’ or z: 


Gir: r’) 


ae ) Ku ( ae 


nab nars\|1,(nar</l) 
— Km (“T ) In | a | T,(nmb/))° (8) 


where /,, and K,, are the modified Bessel functions 
of the first and second kind, respectively, and order 
m, r< is the lesser of r’ or r, and ry is the greater of 
r orr. 

The physical meaning of the two components of 
(4) is easily seen. If we set W(r')=0 for all r’, we see 
that « vanishes so that u=y. Thus y is the potential 
distribution due to the boundary conditions in the 
absence of internal heat generation. If we set u(r’)=0 
for all r’ on the surface of the cylinder, we see that 
y vanishes so that u=w. Thus w is the potential 
distribution due to the internal heat generation with 
the boundaries held at zero potential. 


Equations (7) and (8) are completely equivalent: 
however, in some cases one may be preferred over the 
other, either because the integrations indicated in (5) 
or (6) are easier or because the infinite series resulting 
from these integrations are more rapidly convergent. 
If the surface integration of (5) or the volume integra- 
_tion of (6) cannot be obtained analytically, y(r) or w(r) 
can be determined by numerical integration. 

Since we are considering the case of a homogeneous, 
isotropic material, W(r') would usually be expected 
to be axially symmetric for materials which self-heat 
by radio-active decay. In general #/(r') would not be 
axially symmetric for cylinders which are heated by 
absorption of energy from an external source, e.g.. 


irradiation by neutrons. For the specific case of 
axial symmetry, we assume the following boundary 
conditions: , 

u= fir) 


u= g(r) 
u=f(b)+|g(b)—f(b) | 7tAl2). 


where f(r) and gir) are arbitrary potential distributions 
on the flat surfaces of the disk, and where, in order to 


, Insure continuity, we shall require that the function 


h(z) vanishes at the flat surfaces of the disk. 
h(0)=h(1)=0. We can substitute 


so that 


these boundary 


conditions into (5) and obtain 


y(r, z)=f(b)+4 ib) — fb) ‘| 


PY Jo an i) 1,, sinh (cw, ; | + B,, sinh (a, See 


n 


+y Col) in (=). ae 
1 


: l l 
where a, =a, and the coefficients A, B,. and C;. are 
viven by 


) 
1), 


b2J3(a,) sinh (a,1/b) 


rb 
| { g(r) 
0 


2(b)}rJolanr/b)dr. 


») 


b?J3(a,) sinh (a,//b) 





{fir) —f(b) }rJolanr/b)dr, 


and 
») aw | — 

Cc; => | h(z) sin {| — ) dz. (13) 
Up(kab l) J0 \ l 
In obtaining (10) from (5), we elected to use (7) to rep- 
resent Gir: r’) on the flat surfaces of the cylinder and 
to use (8) on the convex surface. Equation (10) can also 
be written down without recourse to Green’s functions: 
this was done in [1]. 

In the specific case of axial symmetry, we can inte- 
grate over @ in (6) and obtain 


i b 
w(r, z)=27 [ 44 arr Wr, z Gr, z r',2'), aa 
{ 


J 0 ) 
where the Green’s function now takes the two alterna- 
tive forms: 
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JolQnr) b Jolayr'/b) 
ay J H(a,) sinh (a@,1/b) 


x, 


Twhh* — 
” 





< sinh (a, = sinh | @y 


b 


ta (5) Ko ("5 


nib 


j (16) 


Ky | 


(naz Iy(nzr </l) 
(2) a 


Iyinzb/1) 


2.1. Uniform Heat Generation 


For the special case where H (r)= Ho. a constant, 
substituting (15) into (14) and performing the indicated 
integrations yields 

2b7i 0 N J olant b) 
aa a3} (Qn) 
| cosh Ja ({- 2z)/2b | 


cosh (a,1/2b) 
Making use of the identity [7]. 


8 N J Qn) 


a) (Qn) 


eq (17) can be written in the form 


b2W 
1 


(1 


wir, Z) },2 
2b°W 


/ cm a3) (Qn) 


n l 


~S Jo(a@nrlb) cosh [a,(/ - 2z)/26| 


- - (19) 
cosh (@,1/2b) 


The alternative expression for w(r, z), derived from 


(16), is 


Io(nar/l) 
— 20 
Iy(nab 4 : ana 


where the notation n=1, 3, 5, under the summation 
sign indicates that the series is to be summed over 
odd values of n only. Making use of the identity 


— sin (n7rx)=x(1 — x), 21) 
A 


which can easily be derived, eq (10) can be written 
in the form 


L Io(n7r/l) 
— ———— sin 
, . n* Lo(nzb/1) 


Equations (19) and (22) could have been derived 
without recourse to Green’s functions for this specific 
case where H (r)=Wo since the leading terms in 
these two equations are particular solutions to Pois- 
son's equation for this case and are sufficiently simple 
to be written by inspection. In general, however, this 
could not be done. In principle. eq (6) provides us 
with a solution for wir) regardless of the particular 
form of #/(r) and hence is completely general. 

\ practical case in which VW (r) is essentially constant 
(and therefore (19) and (22) are applicable) is that of 
a sample which self-heats by radioactive decay and 
for which the mean free path of a fission fragment is 
much less than the dimensions of the sample. 


3. Calculation of Thermal Conductivity 


The total heat flow through a circle of radius a. where 
a= b, at z=0 is 


Q"=27 [ : rk(v) | _ dr 


fra r 
‘ ou 
ami? {0 (M) a 
0 OZ/ z=0 


From (4), this can be written as 


J0 OZ/ 2=0 


O”"=0+0". 


where 


a ae a 
oy Pe is. 
r ( ~~) dr = 7ra*k* (2) 
\OZ/ 2=0 0Z/ z=0 


Q'=27k* 


“ /aw Sy [OW 
r | | dr= tra7k* ’ 
0 OZ 2z=0 Oz z=0 


From (24), (25), and (26) we obtain 


« _Q"—ma*k* (dw/dz)z=0 
ma” {dy/dz)z=0 





(27) 


This equation gives the thermal conductivity in terms 
of measurable or calculable quantities. Q” is the meas- 
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ured total heat flow through the circle r<a at z:=0. 
The quantity (@y/dz)--o. which is the average tem- 
perature gradient which would exist over the circle 
r<a at z=0 due solely to the boundary conditions 
but with no internal heat generation, can be computed 
from the observed temperature distribution on the 
surface of the cylinder using (5) and (25). The quantity 
(dw/dz)--». which is the average temperature eradient 
which would exist over the circle r<a at z=0 due 
solely to the internal heat generation but with the 
boundaries maintained at zero potential, can be com- 
puted from the (assumed known) rate of heat genera- 
tion, HW (r), in the cylinder using (6) and (26). Note that 
u(r) is inversely proportional to /* so that the quantity 
k* (ou » on the right-hand side of (27) is actually 
independent of A 

If we again consider the specific case of axial sym- 
metry, (10) can be substituted into (25) to obtain 


Oz): 


e(b) ~ fib) 
l 


O=ra?k | 


1, —B,, cosh (a, 


>> Jan) | an 


n i 


ae > Colt | l 


1 
} 


iad |. (28) 


nd (13). 
eration, 


where 4,,. B,. and C;, are given by (11), (12). « 
In the case of constant internal heat ge: 


} Wy. and (26) becomes 


y(Qr) 


OO = 7a lh yO (29) 
where Q=Q(a/b: //b) has two forms corresponding to 
(19) and (22), respectively: 


1)* 


. Jila,a/b) 
— tanh | a, 
ned 1 (Ay) 


== GA 

’ 1 
T,(n7za/l) 
m1 y(nab/l) 


In the special case. //b > 0. (30) reduces t 


a(£: 0) Jana) 
b 


ij 
20 * 
= % 


a = ai Jian) 
where we have made use of the identity [8], 


JilQnx) 


4 


im, AJ (an) 


If //b — 0, (31) also reduces to Q= 1/2. 
In the special case, a/b > 0, the two forms of 2 are 


(9. £\ —26 
0(0; 7) a 


1 


« tanh (a,//25) 


az J (Qn) 


(34) 


a(o:7)=5-4 ¥ : (35) 


b ee n2fo(nab/ly 


For large values of //b, (30) approaches 


Q) 1] ~, 


where 


Jivlana b) 


«Fes ARS : —, 
bak) lan) 


a (a 


Values of ® are given in table 3 of [1]. There 
error in table 3 of [1 ]in the value tabulated for a/b = 
which should be 1.96630. 

Using a digital computer, numerical values for 2 
have been calculated for a range of values of a/b and 
//b and are given in table 1. Values were calculated 
using both (30) and (31) and the resultant two sets of 
values agreed to the number of significant figures 
given in table 1. Representative values are plotted 
in figure 1. 

As discussed in [1]. the thermal conductivity, 4* 
corresponds to the reference temperature, 7*, from 
which the variable, v, is measured. In principle, 7% 
may be arbitrarily selected and the mathematics will 
yield the corresponding value of &*. Thermal con- 
ductivity values may have to be computed by iteration 
as discussed in sections 4.2 and 5.2 of [1]. 


4. Example 


The above analysis is applicable to the several 
methods ef determining thermal conductivity which 


a/bsi.0 


2 
Lib 


FIGURE 1. The coefficient Q for the case of uniform heat generation. 


The values of Q 


shown 


corresponding to intermediate vales of a/b lie between the two curves 
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were discussed in [1] and also to other steady-state 
methods utilizing Dirichlet boundary conditions in 
which the temperatures at all points on the boundary 
are specified. In [1]. we discussed methods of thermal 
conductivity determination in which the investigators 
had assumed parabolic radial potential distributions 
on the flat surfaces of a cylinder and a linear longi- 
tudinal potential on the convex surface. This is a case 
having radial symmetry so that (10) can be used with 
the boundary conditions (9), in which f(r), g(r), and 
h(z) are given by 


and h(z)=0 


where ¥» and Y, are the potentials (corresponding to y) 
at the centers of the faces of the disk, and Ey and E, 
are the potential differences between the edge and the 
center of the disk at z=0 and z=/, respectively. For 
these boundary conditions, (28) becomes 


TABLE l. 


50000 
50000 
50000 
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50000 
50000 
50000 
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zotll Rs 
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.17794 
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.11884 
.10696 


.17573 
.15081 
13201 
11736 
-10562 


.17292 
.14840 
.12990 
.11548 
.10393 


.16890 
.14495 
.12687 
.11279 


.10663 10151 
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S J i(a,a/b) 


VY,=1-16 (41) 


are factors which depend only on the geometry of the 
cylinder and on the fraction of the disk from which the 
heat flow is measured. Numerical values of Vo and WV, 
are given in tables 1 and 2 of [1] for a range of values 


of a/b and L/b. 


In [1] we did not discuss the effect of internal heat 
veneration. For the special case of parabolic radial 
potential distributions on the flat surface of the cylin- 
der, a linear longitudinal potential on the convex 
surface, and uniform heat generation throughout the 
cylinder, (27) becomes 


The coefficient 1 as a function of \/b and a/b for uniform heat generation 
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.21930 
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.18668 
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14777 
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.15559 
.14394 
.13384 


.23153 
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.19505 
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.16788 
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.13406 
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.11530 
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ae OQ" — 7a:lWO 
(7ra?/l )| y; —YVo+ E,V, — EoVol 





(4.2) 


As a specific example of the use of (42), we consider 
the method described by Hoch et al. [9], in which the 
heat flux from the flat surfaces of the disk is calculated 
from the temperature and the total hemispherical 
emittance of the specimen surface. using the Stefan- 
Boltzmann radiation law. That is. the heat flow, Q”, is 
given by 

Q" = za*eoT* (43) 
where € is the total hemispherical emittance, o is the 
Stefan-Boltzmann constant. and 7" is the average value 
of T* on the surface over the circle of radius a. Sub- 
stitution of (43) into (42) results in 


NeoT! —1WQ] 
} _—. Yo T EY, — EyVo ; 4) 


With the simplified boundary conditions assumed by 
Hoch et al., namely ¥;= Yo and £,;=£o,. this reduces 
to 


eo! — IW Q| 


ore (45) 
EI Wo | 

If this is evaluated at the center of the disk (i.e... a= 0), 
as was done by Hoch et al., the thermal conductivity 
is given by 


jp» —lleoh! — WV 
E(¥; — Wo! 





(46) 


where 7, is the (absolute) temperature at the center of 
the disk at z=0, and Vy. V; and 2 are to be evaluated 
at a/b=0. If there is no heat generation (i.e., Wo) =0), 
equation (46) reduces to a form which is exactly equiv- 
alent to eq (13) of Hoch et al. [9]. 

We postulate the following data as having been 
acquired experimentally: 


b=1.00 cm /=1.80 cm 
e—O.17 Wo = 14.3 W/em* 
To = 2360 °K Eo 28 deg. 


From tables 1 and 2 of [1| for the case a/b=0.0 and 
//b=1.8, we obtain 


YW, 0.87364 Wy =e 2.85220. 


From table 1 of the present paper, we obtain 
Q = 0.28749. 


Using 0 =5.6697 X 10-'* W/em? deg*. substitution of 
the above values into (46) yields 


0.388 W/cm deg 


as the thermal conductivity of the sample at the tem- 
perature corresponding to v=0. 
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Equipment and techniques for acoustic cavitation work are described. The test liquid and its 
container form part of a self-sustaining programmable oscillator. Emphasis is laid on the treatment 
of the liquid necessary to yield reproducible results after artificial nucleation, especially by neutrons, 
a-recoils, and fission. With neutrons, for instance, the cavitation rate rises rapidly with acoustic (nega- 
tive) pressure, and at fixed pressure is proportional to neutron flux. The cavitation events are random, 
and no appreciable induction or decay times are cbserved. Each cavitation arises from the action of 
a single neutron (or @-recoil, or fission). The cavitation “threshold” (roughly negative pressure below 
which cavitation is rare) is highly variable from liquid to liquid. Thresholds at about room tempera- 
ture for liquids irradiated with 10 MeV neutrons ranged from 1 to 2 bars for n-pentane, ether, and 


freon 113 to over 50 bars for water. 


Key Words: 


Acoustic cavitation, alpha-particle-induced cavitation, cavitation, cavitation nucleus, 


cavitation threshold, fission-induced cavitation, neutron-induced cavitation, radiation- 


induced cavitation, threshold. 


1. Introduction 


Our interest in cavitation research was aroused 
partly by the work of Lieberman [1959] but mostly 
by that of Sette and Wanderlingh [1962]. Also we felt 
that it would be desirable to develop techniques ca- 
pable of yielding results more reproducible than those 
in the literature. Of course most investigators were 
aware that the day-to-day variability of their results 
reflected that of the population of nuclei in the liquid 
under test:' indeed a cavitation experiment may be 
designed so that the results are best interpreted as a 
determination of the distribution of nuclei, a point 
emphasized by Messino, et al. [1963, 1967]. Still, we 
felt that some improvements in both apparatus and 
technique would be of value, and in this paper we 
describe what we have done so far and give some 
sample results which, although as yet uncoordinated 
by theory, illustrate how data of potential importance 
can be collected by relatively straightforward means. 

We revert to the work of Lieberman [1959]. He was 
able to induce cavitation in acetone and pentane?” by 
neutron bombardment. The results, although very 
interesting and so far as they go in substantial agree- 
ment with those of the present work, were hardly 
surprising in view of what was then known about 
bubble chambers. The situation is different in the 
case of Sette and his co-workers [Sette and Wander- 
lingh, 1962: Messino, et al., 1963] and later, of Finch 
[1964] as well. Their results imply the existence, in 
water, of neutron-generated cavitation nuclei which 
have a long half-life, sometimes an hour or more. 


*Work supported in part by the U.S. Office of Naval Research. 

'The disparity of results of various investigators reflects, in addition, the differences in 
experimental techniques and conditions, and in definitions of “threshold.” 

Probably n-pentane 


The apparatus here described is suitable for de- 
termination of the following: 

(1) Cavitation “thresholds” 
arbitrary definitions. 

(2) Cavitation-event 
sure above threshold. 

(3) Effect of static pressure. 

(4) Influence of various nuclei, 
ficial. The natural nuclei can be filtered out and 
replaced by others. Among these we emphasize 
those which are radiation-induced; by neutrons, 
by disintegration of dissolved a-emitters, and by 
fission, both induced and spontaneous. In the case 
of neutrons, the prevailing view [Lieberman, 1959] 
is that the nucleation is connected with the deposi- 
tion, in the liquid, of the kinetic energy acquired 
by an atomic nucleus of one of the heavier elements 
in the liquid as a result of an elastic collision with 
a neutron. In the case of a@- disintegration, the heavy 
recoil nucleus (a negative ion rather than a positive 
ion as in the previous case), is responsible, accord- 
ing to Riepe and Hahn |{1961|, who have examined 
this process in a bubble chamber. 

So far our apparatus is not well adapted to studies 
of the effect of frequency. 

For an overall view of acoustic cavitation the reader 
is referred to Flynn [1964] whose article contains 152 
references. 


according to various 


rate as a function of pres- 


natural or arti- 


2. Apparatus, Major 


Other things being equal, cavitation occurs in that 
region of the liquid where the (negative) pressure 
is greatest; therefore it is necessary that the pres- 
sure at the walls of the container be substantially 
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less than that in the “test” region, else cavitation 
will occur preferentially on the walls (or on the trans- 
ducer, if it be immersed), and the results will not 
then represent properties of the liquid alone.* 

In the case of traveling waves, they may be focused. 
Sette and: Wanderlingh [1962] used reflectors. The 
radiation from the concave side of a transducer 
which is the segment of a sphere is focused at the 
center: this arrangement was used by Willard [1953], 
and on a massive scale, by Sirotyuk [1962]. 

Standing-wave systems for the most part exploit 
the radial modes of a liquid-filled spherical flask, 
an arrangement devised by Leonard [1946] for another 
purpose and first adapted to cavitation work by Gallo- 
way [1954]. In such resonant systems the power 
requirement is just that needed to supply the losses; 
this advantage over traveling-wave systems is es- 
pecially important where high sound pressures are 
required. 

In our own work we use resonant systems, either 
cylinders or conical flasks. 


2.1. Cylinders 


The modes of vibration of a cylinder of fluid are well 
understood; here we give a brief treatment just to 
fix ideas. Modes that have an azimuthal dependence 
give zero pressure on the axis and maximum pressure 
on the walls; these are undesirable for our purpose 
and we consider radially symmetric modes only. The 
velocity potential, ignoring losses and with the time 
dependence exp (iwt) suppressed, is 


b= Jol krr) [A cos k.z+B sin kz] ‘ (1) 
with 

k2=e2/2=k2 +k: . (2) 
the angular frequencies, w, are thus known in terms 
of the radial wave numbers. k,. and the axial wave 
numbers, k;, and the plane-wave speed of sound, c, 
in the liquid. The ratio B/A depends on the impedance 
at the flat ends of the cylinder. Our practice is to use 
a thin sheet for the ends to make the impedance nearly 
zero; then with the origin at one end, and because the 
pressure is proportional to @¢, 4=0 and kz=n7/L, 
L being the length of the cylinder and n an integer. 
Equations (1) and (2) reduce to 


b@=BJo(k-r) sin (nzz/L). (3) 
and 


= k2 + n?777/L?. (4) 


The pressure, radial velocity, and axial velocity are 
given by 


p=—iwpd, 


uy = 0o/dr=— Bk, J (kr) sin (n7z/L), 


Also, the apparatus may be damaged 


and 
uz=‘db/dz= (Bnz/L)Jo(krr) cos (na/L), (7) 


respectively, with p the undisturbed density of the 
liquid. 

At r=a, the inside radius of the cylinder, the 
ratio p/u, must equal the radial driving-point im- 
pedance, Z(w), at the wall. Then, with k,a= x, 
we get from eq (5), (3), and (6) 


Z\(w). (8) 


iwmpasy(x)/ (x) Ji (x) 
and from eq (4) 


(w/c)? = (x/a)* + (n7/L)?. (9) 
Equations (8) and (9) determine simultaneous values 
of m and x, Z(w) being supposed known.* 

A difficulty common to all resonator systems is 
that the container must be pierced and provided with 
filling tubes for introduction of the sample. These 
upset the sound field and also introduce losses. The 
above analysis shows that the tubes should be lo- 
cated either in the lateral walls near the ends or in 
the ends near the circumference. For at the positions 
r=a, z=0 or L, both the pressure, p, and the radial 
velocity, u,, vanish [eq (5), (3), (6)|. As for the axial 
velocity, uz. the ratio of its value at these positions 
to the maximum value is just Jo(x)leq (7)| which is 
also the ratio of the pressure on the lateral wall to 
that on the axis for any z. To prevent cavitation from 
occurring on the wall, Jo(x) must be small—just how 
small depends on the smoothness of the wall. A 
“good” mode from this point of view therefore implies 
a relatively small disturbance at the filling tubes. 
Such modes are best found by trial; calculations from 
eq (8) and (9) are possible but tedious. 


In the case of a spherical resonator all positions of 
the neck are equally bad. The experiences of Finch 
[1964] illustrate the point. These considerations, 
among others, impelled us to use cylinders instead 
of spheres. We will point out other advantages, and 
many disadvantages, as they come up. Just now we 
remark only that good geometry is easier to come by 
in a hollow cylinder than in a hollow sphere. 


a. Integral Drive 

Cylinders of piezoelectric ceramic are manutac- 
tured to close tolerances in a variety of sizes. We 
have used mostly 7.5 em OD by 7.5 to 10 cm high and 
expect to try larger ones. As shown in figure 1, the 
cylinder serves as part of the container, as the driving 
transducer, and as a pickup. It can be of lead-zirconate 
titanate or of barium titanate; in either case the mate- 
rial should be one especially formulated for low loss 
(internal heating) at high drive. 


‘We here ignore the dependence of Z(@) on z 

From these simple equations we have calculated the resonant frequencies of liquid 
filled cylinders and gotten good agreement with experiment. The pressure distributions, 
measured with a small probe, agree with those calculated, sometimes within a few percent 
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FIGURE 1. Self-driven cylinder, schematic cross section. 


(a) and (b), Alternate terminations 
uration on (b) is the 


and placement of filling tubes. The electrode config 
same as on (a 
1. Ceramic, either lead-zirconate titanate or barium titanate 
2. Metal rings, stainless steel or brass 
3. Sheet material. stainless steel or glass 
Cemented joints 
Filling tubes 
Electrode, ground 
Electrode and connection, input 
Electrode and connection, output 


The metal rings allow for filling holes and connec- 
tions. The flat ends are either stainless steel, 12 to 75 
wm thick, or glass, 250 wm thick. The glass, which is 
cut from microscope cover-glass stock, is fragile and 
difficult to handle, but it allows one to see what goes 
on. The thinner material provides a better pressure 
release than the thicker, but it is more easily excited 
(parametrically, we think) into radial flexural vibrations. 

The joints shown in figure 1 are best cemented with 
a silicone-rubber adhesive which cures at room tem- 
perature. Many organic liquids attack this rubber, 
however; for these, some other adhesive, such as an 
epoxy, may serve. We have found no universal adhesive 
nor one which lasts indefinitely with any single liquid. 

Manufacturers supply cylinders electroded on both 
the inside and outside with a fired-on silver-glass 
mixture. The outside electrode configuration shown 
in figure 1 is produced by dividing the electrode as 
furnished and the inner ground electrode is connected 
with the outer ground bands by a conducting paint. 
In some cases we have removed the inner electrode 
and replaced and extended it over the edges with 
fired-on gold or platinum.® 

These cylinders have two main limitations. The 
first, namely limited life of the adhesives, we have 
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already mentioned. The second has to do with the 
roughness and friability of the inside walls. We would 
prefer to operate at the “breathing” mode [(r, 6, z) 
=(1, 0, 1)] because the cylinder walls are for this 
mode more nearly pressure-release than for any other. 
Also, a larger fraction of the liquid is contained in the 
high-pressure region. However, unless the pressures 
are low, the wall motion in the (1, 0, 1)-mode is so large 
that the inside electrode (and sometimes the inner 
surface of the ceramic) disintegrates and contaminates 
the liquid. The overall performance is better in the 
(3, 0, 1)-mode and higher pressures are obtained. 
There is a tendency for cavitation to occur on rough 
spots on the wall, and this limits the maximum useful 
negative pressure obtainable to the range 20 to 30 
bars, depending on the particular cylinder and the 
particular liquid. Such pressures are more than ade- 
quate for many organic liquids exposed to fast 
neutrons. 


b. External Drive 
“Clean” water exposed to fast neutrons cannot 
be cavitated in the ceramic cylinders described above. 
\ special cylindrical vessel was fabricated from stain- 
less steel tubing. The schematic, figure 2, is mostly 
self explanatory. 

The cylinder is operated in the (3, 0, 1)—mode, 
and the driver is adjusted to resonate at the same 
frequency at the temperature of test (in this case 
about 42.9 kHz at about 30 °C) in order to minimize 
the stress in the epoxy joint between tiie driver and 
a small milled flat on the cylinder. 


* The cylinder is then repoled. 
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FIGURE 2. Externally driven cylinder, schematic cross section. 

1. Stainless steel cylinder, 8.9 em O.D. by 10.2 cm high, wall 0.16 em thick. Inner surfaces 
polished 

2. Stainless steel sheet, 75 wm thick. 

3. Silicone-rubber adhesive joints. 

4. Filling tubes. 
5. Composite, prestressed driving transducer, 13 mm O.D. The center section is lead 
zirconate-titanate: the ends and prestressing screw are of brass. 

6. Wire strain gage. 

7. Epoxy joints. 
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FIGURE 3. Conical flask. 
1. Erlenmeyer flask. neck removed 
2. Solid glass plug with 2 holes 
3. Filling tubes 
4. Ceramic transducer 
>». Epoxy joint 


At some temperature which fluctuates from day to 
day but is always within the range 29-3] °C, the over- 
all Q is a maximum at about 12 000. Negative pres- 
sures of 160 bars are sustained by the water for a 
minute or so at a time, and of 210 bars for a few sec- 
onds. The cavitations are presumed to be nucleated 
by either residual motes or adventitious neutrons. 

The readout is from a wire strain gage noninductively 
wound around the circumference of the cylinder near 
midheight. It consists of three double turns of 75 uwm- 
diameter Karma wire cemented into a double V- 
thread, 50 wm deep, cut into the cylinder. 

This cylinder could be used with various liquids 
and at various temperatures for work at low pres- 
sures; however, we have used it only for high pres- 
sures in water at one temperature. A change in tem- 
perature detunes the system and requires the trans- 
ducer to be retuned. So far we have not found it 
worthwhile to do this. 


2.2. Conical Flasks 


Cylindrical geometry allows calculation of the pres- 
sure field and experimental verification (see footnote 5); 
this was one motivation for our choice. Nevertheless, 
we have made considerable use of the arrangement 
shown schematically in figure 3 for which these 
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FIGURE 4. Circuit for use with cavitation cell. 


Schematic block diagram lhe cell is the fre quency-determining element of an osc illator. 


advantages are absent. The container is made by re- 
placing the neck of an ordinary Erlenmeyer flask 
(125. 250. 300. 500, or 1000 ml) with a massive glass 
plug, typically about 25 mm OD by 25 mm long. 

It was found empirically that each assembly has 
several modes for each of which there is a pressure 
maximum on axis and almost no motion in the plug. 
The flask can therefore be held by the plug without 
much loss of Q; the Q may reach 10 000. The resonant 
frequencies and the Q’s depend on the particular 
flask, transducer, liquid, and on the temperature, 
and cannot be predicted. 

The advantages here are (a) that the vessel is of 
one-piece, all glass construction with no joints and 
(b) that the assemblies are cheap once the glass- 
blower has taught himself how to eliminate rough 
spots on the inside. These advantages often outweigh 
the rather haphazard nature of the entire scheme, 
and we have collected a considerable amount of 
valuable data using it. 

The drive is the voltage on, and the readout the cur- 
rent into, the transducer. See section 2.3. 


2.3. Electronics 


In all of our setups the cavitation cell, together with 
its associated reactive elements, such as transducers, 
is the frequency determining element of a self-sustain- 
ing oscillator. The scheme is shown in figure 4. The 
wave filter, a simple combination of L,C-elements, 
is adjusted so that its pass band is roughly centered 
on the frequency of interest. For some setting of the 
phase shifter the system will oscillate at the desired 
frequency if the loop gain is large enough. Small 
changes in the resonance frequency of the cavitation 
cell, caused by temperature drifts, etc., are thus 
automatically tracked. 

The amplitude of oscillation is adjusted to any de- 
sired value by variation of the supply voltage to a 
transistor stage (“level set” in fig. 4). The level can 
be set manually or programed in any desired way. 
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FIGURE 5. Circuit for use with cavitation cell. 
Variant featuring easy portability 


For instance, the pressure in the cell can be caused 
to increase linearly with time by application of a linear 
voltage ramp to the “level set.” 

The “cavitation detector” rectifies and differen- 
tiates the amplified output of the cell. A cavitation 
event drops the Q of the system and therefore also 
the output amplitude; the resulting pulse can be made 
to trip a relay and interrupt the supply to the “level 
set,” stopping the oscillation. The cavitation detector 
has adjustable sensitivity and can be made to trigger 
on mild events or on catastrophic events only. It can 
also be disabled. The rectified signal is amplified and 
applied to a small loudspeaker so that the cavitation 
events can be “heard.” 

After the relay has been tripped it can be reset 
after any desired elapsed time which exceeds the 
minimum delay set by the Q of the system. Delays of 
50 sec or less are obtained from the sweep voltage 
available from the oscilloscope; longer delays are 
provided by a clock motor-switch arrangement. 

The arrangement shown in figure 5 is not so versatile 
as the first, but is very simple and therefore well 
adapted to portable equipment. The amplified output 
synchronizes a Wien-bridge type oscillator. The os- 
cillator functions as amplifier, filter, phase-shifter, 
and amplitude control (level set). The main disadvan- 
tage is that each time the amplitude is changed the 
frequency dial must be changed to restore the proper 
phase shift, for at a given setting of the frequency 
dial, the phase shift between the synchronizing and 
output’ voltages depends on the amplitude of the 
former [Owen, 1967]. A limiter at the synchronizing 
jack would perhaps improve the performance. 

In all of the above, the term “cell output” refers 
to an electrical signal proportional to the maximum 
sound pressure in the liquid. In figure 1, the output 
is the voltage on the receiving bands (8). In figure 2, 
it is the voltage across the strain gage (6) which has 
a resistance of 510 ohms, carries a bias of about 14 mA, 
and is coupled to an amplifier through a 1 to 14 step-up 
transformer. The signal-to-noise ratio is poor so that 
filtering is necessary. In this case the amplifier is 
stagger-tuned; the 3-dB bandwidth is about 2 kHz. 

For the case of figure 3, the input voltage is applied 
to the electrodes of the transducer; the output is the 
current into the transducer. As shown in figure 6, 
the current is sensed with a “‘clip-on probe,” a common 
accessory for an oscilloscope. The probe measures 
the algebraic sum of the current to the transducer 
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FIGURE 6. Output system for conical flask. 


Can also be used with the cylindrical cells of figures 1 and 2. Especially convenient for 
remote operation. The capacitance of (3) is made equal to that of the leads plus that of the 
transducer: the probe then responds to the motional current only. 


and that to the capacitor (3) with its leads reversed. 
The capacitance of the auxiliary capacitor is made 
equal to the sum of the “free” electrostatic capaci- 
tance of the transducer and the capacitance of the 
leads: only the motional part of the current is meas- 
ured. This arrangement makes it possible to have the 
cell a hundred meters or more from the rest of the 
equipment with only a single coaxial cable connection, 
provided that the power amplifier is capable of deliver- 
ing the necessary in-quadrature current. Where a 
radiation hazard exists, this is a valuable feature. 


2.4. Calibration 


a. Probes 

A principal object of a typical experiment is the 
determination of the maximum negative pressure in 
the liquid at the onset of cavitation. Those who work 
with spherical flasks usually make use of a probe 
microphone (probe). The active element is in most 
cases a small ceramic tube or disk; it is mounted on 
the end of a metal tube of diameter 1.5 to 3 mm which 
protrudes through the neck of the flask so that the 
active element is in a subsidiary pressure maximum. 
The designs are numerous: for a complete descrip- 
tion of an especially clever one, see Barger [1964]. 

We have built many such probes, a few novel but 
none remarkable. In our work we find it essential 
not to have the probe in the liquid during test. In- 
stead, the output of the cylinder is compared with 
the reading of a calibrated probe which passes through 
a small hole in the flat top of the cylinder during a 
separate run during which the sound pressures are 
not so great as to produce cavitation on the probe. 
Such a comparison has to be made for each liquid 
at each frequency and at each temperature. The pro- 
cedure is tedious and in any case cannot be used in 
conical flasks. Also tedious is the calibration of the 
probe itself. For instance Barger [1964] employed 
a free-field reciprocity method, in water. We shall 
describe a novel method later on in this section. 
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b. Overpressure 

We describe a method of calibrating a cell, that is, 
of relating the sound pressure in the liquid to the 
reading of the output device, which is obvious but 
subject to some limitations. Suppose that under some 
“specified condition a cavitation occurs at some nega- 
tive pressure, T (threshold, a posiiive number); then 
if ps is the static, and p, the amplitude of the acoustic, 
pressure, 


T= pa Ps- (10) 


The output is some voltage V, 
to Pa, SO that eq (10) is 


supposed proportional 


V=(ps+T)/m, (11) 
in which pa=mV. Two measurements of V for two 
values of ps determine T and m. The method fails if 
T depends on py. For instance, if the cavitation nucleus 
is a gas or vapor bubble, somehow stabilized at (say) 
atmospheric pressure against dissolution or condensa- 
tion, respectively, then an increase in the static pres- 
sure ps will shrink the bubble and the liquid will be 
“strengthened,” i.e., T will change with ps, so that 
eq (10) is no longer linear. This phenomenon has been 
discussed by Strasberg [1959]. 

On the other hand, we find that if the cavitation is 
induced by neutron bombardment or a-disintegration, 
then the threshold, 7, is independent of the static 
pressure ps. The cell is enclosed in a pressure vessel 
and the output, V, at threshold is determined for 
various values of ps. A plot of V as ordinate versus p, 
as abscissa yields a straight line, of which the p-- 
intercept, with reversed sign, is the threshold, T. There 
is no measureable hysteresis. In practice, T is calcu- 
lated from the data by least squares. The vessel is 
pressurized with nitrogen and the value of ps read on 
a high-quality bourdon-tube gage. A piston arrange- 
ment (actually a hypodermic syringe) transmits the 
pressure to the liquid, as will be described later. 

A major advantage of this “‘over-pressure” method 
is that the nature of the pressure field in the cell does 
not need to be known. 


c. Calibration of Probes 


As a byproduct, the overpressure method provides a 
simple way of calibrating the small, rather insensitive 
probes already discussed. 

Some threshold, T, which as always has an arbitrary 
definition, but which must be reproducible, is meas- 
ured by the overpressure method in a cylindrical cell 
such as shown in .figure 1. The maximum acoustic 
pressure amplitude (at the center of the cell) is now 
known in terms of the output voltage. A probe inserted 
through a small hole in the top is then easily cali- 
brated against the output voltage; the result applies 
at the frequency of test only. Other frequencies may 
be gotten by changes of mode and/or liquid. For probes 
which are sensitive on the flat end only, the response 


may be extended over any frequency range by means 
of an electrostatic actuator.’ 


3.. Apparatus, Auxiliary 


3.1. Cleanup of Liquid 


The exact nature of the natural nuclei in liquids 
is unknown and almost certainly variable from sample 
to sample and: from time to time. Presumably they 
are either bubbles (gas or vapor) somehow. stabi- 
lized upon motes or motes which are poorly wetted; 
in any case those larger than a specified size can be 
removed by filtration. Membrane filters of large 
area, with uniform pores and relatively large flow 
rates are long since a common article of commerce. 
General purpose types have a pore size as small as 
0.01 xm. These are attacked by some organic liquids: 
for many of these, other types having a pore size as 
small as 0.2 wm are available. 

The degree of cleanup required depends on the 
nature of the case. For instance, suppose one wishes 
to measure some cavitation threshold, 7|, under ir- 
radiation by neutrons from a_ particular source. 
Then it is necessary, in order to get reproducible 
results, to remove from the liquid all nuclei upon 
which cavitation can occur at than 
about 2'7,. Because the rate of flow of a liquid through 
a filter is smaller the smaller the pore size, a cleanup 
more drastic than necessary is undesirable. 

In these respects, so far at least as 7, is concerned, 
the most refractory liquid, among those we have so 
far tried, is water. If all motes retained by a filter of 
pore size 0.2 ym are removed, water will sustain 
negative pressures in excess of 200 bar for a short 
time, and of 100 bar almost indefinitely, so that 
conditions are satisfactory for the determination of 
a value of 7, near 50 bar. However, it does not suffice 
merely to pour distilled water through a filter into 
the test container, for the latter contains dirty air; 
even if it is first evacuated, the water will become 
contaminated by material on the walls. Rather it is 
necessary to make use of a scheme equivalent to 
that shown in figure 7. After the assembly shown is 
filled, the valves (2) may be set to form a closed sys- 
tem so that the water can be continuously circulated 
by a bellows or diaphragm pump. Alternatively, 
if a still is available which can produce water as 
fast as it flows through the filter, the pump can be 
dispensed with and the valves (2) set to admit and 
discard the water continuously. The cleanup time is 
greatly shortened by acoustic agitation of the cell 
during the procedure and by the addition of detergent 
to the inflowing water during the early stages. The 
progress of the cleanup can be observed through 
a dark-field microscope focused on the effluent 
in a piece of precision-bore rectangular glass tubing, 


stresses less 


7 The electrostatic actuator of Ballantine [1932] is described by Beranek [1949], p. 173. 
In our case it is not feasible to use it as an absolute device 
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FIGURE 7. Cleanup system, schematic. 
>of the components are used « 
valves 


nly in special cases 


optically finished on the outside. The microscope 
arrangement is rather crude, but it works well so 
long as the effluent is still rather dirty; however the 
field of view is so small that after a few minutes of 
operation the waiting time for passage of a mote be- 
comes long enough to fatigue the observer and render 
his perception unreliable. Its main virtue is that 
it makes it possible to determine whether or not the 
shock wave caused by the collapse of a cavitation 
bubble has dislodged additional motes from the 
wall. In the case of water in a stainless steel cell, 
such dislodgment happens almost invariably, so that 
meaningful observations can be taken only at long 
intervals, say 5 to 15 min. 

A much less stringent procedure is required for 
organic liquids, at least the dozen or so we have so 
far used. Two possible reasons come to mind; first, 
the thresholds, 7,, are lower than for water and sec- 
ond, the motes are either better wetted by, or more 
soluble in the organics. Whatever the reasons, it turns 
out that 7; is. in all cases so far, easily measured after 
only 10 or so minutes of circulation in an apparatus 
like that in figure 7. (The microscope is omitted.) Valves 
(1) are then closed and the connecting tubes detached 
so that the cell and relief element can be removed for 
use. The relief element, which js a hypodermic sy- 
ringe, maintains the ambient pressure, or the over- 
pressure, on the contents. 

For some liquids pouring or pumping through a 
filter into the flask, previously cleaned and rinsed 
with the liquid, is sufficient. In others, even the filter- 
ing may be omitted, especially if the sample is from a 
freshly opened bottle. 


3.2. Gas Content of Liquid 


Under some circumstances the cavitation threshold 
is sensitive to the gas tension in the liquid. This ten- 
sion can be adjusted by allowing the liquid to come into 





FIGURE 8. Apparatus for determination of oxygen, schematic. 
Sensor, active tip in liquid. 
Metering syringe, 5 cm* Luerlok hypodermic. 
. Rack, hinged to plunger, engages pinion (not shown) on clock motor (not shown). 
. Sample or standard syringe, 10 em* Luerlok hypodermic 
Stopcock, one-way, attached to syringe (4). 
. Stopcock, 3-way, attached to syringe (2). 


Waste 


) 


equilibrium with a gas at a partial pressure which is 
itself adjusted. For example, a liquid is easily satu- 
rated by bubbling the gas through it at the particular 
temperature and pressure for which saturation is 
desired. Also most of the dissolved gas can readily 
be removed by pumping on.the space above the liquid 
in a closed container while stirring violently. A com- 
bination magnetic stirrer and hot plate is invaluable 
for this purpose. 

The gas tension may be measured by any of several 
methods known to physical ¢chemists. Strasberg [1959], 
and others since, made use of a Van Slyke blood-gas 
apparatus, available from chemical supply houses. If 
the gas is air, which is usually the case,® one can use 
techniques peculiar to oxygen. We have found a polar- 
ographic technique based on the oxygen cathode to 
be very convenient. For a comprehensive review see 
Davies [1962]. We describe only the modifications 
peculiar to our work. 

In figure 8, (1) represents a sensor, about 1 mm in 
diameter, which is embedded in the wall of a tube 
about 3 mm in diameter past which the test liquid flows. 
The seasor has a current output (a few microamperes) 
which is proportional to the rate at which O. is reduced 
on an inert (gold, in our case) cathode separated from 


the test liquid by a thin film of electrolyte and a teflon 


* An exception is the study of sonoluminescence which is sensitive to the nature of the 
dissolved gas. 
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membrane. The QO» reaches the cathode by diffusion 
through the membrane so that the output current is 
proportional to the concentration of O. just outside it. 
The electrolysis depletes the liquid of O2, so that the 
current output depends on the speed of flow in the 
tube. The device shown in figure 8 standardizes the 
flow rate and allows compensation for other variables 
to be described. But first we note that the sensor, 
which is a modification of the “Clark cell’ is easily 
made. Ours is a modification of a commercial model. 
The original Ag anode has been replaced by Ag—Ago0, 
and instead of a KCI gel.we use, as electrolyte, 0.5N 
KOH saturated with Ag.O and thickened with hydroxy- 
ethyl cellulose. This arrangement is more stable and 
has a longer life than the original; it is based on the 
work of Connelly and of Carrit and Kanwisher [see 
Davies, 1962, p. 172]. 

The metering syringe in figure 8 is driven by a clock 
motor through a rack and pinion; the flow is 5 cm?/min 
which amounts to about 2 cm/s past the sensor. Be- 
fore making a measurement, a sample is withdrawn 
from the apparatus into a 10-cm*® syringe provided 
with a stopcock (5), taking the usual precautions to 
exclude air. At the same time a sample is withdrawn 
into another, similar syringe from a bottle of the same 
liquid in the oven and thus at the same temperature 
as the test sample. This “standardizing” liquid is 
maintained saturated with air at the temperature of 
the oven and the pressure of the atmosphere. In each 
case the first withdrawal is discarded so that the 
syringe is purged of air and residual liquid. The two 
syringes, “standard” and “sample” with stopcocks 
(5) closed are then allowed to come into temperature 
equilibrium with the ambient at the location of the 
apparatus in figure 8; the contents are then succes- 
sively run through. The ratio of the “sample” reading 
to the “standard” reading gives the oxygen tension 
in percent relative to saturation at 1 atm total air 
pressure. Before taking a reading, the metering syringe 
and associated tubing are purged with “test” or 
“standard” liquid by an obvious manipulation of the 
stopcocks; the rack is hinged so that it can be lifted 
from the pinion to allow hand operation which facil- 
itates the process. The output current reaches a steady 
value after about 1 cm? of liquid has been expelled. 

It will be seen that the scheme described eliminates 
all of the important sources of error. These are (a) 
changes due to changes in flow rate past sensor, 
(b) drifts due to aging of the cell, and (c) changes in 
sensitivity which arise from the high temperature 
coefhicients of diffusion in the liquid and in the Teflon. 


2:3: Radiation Sources 


As was mentioned in the introduction, cavitation has 
been induced by neutron bombardment, by a-disinte- 


gration, and by fission. 
a. Neutrons 


Most of the work so far done has been with (a, n) 
9. 
sources. In all cases the target atom has been Be’: 


the a emitters have been Ra [Sette and Wanderlingh, 
1962], Po?! [Lieberman, 1959], Am**!, and Pu? 
[Hahn and Peacock, 1963]. The main disadvantage of a 
Ra source is its high y-ray background, about 6500 
y’s/neutron. A Po. source is small, but has a short 
half-life (138 days). Of those sourcés we have had an 
opportunity to use, the best from the standpoint of 
size, half-life, and safety is Am-Be. The Pu-Be source 
is about as good but larger; we now use Pu-Be because 
of its lower price. 

The characteristics of (a, n) sources are listed in 
NBS Handbook 72. The neutron energy spectra of such 
sources are variable even among those nominally the 
same, and they are known only approximately. The 
energy varies from zero to a maximum of 13.1, 10.9, 
11.4, and 10.7 MeV for Ra, Po, Am, and Pu, respec- 
tively, and the average energy is about 4.5 MeV. 

Most of our work has been done with a 10-curie 
(160g Pu) Pu-Be source having a total output of 
1.47 X 107 n/s. 

The data would be more easily interpreted if the 
more nearly monoergic neutrons produced by particle 
accelerators were used. Finch [1964] worked with 
14 MeV neutrons produced by the [T(d, @)n] reaction. 
We have made a few runs using this and also the 
[Did, n)He*] reaction (2.5 MeV neutrons), but the 
neutron fluxes were rather unsteady. 


b. a-Disintegration 


The short range of the heavy recoil nuclei from 
disintegration of a dissolved a-emitter is a valuable 
property; the cavitation data are more easily inter- 
preted. In order that the dissolved emitter not affect the 
properties of the liquid too much, it should have a 
high activity, so that only traces are required. Por? 
perhaps represents a good compromise between high 
activity and reasonable half-life: we have not yet tried 
it for lack of facilities for safe handling. A possible 
difficulty is that for some liquids too much of the Po 
compound may be adsorbed on the walls of the con- 
tainer, leaving too little in the bulk of the liquid. 

Ordinary thorium nitrate, Th(NOs); -4H2O is 
soluble in the lower alcohols and ethyl ether. We have 
used concentrations of 1 to 10 g(Th) per liter in methyl, 
ethyl and isopropyl alcohol with good results. Urany] 
nitrate UO.(NOs;).-6H2O, is also soluble in these 
liquids. Although about three times as active as tho- 
rium, uranium is not as well suited to this work because 
of the high spontaneous fission rate ® of U?**. 

It should be remembered that laboratory reagent 
grades of thorium and uranium salts are contami- 
nated with daughter products, many of which are 
a-emitters. 


c. Fission 


Cavitation occurs at relatively low stresses on 
fission events in liquids. Any salt of the fissionable 
material, soluble in the liquid, may be used. Slow, fast, 
and spontaneous fission have been observed. 


* About 25 n/g hr. For Th***, < 0.15 n/g hr. 
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d. Natural Background 


Cosmic-ray neutrons are a potential source of 


cavitation nuclei and should be reckoned with under 
some circumstances. Watt [1967] has reviewed some 
of the literature; as an average he takes the total 
flux density to be 65 n/cm2/hr at sea level in the middle 
latitudes. Of the total, most have low energies and 
only 7 percent (4.5 n/em2/hr, are above 10 MeV. Recent 
measurements by Tajima, et al. [1967], indicate a 
much lower flux density, at least in the range 2-10 
MeV, and an altitude dependence like exp (— 0.0063 x), 
x being the atmospheric depth in g/cm?. 

In total flux density, the cosmic-neutron background 
is equivalent to the neutron source in our laboratory 
(1.47 107 n/s into 47) at a distance exceeding 80 
meters, and the former has a smaller fraction of neu- 
trons at the higher energies. 


4. Some Results 


4.1. General 


To fix ideas, we consider water at atmospheric pres- 
sure. The frequency of the sound pressure is in the 
range 25-50 kHz and as a threshold we take the nega- 
tive pressure at which the average waiting time for 
cavitation is (say) 30 sec. The threshold in ordinary tap 
water or in distilled water, both of which contain 
numerous motes, varies markedly with the amount 
of dissolved air, being higher the less air; furthermore 
the threshold is variable from sample to sample and 
from time to time in the same sample. Values range 
from 1 bar or less at saturation (with air) to perhaps 
10 bars at 10 to 30 percent saturation. Barger [1964| 
studied these effects for a wide range of many of the 
variables and Messino et al. [1967], added to the 
water both wettable and unwettable their 
results well illustrate the possibilities. 

We have found the situation to be different, however, 
if the water is cleaned up taking the precautions 
described in section 3.1. In this case the threshold 
rises as more of the motes are removed, no matter 
what the air content, so long as it is below saturation 
[Greenspan and Tschiegg, 1966a, b|. It seems that no 
important purpose is served by protracted cleanup: 
the threshold continues to rise unless (presumably) 
the experiment is terminated by an adventitious 
neutron. 

Table 1 summarizes a few results from the literature. 
The three sets are not strictly comparable although 
Finch duplicated Galloway’s apparatus so far as pos- 
sible. Note the marked effect of the amount of dis- 
solved air. 

Our measurements were taken at 43 kHz. At 46 
kHz, Barger [1964] reports thresholds ranging from 1.2 
bar at 760 torr partial pressure of air to 12.1 bar at 
10 torr. These contrast with our thresholds of about 
160 bar (several minute basis) or 210 bar (several 
second basis) independent of air content from 120- 
650 torr. Presumably the air content has less effect 
the cleaner the water. 


motes; 


TABLE 1. Some cavitation thresholds in water 








Threshold,* bar 





Partial pressure of dissolved air, torr. 





0.2 100 250 700 





Galloway [1954] »:¢ 
Barger [1964]© 
Finch [1964] 


200 * : 1.1 























* Sound pressure amplitudes; static pressure 1 bar. The three 
threshold definitions are roughly equivalent. Frequency ~ 26 
kHz. 

» Interpolated from a curve. 

“Spread about + 12 percent. 

“ According to a private communication, Galloway to Finch 
[1964]! high value gotten sporadically. More often 15 to 20 bars. 

_* Spread not given. 
' Another mode of operation gave 50 to 80 bars. 


4.2. Neutron-Induced Cavitation 


This effect was first found in acetone and in pen- 
tane by Lieberman [1959]. We have observed it 
since in water, isopropanol, and freon 113 (F2CICCCI.F) 
[Greenspan and Tschiegg, 1966a], in ethanol [Green- 
span and Tschiegg, 1966b] and later in another dozen 
or so organic liquids. 

Two results stand out. First, the air content of the 
liquid affects the threshold little if at all. Lieberman 
also found this. Second, there is no appreciable induc- 
tion or decay time for the effect. West and Howlett 
[1967] have since investigated this point in detail; 
in the liquid they used, tetrachloroethylene, the life- 
time of a bubble is about 24s. These results were 
obtained in an elegant experiment using a pulsed 
neutron source. We have already alluded to the results 
of Finch [1964] and of Sette’s group, which contrast 
markedly with, e.g., ours. They find a cumulative 
effect; many minutes or even an hour of exposure to 
neutrons are required to “‘weaken”’ the liquid (water 
in both cases) and the “strength” returns at a similar 
interval after exposure is terminated. Sette and his 
co-workers [Bertolotti et al. 1965, p. 109] feel that 
long lifetimes are characteristic of gassy (i.e., sat- 
urated with air at the static pressure) water; however 
Finch’s water, which showed the effect, was only 
0.5 percent saturated, and our water does not show 


the effect even if gassy. Perhaps dirt is more important 


than gas. 
a. Statistics 


The question which obtrudes itself most forcefully 
is whether or not the cavitations are single-particle 
events, that is, whether each cavitation is traceable 
ts a single neutron or whether each is due to the 
cooperative effect of two or more neutrons. In the 
former case the cavitation rate would (a) follow a 
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Poisson distribution, and (b) satisfy the inverse-square 
law for various distances of the neutron source from 
the cavitation cell. 

Lieberman’s [1959] results in pentane were “rea- 
sonably close to a Poisson distribution.”” We have 
investigated this point in detail; the test liquid was 
isopropanol (table 2). Two independent runs were 
made using the (3, 0, 1)-mode of a lead zirconate 
titanate cylinder of the type shown in figure l(a). 
The temperature in each case was about 31 °C and 
the frequency about 40 kHz. An Am-Be source radiat- 
ing 1.6 X 10° n/s into 47 was used. For the first run. 
the interval distribution was analyzed, and for the 
second, the events per unit time. 

Run |: With the source about 1 meter from the cell 
the apparatus was run at constant negative pressure 
until 245 cavitation events had occurred; the times 
were automatically recorded on a strip chart. On 
account of the dead time (1.31 s) 24 short intervals 
were missed. The total of 268 intervals were arranged 
in 8 groups and analyzed as shown in table 2 using 
the x? test [Evans, 1955, ch. 27]. No statistically sig- 
nificant departure from an exponential ‘distribution 
of interval lengths was observed and hence we may 
regard the cavitation events as occurring individually 
and collectively at random. 


TABLE 2. Test of randomness (interval distribution), in isopropanol 


(~“—=« 


X, 


mm mm 
0 ] 

] 2 
2 4 
4 6 
6 8 
8 

10 13 
13 24 

















| } P 
* Time measured by length / on chart. 
» Observed number of intervals of le ngth / such that /; <1 < Js. 








Expected number of intervals from x, 
— exp (— 0.242 /.)]: 0.242 


Run 2: With the source at about 0.6 meter from the 
cell, the negative pressure was adjusted so that about 
14 events occurred per minute (gross). Events were 
counted during 23 6-min runs. After adjustment for 
dead time (1.44 s/event) the counts were subjected 
_to x” analysis. The results were 


= 268 [exp (— 0.242 /,) 
=mean number of events/mm. 


xX” = 24.6 


F=22 
P= 0.30. 


The average count rate was 21.5/min with,a stand- 
ard deviation from residuals of 6.1/min and an expected 
standard deviation of (21.5)'2=4.6/min. The con- 
clusion is the same as for run 1. 

The inverse-square test is more difficult to do, 
especially in a small room, because of the reflection 
of neutrons from the sides and the furniture. Using 
the same equipment and conditions as before, some 
tests were made on isopropanol and more on freon 
113 [Greenspan and Tschiegg, 1966a|. We indicate 
the results for the freon 113. At a fixed negative pres- 
sure 3 runs were made at each of the distances 0.91, 
1.09, 1.28, 1.52, 1.83, and 2.1] meter. The products, 
counts/unit time by square of distance, were averaged 
for the 18 runs and used to calculate the expected 
count for each run. The x? test gave P=0.18. The 
experiment was repeated at a pressure amplitude 2 dB 
higher with the result P=0.22. It thus appears that 
the cavitation rate at fixed pressure is proportional 
to neutton flux. West and Howlett [1967] have per- 
formed a similar experiment on tetrachloroethylene 
using a Po-Be source. They feel the inverse-square 
law to have been verified within the experimental 
error. 


b. Threshold 

As a preliminary to the definition of a neutron thres- 
hold we consider how the cavitation rate varies with 
the negative pressure, —p. Lieberman [1959] intro- 
duced the concept of “effective sensitive volume,” 
V.: his argument may be paraphrased as follows. In 
the hypothetical case of a liquid populated by identical 
nuclei, upon each of which cavitation occurs when- 
ever it is in a region in which the negative pressure, 
—p, exceeds a fixed value T, called the threshold, the 
“sensitive volume,” V,, is that volume in which, 
everywhere, — p > T. During each period of the pres- 
sure oscillation, there is a “sensitive time,” ts, during 
which —p>T: this time varies with position in Vx, 
being greater the greater the pressure amplitude and 
zero at the boundary of V,. Then he defined 


V. [ ars 
aD 


~ +e 


(12) 


(7 is the period), to which the cavitation rate might 
be expected to be proportional under the assumed 
conditions. Lieberman calculated, by numerical in- 
tegration, V. as a function of —p/T for the case of a 
spherical resonator (sound pressure amplitude varies 
like sin kr/kr).'° 
His curve, which we have extended downward, 

shown in figure 9. Shown also are our results!' in 
freon 113 for the conditions given in table 3. It is seen 
that the cavitation rate increases more rapidly with 
negative pressure than does the effective volume. 
From a similar experimental result Lieberman was 
able to draw various conclusions about the action of 


® Rough calculations show that the shape of the curve is nearly the same for other pres 
sure distributions 

"The results of Lieberman |1959]. in pentane, 
and Howlett [1967]. in tetrachloroethylene, 


are qualitatively the same; those of West 
very different 
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CAVITATION EVENTS/ min 
Ve , ARBITRARY UNITS 











4 6 
—p, dB reT 


FIGURE 9. Cavitation rate negative for freon 113; 


neutron irradiated. 


versus pressure 


See table 3 for conditions. Lieberman's effective sensitive volume is shown in arbitrary 
units 


neutrons in nucleating cavitation. These matters are 
beyond the scope of this paper: we only remark that 
we have made some rough measurements using nearly 
monoergic neutrons which support Lieberman’s ideas. 
Part of the freon 113 data in figure 9 is shown again 
in figure 10 together with data on methanol, isopro- 
panol, and ether. The ordinates have been normalized 
as shown in table 3, which also gives the test condi- 
tions. Although the actual cavitation rate spans the 
range 0.37 per min (in freon 113) to 680 per min (in 
methanol), the points fall reasonably well on a universal 
curve although systematic differences from liquid to 
liquid are unmistakable (compare isopropanol with 
ether). The result is consistent with the following: 
(1) the cavitations are “single-particle” events: 
(2) the range for — p/7, about 6 dB, is too small to 
show differences between the neutron spectra 
of Am-Be and Pu-Be at their upper ends: 
(3) over the same range, geometrical effects (V_) as 
between cylinders and cones are not large. 


TABLE 3. Test conditions for figures 9 and 10 


Liquid Methyl 


alcohol 


Ethyl 


ether 


Isopropyl 


Freon 
113 alcohol 





Source Am-Be Pu-Be Pu-Be Am-Be 





n/sec into 47. 1.6 x 10° 1.5 x 10° 1.6 x 10’ 
Distance, eni 91 91 25 
CMs. <. (°) ‘ (”) (") 
Temperature, °C 31.5 31.5 28.5 
Frequency, kHz... 34.7 39.5 14.8 
Normalization for figure 10° 0.07 0.37 0.47 














* Source adjacent to cell. 

"Ceramic cylinder in (3, 0, 1) mode, see figure 1. 75 em OD by 100 cm high. 
© Conical flask, mode unknown, see figure 3. 125 ml. 

* Cavitations/min corresponding to the ordinate “1.” 


























CAVITATION RATE, NORMALIZED 











-p, dB reT 


FIGURE 10. Cavitation rate versus negative pressure for four liquids; 
neutron irradiated. 


See table 3 for conditions. The ordinate “1” corresponds to 0.07/min for freon 113 (A), 
to 10/min for methyl! alcohol (@), to 0.37/min for isopropyl alcohol (©), and to 0.47/min for 
diethyl ether (xX) 


As can be seen from figures 9 and 10, if a neutron 
source be placed adjacent to the cell, the operator 
can find a value of the electrical readout below which 
cavitations occur rarely and above which they occur 
at a rate of several per minute. The procedure is 
subjective and the reading obtained depends some- 
what on the neutron output of the source and on the 
effective sensitive volume in the liquid, but neverthe- 
less it can in most cases be reproduced to within + 1 
dB at worst and often to within + 0.2 dB. Such readings 
are arbitrarily taken to correspond to thresholds, 
and from a series of them taken at various over- 
pressures, a “neutron threshold” corresponding to 
the particular source is gotten as described in sec- 
tion 2.4b. The main advantage of the procedure is its 
speed and ease. It is desirable that our thresholds 
be comparable to those obtained by others. Unfor- 
tunately very few thresholds for (a, n) sources have 
been measured. For n-pentane, Lieberman [1959] 
got 2.1 bars in the room temperature region; we got 
about 1.3. For acetone his values average to 5.5; 
ours is 4.8. We have a private communication from 
C. West [Harwell] giving 6.5+1.2 bar in tetrachloro- 
ethylene, not far from our 7.8. These are all we can 
find, and although the agreement seems marginal, 
those familiar with cavitation data may find it rather 
good. Both authors cited worked at about 20 kHz, 
and our results are for 28 and 35 kHz in n-pentane, 
for 40 kHz in acetone, and for 49 kHz in CoCl,. We 
have made this comparison because there is some 
evidence that the neutron threshold depends little 
on frequency in the range 20-60 kHz, at least; never- 
theless this. point requires investigation. 
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TABLE 4. Some neutron thresholds, Pu-Be source 


Threshold 


(absolute) 


Normal Fre- 
boiling 
point 


Liquid Temp. quency 








Methanol, CH,O 
Methanol, CH,O 
Ethanol, C.H,O 
}-Propanol, C,;H,O 
1-Butanol, C,HwO 


2-Propanol, C,;H,O 








n-Pentane, 

n-Pentane, C; 
n-Hexane, Cy 
n-Hexane, Cy 


Diethylether, C,H iO 
Diethylether, Cs;HwO 


Freon 113, ¢ F 
Freon 113, C.ClyF 
Freon 113, Cy F 
Freon 113, ¢ F 
Freon 113, CsClyF 


Freon 11, CCIsF 
Freon 11, CC 





Ethyl bromide, C.H;Br 
Ethyl bromide, C.H;Br 
Carbon tetrachloride, CCl, 
Trichloroethylene, C,HCI 
Tetrachloroethylene, C.Cl, 








Methylene chloride, 
Methylene chioride, 
Methylene chloride, 
Methylene chloride, 

















‘In Ba TiO, cylinder, 7.5 cm D by 7.5 em high, (3, 0, 1)-mode. See figure l(a) 

300 ml conical flask, mode unknown. See figure 3 

125 ml conical flask, mode unknown. See figure 3 

“2-propanol is out of place, but is shown for companion 

"In lead zirconate-titanate cylinder, 7.5 cm OD by 7 
l(b) 

' As (e), but mode unknown. 

* 250 ml conical flask, mode unknown. See figure 3 


Some neutron thresholds obtained by the method 
described are given in table 4. These are intended to 
be illustrative only, but some possible trends emerge. 

1. There is some evidence for frequency inde- 
pendence (part 1, methanol: part 2, freon 113: part 4, 
ethyl bromide). 

2. In a homologous series, the neutron threshold in- 
creases as the series is ascended. Part 1 lists 4 nor- 
mal alcohols and part 2, two normal alkanes. 

3. In general the threshold decreases as the tem- 
perature increases. The dependence is stronger for 
temperatures closer to the normal boiling point. 
See parts 2 and 3. To interpret the results for freon 113 
(part 3) this way we have to suppose that these thresh- 
olds are subject to an error of at least 0.2 bar. 

4. One gets the impression that, other things being 
equal, low boiling liquids have a lower threshold than 
high. However, freon 11 and freon 113 contain the 
same elements, and have approximately the same 


5 cm high, (3, 0, 1)-mode. See figure 


TABLE 5. Neutron threshold versus neutron energy: 


T=29 °C, F=44-52 kHz 


— ae 








Neutron Threshold, bar 





Freon 113 Methanol Ethanol 
Boa Se 


2.5 [Did, n)He*] 13.4 


Pu-Be (10.4?) - 6.7 











density, thermal expansion coefficient, surface ten- 
sion, specific heat, and latent heat (per mole);'* yet 
they have about the same threshold at (say) 27 °C 
although the vapor pressure of freon 11 at this tem- 
perature is almost twice that of freon 113. 

5. Note how the effects of atomic composition can 
be separated from those of gross physical proper- 
ties. Examples are the isomers 1-propanol and 2- 
propanol in part 1, and the pair 1-butanol in part 1 
and diethyl ether in part 2, both of which are C,H. 

6. Part 5 is included to show a liquid, methylene 
chloride, in which the neutron thresholds as measured 
vary irregularly with frequency even though the 
original data look good. We suspect that the sensitive 
volumes were so small in relation to the available 
neutron flux, and varied so much from frequency to 
frequency, that the measured “thresholds” are not 
comparable. The neutron thresholds as defined would 
be much more determinate if the source were either 
larger, more nearly monoergic, or both. 

Some miscellaneous results are given in part 4. 


c. Threshold Versus Neutron Energy 

We have made some preliminary measurements of 
threshold as a function of neutron energy. In table 5, it 
is tentatively presumed that what are being seen in the 
Pu-Be source are the neutrons of maximum energy. 
The 2.5 MeV neutrons are from the D—D, and the 
14 MeV, from the D—T, reactions."* The results are 
probably only rough because the neutron generator 
was operated near its minimum output where the flux 
is both unknown and irregular. The implications of 
these results with respect to a principle for neutron 
spectrometry is obvious. 


4.3. Alpha-Recoil Induced Cavitation 


An a-emitter dissolved in a liquid provides cavita- 
tion nuclei in any concentration desired [Greenspan 
and Tschiegg, 1966 a, b]. As Riepe and Hahn [1961] 
who worked with a-emitters in a bubble chamber have 
pointed out, the heavy recoil nucleus, although it 


" The properties listed are those which are considered significant in various theories. 
The viscosity has been omitted, but this is larger for freon 113 than for freon 11, which 
would be expected to further increase the threshold for 113 compared to that for 11. 

'8 We are indebted to Frank Sentfle and his co-workers of the U.S. Geological Survey 
for making available and operating the neutron generator. 
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carried only about 2 percent of the total energy, is the 
responsible agent, because of its very small range rela- 
tive to that of the a-particle. Indeed the range is in 
most cases so small that all of the energy is deposited 
in the immediate neighborhood of the impending cavi- 
tation, so that if all disintegrations were of the same 
energy all cavitation nuclei would correspond to the 
same threshold. In this case the function, cavitation 
rate versus negative pressure, would be expected to be 
nearly the same as Lieberman’s effective sensitive 
volume function already discussed and shown in 
figure 9. 


We have made measurements on methanol and on 
diethyl ether: the a-emitter was ordinary reagent- 
grade thorium nitrate at a concentration of 3 g/l, 
which corresponds to about 300 disintegrations per 
minute (dpm) for the parent and each daughter in 
equilibrium with it. The age of the salt is unknown, 
but rough a-range measurements in air show that a 
considerable amount of the last a-emitter in the chain 
(ThC’) is present. The ThC’ has the highest a-disinte- 
gration energy, 8.95 MeV, in the chain, and if in 
equilibrium would yield 200 a’s/min on account of 
branching in the decay of its parent, ThC, which would 





CAVITATION EVENTS/ min 
Ve , ARBITRARY UNITS 








1 


4 
-p, dBreT 





FIGURE 1]. Cavitation rate versus negative pressure for methanol 
and ethyl ether; a—disintegrations. 
x methanol at 26 °C. in 125 ml conical flask (see fig. 3). f= 37 kHz. 
0 ethyl ether at 28 °C, in similar cell, (see fig. 3) f=45 kHz. 
The solid line is Lieberman's effective sensitive volume curve. 
The dashed lines show it shifted +0.2 dB. 


provide about 100 a’s/min of energy 6.05 MeV. Five 
other a-emitters are also present; their a energies 
range from 4.00 to 6.77 MeV. 


The results are shown in figure 11 superimposed on 
Lieberman’s V. curve in such a way as to give the 
best fit by eye. The dashed curves are drawn to 
show the errors in the cavitation rates produced by a 
+0.2 dB drift in the apparatus. The methanol results 
show good agreement as would be expected if only 
the ThC’ disintegrations were effective and the 
ether results indicate that the lower-energy disinte- 
grations come into play as the negative pressure 
increases. It is hardly likely that one situation prevails 
for methanol and the other for ether, so the results 
must be considered as preliminary to further work 
with single-isotope solutes. The “‘a-thresholds” were 
11.0 bars for methanol and 4.0 bars for ether. 


4.4. Fission 


Messino, et al. [1963], reported that the cavitation 
threshold of water was lowered from 1.6 to 0.55 bar 
by the addition of 5 percent uranyl nitrate, and as- 
cribed the result to fission of the U** fraction (0.7%) 
by cosmic neutrons thermalized by the water. In 
spite of the low cosmic-neutron flux the hypothesis is 
consistent with the authors’ view that the cavitation 
nuclei have a long lifetime, although it is not clear 
that the authors gave due consideration to the com- 
peting sources of spontaneous fission in U*** and of 
the laboratory neutron source (see section 3.3). 


We have observed fast fission in solutions of thorium 
nitrate in isopropanol and in methanol, at thresholds 
not very well defined, but about 10 dB below the 
neutron threshold. The spontaneous fission rate of 
Th? is about 0.15/¢ hr and that from cosmic neu- 
trons < 0.003/g hr [Segre. 1952]. 

A solution of natural uranyl nitrate in clean water 
exposed to neutrons from our Pu-Be source cavitates 
at negative pressures in the range 10 to 20 bars. 
The events appear to be as abrupt as those induced 
by neutrons or a-recoils and cavitation ceases im- 
mediately as the source is encased. No such effect is 
observed at much higher stresses if other heavy 
metal salts are substituted, and without doubt both 
thermal and fast fission are observed in the natural 
uranium, for which the fission cross sections are 
about 4 and 0.5 barns, respectively. 


Uranium in methanol exhibits some peculiar fea- 
tures. Many of the events leave large bubbles which 
are trapped in the pressure nodes, and which persist 
until the sound field is removed, even though the liquid 
has been first degassed. It is difficult to say anything 
about the effect of cosmic neutrons on the fission, 
because we can remove our encased source (23 em 
of paraffin) only to about 40 meters from the cell. At 
this distance the output of the shielded source is still 
comparable to the cosmic flux. Only a few fissions per 
hour are observed. 
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5. Remarks 


We have described several types of apparatus and 
procedures for their use which allow the acoustic 
cavitation properties of liquids to be determined 
with much greater reproducibility than hitherto 
attained. They are susceptible to further improve- 
ment, so that eventually it should be possible to 
make measurements as precise as the statistical 
nature of cavitation properties, and the imperfect 
characterization of the test liquids, permit. 

Contaminants can be removed, with some ease or 
with considerable difficulty, depending on the nature 
of the liquid, so that over a range of (negative) pres- 
sure levels cavitations occur only when nucleated by 
certain ionizing radiations. The cavitations are always 


single-particle events, at least under the conditions of 


our experiments. The results presented are mostly for 
illustration: more systematic investigations, which we 
hope refinements in the apparatus will make possible, 
will yield a body of empirical data which can be 
expected to support some existing theories of the 
process and cast doubt upon others, and perhaps to 
point the way to their modification. Presumably the 
theories of ultrasonic cavitation, of the bubble cham- 
ber, and of the static strength of liquids, are inter- 
related, and some features of these 


processes are 


doubtless shared by hydraulic cavitation as well. 
lonizing particles provide reproducible nuclei for 
studies of ultrasonic cavitation per se, but the pos- 
sibility of using radiation-induced cavitation in nuclear 
instrumentation ‘has intrigued many, beginning with 
Hughes [1960]. Hahn and Peacock [1963] proposed a 


neutron detector discussed at greater length by 
Bertolotti, et al. [1965]: more definitive experiments 
were performed by West [1965]. Our work extends 
these notions to a wide variety of liquids and opens 
up the possible development of a neutron spectrometer. 
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A Kossel pattern generator, designed and built at the National Bureau of Standards, is described 
in detail. The unit is modular and consists of an electron beam column, vacuum system, light micro- 
scope, film cassette, and Kossel camera. The camera component includes microgoniometric capa- 
bilities. The Kossel x-ray technique enables the investigator to obtain lattice spacing data precise 
to two or three parts per million, and orientation of crystals to 0.1° of arc. The Kossel pattern gen- 
erator described permits data for such determinations to be obtained quickly and precisely. 
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instrumentation. 


1. Introduction 


A divergent x-ray beam, emanating from an effective 
point source of x-radiation, produces diffraction cones 
on propagating within a single crystal. These cones 
have a semiapex angle of (90°—6), (6 the Bragg 
angle); their directions are governed by the crystal 
orientation and by the lattice symmetry. When such 
cones are intercepted by a flat photographic film, 
conic sections corresponding to each monochromatic 
component (such as Ka,, Kaz, and KB) of the pri- 
mary beam will be recorded superimposed upon a 
general background blackening. These conic sections 
have been referred to as Kossel lines in honor of that 
pioneer of the divergent beam x-ray technique []].! 
The entire photograph is often called a Kossel pattern, 
or, if the x-ray source happens to be external to the 
diffracting crystal, the photograph may be called a 
pseudo-Kossel pattern [2]. 

Information which can be obtained from the photo- 
graph includes lattice spacings, crystal orientation, 
and data concerning the degree of crystal perfection. 
In particular, lattice spacing data precise to two or 
three parts per million can be obtained. Orientation 
may be determined to accuracies approaching 0.1° of 
arc. Finally, a complete stress-strain analysis may be 
obtained from line profiles and intensities. 

It was known 30 years ago that such information 
could be obtained; nevertheless, the divergeni beam 
method was not widely used until about 1962. In the 
past 5 years, new procedures allowing the investi- 
gator to obtain all of these data in a straightforward, 
rigorous fashion have been developed [3-10]. How- 
ever, a drawback to the use of the divergent beam 
technique was the fact that the development of suit- 
able instrumentation proceeded at a slow pace. This 
paper describes the design of a Kossel pattern gen- 
erator (KPG) which enables the investigator to take 


' Figures in brackets indicate the literature references at the end of this paper. 


full advantage of the divergent beam technique. The 
KPG is a self-contained instrument of modular design 
whose camera component meets all of the preferred 
design criteria described by Vieth and Yakowitz [11]. 


2. Evolution of Kossel Instrumentation 


The first successful experimental demonstration of 
a divergent beam diffraction technique was made by 
Rutherford and Andrade in 1914. These workers used 
an external radium source of y-rays and a cleaved 
rocksalt crystal as the target. From the photograph- 
ically recorded diffraction pattern, they were able to 
determine the monochromatic radiation components 
of the source [12]. 

Three years later Seemann [13] developed a special 
wide angle tube which allowed x rays produced over 
a large area to converge and then to exit through 
a small circular aperture. In 1921, Gerlach [14] de- 
veloped a system in which a fine beam of high energy 
electrons impinged upon a thin foil target. The x rays 
which were produced in this thin foil diverged through 
it as well as the specimen crystal which was placed 
in the opposite side of the foil. This was the fore- 
runner of the modern capillary x-ray tube device 
which is still commercially marketed. 

In 1929, Linnik |15] tried a novel approach using 
conventional x-ray equipment. A fine beam of mono- 
chromatic x rays was allowed to pass through a thin 
crystal. A given crystal volume was made to suc- 
cessively occupy all possible positions within a defined 
solid angle. The rigidly connected crystal-film assembly 
was slowly rotated around an axis which was parallel 
to the x-ray beam and which passed through the 
crystal while remaining perpendicular to the film. 
Simultaneously, the assembly was rocked about an 
axis orthogonal to this which passed through the 
crystal and intersected the rotation axis. 

Kossel [16] designed a system in which the crystal 
of interest was made the anticathode of the x-ray 
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tube. The system was maintained in vacuum and a set 
of apertures was used to produce a narrow beam of 
electrons. Later an electromagnetic lens was used 
to decrease the size of the beam of electrons striking 
the crystal. An aluminum window provided the exit 
port for the diffracted x rays. 

Borrman [17], using conventional x-ray equipment, 
employed a narrow beam of white radiation to cause 
fluerescence of the elements of a crystal. The result- 
ing divergent secondary radiation produced the desired 
patterns. In this method, Laue spots appear in addition 
to the divergent beam pattern.’ 

In addition to being cumbersome and complex, all 
divergent beam devices prior to 1951 suffered the 
handicap of relatively large sources of x-radiation. 
This led to several problems: a large primary source 
causes the Kossel lines to be broadened and impairs 
the inherent precision of the method: in addition, 
exposure times are prohibitively long (24 hr in some 
cases) and photographic contrast is poor. Further- 
more, large single crystals are required to accommo- 
date the large x-ray sources. 

These problems are removed by the use of the 
electron probe microanalyzer. This instrument pro- 
vides a micron sized source of x rays, as well as an 
optical microscope for viewing the area undergoing 
irradiation. Power densities five to ten times greater 
than those obtainable with earlier equipment are 
available. The value of the electron probe micro- 
analyzer as a source for Kossel patterns was recog- 
nized by Castaing as soon as his first instrument was 
completed. He presented Kosse! patterns in his Thesis 
which clearly demonstrated the superior characteris- 
tics of the microprobe in this area [18]. 
Modern 


3. Economics Kossel 


Applied to 
Instrumentation 


All present day Kossel instrumentation utilizes a 
focused electron beam to excite the required micro- 
source of x-radiation. While commercial Kossel 
camera attachments for electron probe micreanalyzers 
represent the most readily available instrumentation, 
investigators have successfully prepared Kossel 
cameras for use in point-projection x-ray microscopes 
and in electron microscopes, possibly since all of these 
instruments provide convenient electron optics. 

However, adapting such instrumentation to electron 
optical columns primarily designed for other purposes 
is usually inefficient. Consider a commercial electron 
probe microanalyzer as an example. Normally, when 
Kossel patterns are being prepared, equipment includ- 
ing x-ray spectrometers, (multichannel analyzer). 
scanning displays, and x-ray read out channels is 
inoperative. Furthermore, almost all Kossel cameras 
built for insertion into existing electron columns con- 
tain serious compromises with preferred design 
characteristics [11]. Hence, it is not only uneconomical 
to tie up a microprobe or electron microscope, but in 
addition, research may be hindered by the camera 
design. 


Another important consideration is that the electron 
optical requirements for the vast majority of Kossel 
line based research are less stringent than those for 
electron probe microanalysis and electron or x-ray 
microscopy. A single electromagnetic lens combined 
with an inexpensive flat grid election gun makes an 
entirely satisfactory electron column.-A relatively inex- 
pensive power supply can be used since gun voltage 
stability requirements are less demanding than for 
electron microanalysis or microscopy. The same is true 
for the single lens power supply needed. To complete 
the system (exclusive of the Kossel camera), a small 
vacuum system, vacuum gauges, suitable light optics, 
and a current meter whose range is 10 to 1000 nano- 
amperes are required. Finally, a modular design can 
be modified at will to incorporate improvements. On 
the other hand, it may prove difficult to modify an 
electron microprobe or microscope for many reasons. 
Based on the foregoing discussion, the decision was 
made to design and build a separate KPG. 


4. Design of the Kossel Pattern Generator 


4.1. Basic Considerations 


Kossel patterns can be obtained either in transmis- 
sion or in back reflection. In the transmission mode, 
the pattern is representative of the entire sample 
thickness traversed. Sample thicknesses are usually 
chosen for wX=2, w being the linear x-ray attenua- 
tion coefficient for the characteristic radiation of 
interest and X being the thickness [19]. In back 
reflection, only regions near the surface are sampled 
so the specimen can be thick. 

All of the information about orientation, lattice 
spacings and strain configuration can be obtained 
from either type of pattern. However, from detailed 
considerations of the two modes, we concluded that 
the transmission pattern is best for precise lattice 
spacing measurements. Back reflection is probably 
superior for most orientation studies [11]. Since 
much of our work is concerned with the study of 
small changes in lattice spacings as a function of 
applied stress or other perturbation, we decided to 
build the KPG suitable for work in the transmission 
mode. This choice only affects the camera module 
design. The electron beam column to be described 
can be used with a back reflection camera as well. 
The unit consists of (1) Kossel camera, (2) Light 
microscope, (3) Film cassette, (4) Electron beam 
column, (5) A suitable vacuum system. Figure 1 shows 
the assembled instrument. 


4.2. The Kossel Camera 


The most important design problem concerns the 
camera module. The remainder of the KPG is designed 
around the camera. The design concepts for a success- 
ful Kossel camera have been detailed previously [11]. 
They may be summarized as follows: 

(1) It should be possible to locate inclusions, twins. 
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FIGURE 1. General view of the KPG fully assembled. 


slip bands, grain boundaries, etc., and place them under 
the beam. This requires capabilities for movement of 
the specimen in orthogonal directions in its own plane 
as well as in a direction perpendicular to the specimen 
plane, and for an optical microscope to view the area 
of interest, 

(2) To eliminate the distortion effects of geometric 
projection on the flat film, carry out orientation studies, 
and for almost all lattice spacing studies, the capability 
of orienting the specimen by means of a micro- 
goniometer controlled from outside the vacuum 
system is necessary, 

(3) The pattern center position must be found for al- 
most all lattice spacing studies and as an origin for 
orientation studies. Therefore, one must be able to 
mark the position of the pattern center precisely, 

(4) Small values of the specimen-to-film distance 
may cause contrast losses on the film while large values 
may increase exposure times to prohibitive values. A 
compromise is a specimen-to-film distance of 5 to 
11 cm, 

(5) The angle between the important cubic poles 
(100) and (111) is about 55° so that an x-ray divergence 
angle of at least 56° is desired on the film, and 

(6) Convenience and control of operation are 
greatly increased by designing the camera so that the 
specimen is in the vacuum system but the film is in 
air. 

The entire system was designed so that it could be 


dismantled and repaired with a minimum of incon- 
venience while still satisfying all of the design concepts 
listed. The camera is a separate module removable 
from the electron beam column by means of three bolts. 

The ,stage allows translation of the specimen 
amounting to 4/2 in in orthogonal directions. Each 
motion is spring loaded to minimize backlash and 
each is connected by means of vacuum shaft seals to a 
micrometer shaft. Thus, the coordinates of a given 
point on the specimen can be found and recorded so 
that the point can be found again later. The vertical 
motion of the stage which makes it possible to focus 
the optical microscope is controlled by a worm and 
wheel arrangement. The wheel has 230 teeth and ad- 
vances a four lead thread 5/s in per revolution. This 
gives a rapid vertical motion but retains enough sensi- 
tivity to focus easily in the light microscope. 

A two-axis goniometer is mounted on this mechanical 
stage. Each tilt motion can be controlled from outside 
the vacuum system. An indicating dial is used to index 
the angle of each axis. The tilt is controlled by a worm 
and a worm-gear: each tilt mechanism is spring loaded 
to minimize the effects of backlash. We were forced to 
use a small worm-gear in view of space limitations. 
This limited the number of teeth and. hence the angular 
sensitivity of the goniometer to a minimum readable 
motion of about 0.2°. The goniometer may be tilted 
+12° about each axis. At this angle the edge of the 
stage strikes the light optical objective lens. 

The essential part of the goniometer design as well 
as of the vertical specimen motion is the linkage con! 
necting the external control device to the actual mov- 
ing part. Each linkage can move through an angle up 
to 15°, traverse a distance of 0.70 in, and rotate freely. 
This is accomplished by attaching small universal 
joints on the ends of standard sliding tubular linkage. 
In order to prevent galling of the sliding mechanism 
in vacuum, we lubricate each part with a suspension 
of MoS, in diffusion pump oil. These linkage com- 
ponents are shown in place in figure 2. 

Finding and marking the pattern center is carried 
out by a method devised by Ogilvie [20]. A small 
aperture at the specimen plane is brought to the 
coordinate point of the beam. Light from the hot 
filament is allowed to strike the film; the high voltage 
is switched off during this operation. A small spot is 
impressed at the pattern center by this means. 

The specimen holder will accommodate a ¥/4 in 
diam sample. Holders are of brass and the specimen 
is sandwiched between two brass discs. A source foil 
for pseudo-Kossel studies may be placed on the 
top disc where it is not in thermal contact with the 
actual specimen. The entire assembly fits into an 
insulating ring: a lead wire connects to a vacuum 
electrical feedthrough. In this fashion, electron cur- 
rent from the specimen can be read on an appro- 
priate electrometer. 

For transmission studies, the specimen-to-film dis- 
tance has little effect on the photographic contrast, 
which is controlled by the operating voltage, specimen 
thickness and crystal perfection [19]. We chose a 
specimen-to-film distance of 12 cm as a compromise be- 
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FIGURE 2. Kossel camera showing linkage for X-Y-Z and tilt motions. 
Specimen in place. 


tween exposure time, angular divergence to be re- 
corded on the film and Kossel line distribution on 
the pattern. With this distance, a maximum divergence 
angle of 70° is recorded on the film. This value is 
almost always satisfactory for obtaining Kossel 
pattern geometry suitable for precisely determining 
lattice spacings [9]. Furthermore, sufficient informa- 
tion is provided for indexing and orientation de- 
termination of most crystals with the 70° divergence 
angle. 


Patterns are recorded on standard 5X7 in x-ray 
film or glass plates. The film cassette is an ordinary 
plate-holder which fits into an aluminum frame bolted 
onto the camera bottom. The film is in air so that 
many exposures may be taken without disturbing the 


crystal or the vacuum. This is accomplished by 
placing a polyester window 5 mils thick directly 
below the specimen (fig. 3). The window diameter is 
17/3 in. The window sits in a removable cap and is 
sealed by O-rings inside and outside. Window failure 
is rare provided normal care is exercised. 


4.3. Optical Microscope 

We found refracting optics viewing the specimen 
at a 35 deg angle to the horizontal by means of appro- 
priately placed mirrors to be unsatisfactory. Image 
distortion was severe and the focusing characteristics 
ltoo critical for convenient use. Therefore, the system 
chosen consists of a modified commercially available 
reflecting objective (50, numerical aperture 0.56, 
working distance 5.5 mm), vertical illuminator and 
ocular. The objective is placed between the pole piece 
of the electromagnetic lens and the specimen. Con- 
trols are provided so that the axis of the light optical 
objective can be made to coincide with the axis of 
the electron beam. Figure 4 shows the light optical 
system. 


FIGURE 3. Bottom view of Kossel camera showing polyester window 
with specimen visible beyond. 


FIGURE 4. View of the column showing the light optical system: 
objective, transfer mirror, centering devices, light source, vertical 
illuminator and ocular. 
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4 
4.4. The Electron Beam Cclumn 


The column design is based on the general discus- 
sion of single lens systems given by Cosslett, Nixon, 
and Pearson [21]. A magnetic lens was chosen because 
its spherical aberration constant, Cy, is smaller than 
that of a corresponding electrostatic lens. The smallest 
beam that can be obtained with an appropriately 
placed aperture, which accepts the full beam current, 
requires C, to be minimized. In order to accommodate 
the light optical system, the specimen must be placed 
1.8 in (4.6 cm) below the nose of the lens. This results 
in a focal length of 5 em for the lens. This long focal 
length gives a value of Cy, equal to 67 cm with the 
asymmetric pole piece geometry of our lens. Never- 
theless, calculations showed that a final beam diameter 
of about 1 yu could be obtained at the specimen plane. 
This required a 25 yu diam source to be limited by an 
aperture to 4.6X 10-4 radians. Placing an aperture 
100 yw in diameter some 5¥2 in (14 cm) below the 
electron gun filament achieved the desired result. 

The electron gun chosen is a commercially available 
flat grid cap type usually supplied with small electron 
microscopes. It is mounted on the axis of the lens. 
The gun bias is operated nearly at the cutoff point in 
order to reduce beam angle and to operate where the 
source brightness is greatest. The filament is a hairpin 
loop of 8 mil diam tungsten wire. Centering of the 
filament within the grid cap as well as the filament 
tip height adjustment with respect to the cap can only 
be carried out with the gun removed from the column. 


FIGURE 5. 


This has not proved to be a great liability. We have 
provided controls for centering the gun assembly and 
the aperture with respect to the axis of the electro- 
magnetic lens. 

The power supply is capable of providing a voltage 
between 5 and 50 kV. The filament heater, bias re- 
sistance and lens current supply complete the require- 
ments for obtaining the electron beam. A suitable 
vacuum system consisting of forepump, 2-in diffusion 
pump, appropriate valves, vauges and piping and a 
cold-bafle completes the basic column. A circulating 
refrigerant cold trap is used so as to eliminate re- 
charging problems. 

The electron beam is monitored by an indicating 
nanoammeter (range 1—10,000 nA) equipped with a 
charge integrator. This device is connected to the 
specimen: the specimen is isolated from ground so 
that the electrometer can read the electron charge 
“absorbed” in the specimen. This current, generated 
by the electron flow, is called the specimen current. 
Several advantages result from being able to monitor 
the specimen current. First, we use it to know when the 
gun is saturated. After electron beam alinement within 
the column, variation of the specimen current is a 
measure of system instability. It is also helpful to 
use the integrated charge to obtain the proper photo- 
graphic exposure. This is done by using the equations 
relating exposure and incident electrons given’ by 
Yakowitz and Vieth [22]. The advantage of this pro- 


cedure is that slight variations in the specimen cur- 
rent can be tolerated by using the total charge. 


Transmission Kossel pattern of Be-4.4 w/o Cu alloy using Co-K radiation. 


he position of the KB circle, only faintly visible on the print, is shown by arrows. The Ka circle is also marked. 


317 


275-286 O-67—5 





5. Performance 


The KPG has been used to study SrTiO;, nearly 
perfect silicon, Fe-3 wt. pct. Si alloy, indium, beryllium 
and beryllium alloys. The case of Be-4.4 wt. pet. Cu 
will serve as a brief example. This alloy is interesting 
since conventional x-ray methods for determining its 
lattice constants are difficult to apply. Furthermore 
the alloy has a hexagonal close-packed structure: 
traditionally, the Kossel method has been applied to 
cubic structures. 

Witt et al., determined a general method for the 
solution of lattice spacings in noncubic structures by 
the Kossel technique [10]. This method requires one. 
two, or even three different crystal orientations for 
success. Thus, the necessity for a goniometer and the 
ability to take many photographs rapidly are underlined. 

The method was applied to the Be-4.4 wt. pet. Cu 
alloy: Co-K radiation was used in the pseudo-Kossel 
experiment (for details, see [10]). All patterns were 
prepared in the KPG: orientation was carried out so 
that the (00.1) pole would be exactly perpendicular to 
the film. The resulting pattern is shown in figure 5. By 
measuring the ratio of the Ka, to KB circle diameters, 
the lattice constants of the alloy were determined as 
c=3.5908 A and a=2.2838 A giving c/a=1.572. The 
entire experiment to obtain figure 5 required less than 
2 hr. Without the goniometric capabilities, the experi- 
ment could not have been properly performed. With- 
out the capability for preparing many photographs 
without disturbing the crystal or the vacuum, we esti- 
mate the time to obtain figure 5 at about 10 working hr. 


Special thanks are due to L. Schneider and C. Taylor 
of the NBS Instrument Shops for their cooperation in 
the manufacture of the KPG. 
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An improved formula is derived for accurately computing the near-zone magnetic field of a small 
circular transmitting loop antenna. Such a field can serve as a reference standard for calibrating 
field-strength meters employing small receiving loop antennas in the frequency range 30 Hz to 30 MHz. 

This formula includes correction terms for frequency (due to the finite time of propagation), as 
well as corrections for the finite radii of both the transmitting and receiving -loops. Other formulas 
appearing in the literature often fail to include such corrections which can result in errors of up to 
20 percent and more in computing standard-field values. 

The NBS formula is derived by expanding the integrand of the retarded vector potential into an 
infinite series of spherical Hankel functions of increasing order. The resulting series expression is in 
error by less than 0.2 percent, is rapidly converging and simple to use without recourse to a table of 


functions or a computer. 
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1. Introduction 


The purpose of this paper is to present a derivation of 
an improved high-frequency formula intended for use 
in accurately computing the near-zone magnetic com- 
ponent of the electromagnetic field of a small circular 
transmitting loop. Such a field can serve as a reference 
standard for calibrating field-strength meters employ- 
ing small receiving loop antennas in the frequency 
range 30 Hz to 30 MHz. The new formula includes 
correction terms for frequency (due to the finite time 
of propagation), as well as corrections for the finite 
radii of both the transmitting and receiving loops. 
Formulas used by other workers in the field often do 
not include such corrections. This new formula is a 
rapidly converging infinite series. It is simple in form 
and easy to use by workers in the field without recourse 
to a table of functions, or to a computer. 

An earlier formula, which prompted this derivation, 
has been used during the past 15 or 20 years at the 
National Bureau of Standards [1].' It has served as the 
basis for the NBS calibration service for commercial 
and military CW field-strength meters over the above 
frequency range. This formula was a single-term ex- 
pression and was first obtained from a detailed numeri- 
cal analysis of the existing low-frequency formulas, of 
which there were many in the literature [2]. Most of 
these were in the form of an infinite series. Many were 
slowly converging and required the use of at least two 
or three terms of the series to yield the accuracy 


*Radio Standards Engineering Division, National Bureau of Standards, Boulder, Colo. 
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' Figures in brackets indicate the literature references at the end of this paper. 


achieved in the single term of the earlier NBS formula, 
over the usual range of the parameters involved. 

From the form and symmetry of this single-term 
expression, it became apparent that it might be the 
leading term of a rapidly converging infinite-series 
expression for the magnetic field of the loop. Efforts 
were then directed toward finding an analytical solu- 
tion in series form having the earlier formula as its 
leading term. A number of methods were tried without 
success. The one that finally achieved the desired re- 
sult is based on expanding the integrand of the re- 
tarded-vector-potential into an infinite series of spheri- 
cal Hankel functions of increasing order. While this 
approach must have occurred to other workers in the 
past, the author has been unable to find any record of a 
high-frequency derivation in this particular series 
form in the literature. 

All time-varying quantities used in this paper are 
expressed in terms of their rms values. Rationalized 
mks units are used throughout. 


2. Calibration and Use of a Small Loop 


The physical arrangement used at NBS for cali- 
brating small receiving loop antennas is shown in 
figure 1. The standard magnetic field is produced by 
a small, single-turn, circular transmitting loop of 
radius, r; meters, in which a current, / amperes, is 
flowing. The small receiving loop being calibrated 
has a radius, r2 meters, and is positioned coaxially 
with respect to the transmitting. loop. The two loops 
are spaced an axial distance, d meters, apart during 
the calibration. It is assumed that the distance of 
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FIGURE 1. Coaxial ~arrangement of transmitting loop antenna 
(radius r; meters), and receiving loop antenna (radius r2 meters), 
spacing d meters, in the NBS calibration setup. 

An RF 


current, / amperes, flows in the transmitting loop 


separation, d, is small compared to the wavelength, 
\, and that near-field conditions prevail. 

The response or indication of the field-strength 
meter is directly proportional to the voltage induced 
in the receiving loop. The induced voltage is, in turn, 
proportional to the time-rate of change of the total 
flux cutting the loop. This is given by Faraday’s law 
of electromagnetic induction, written here for the 
sinusoidal time-varying case (with the time factor 
e/ understood), 


$E- T.=—jo | | B- Be (1) 


The voltage induced in the loop is given by the left- 
hand side of eq (1), which is the line integral of the 
tangential component of the incident electric field, E, 
around the periphery, s., of the receiving loop: B 
is the magnetic-flux density at any point on the plane 
surface, S., bounded by the circular receiving loop: 
w=27f, andj=V~—!. 

When making a field-strength measurement, one 
does not know the spatial distribution of B over the 
area, So, and it cannot be resolved by the loop. The 
response can therefore only be expressed in terms of 
the average field over the area of the loop, and eq (1) 
becomes 

V=—jwB,, * So, (2) 
where }’ is the induced voltage and B,, is the normal 
component of B averaged over the area, So, of the 
receiving loop. This assumes that the relative orienta- 
tion of the receiving loop with ‘respect to the electric- 
field vector, FE, is such that the electric-dipole response 
is negligible [3]. 


3. Retarded Vector Potential 
Transmitting Loop 


of a Circular 


The basic far-field properties of a small transmitting 
loop antenna are usually derived for conditions under 
which second-order corrections to the field due to the 
finite radii, r; and rz, of the transmitting and receiving 
loops, can be neglected. In near-field use, however, 
r; and rz are usually not negligible compared to the 


distance of separation, d, between the loops. There- 
fore, corresponding corrections must be provided in 
the field magnitudes involved by including both of 
these radii in the derivation. 

The expression for the axial magnetic component 
of the near field of the electrically-small (2ar,; <A) 
transmitting loop antenna will now be derived. It will 
be assumed for the present, that the loop current is 
constant in amplitude and phase around the loop, and 
that the radius of the wire with which the loop is wound 
is negligible compared to the loop radius, r;, itself. 
The effect of a non-uniform current distribution on the 
accuracy of the resulting field formulas at frequencies 
above 5 or 10 MHz will be discussed later in section 5. 
The effect of the finite time of propagation will be 
included in the derivation. uh 

The magnetic-flux density, B, at any point, P, on 
the surface, S2, of the receiving loop, is given by the 
curl of the retarded vector potential, 4, at the same 
point [4], i.e., 


B=V XA. (3) 


The geometry involved is shown in figure 2. The aver- 
age value of the normal component of B will be given by 


st] 
_ 


pf=permeability of free-space (u=47- 
meter) 


Bay = <A + dS», (4) 


where [4] 


and 


henry/ 


/=transmitting-loop current, rms amperes 
B= wavelength constant (B= 277/A) 

A= free-space wavelength, meters 

f =frequency, hertz 

j=V-1. 


It can easily be shown that the vector potential will 
have only a component, Ag, which is azimuthal at the 
point, P. Let P be chosen arbitrarily to lie in the xz 
plane where 6=0. Equidistant line elements of length 
ds; at +q and —@ are paired. These elements can 
each be resolved into a component, ds; sin @, parallel 
to the xz plane and a component, ds; cos @, normal to 
the xz plane. In taking the resultant, the former com- 
ponents will cancel, while the latter (azimuthal with 
respect to the point P) will add. Equation (5) can 
therefore be rewritten [5], letting ds; = ridd, 


T p—JjBR 
Ay= 4 | —— cos o dd. (6) 
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By Stokes’s Theorem [6], the integral of the normal 
component of the curl of 4 over the surface S2 bounded 
by the receiving loop can be transformed into the line 
integral of the tangential component of 4 around its 
periphery, i.e., 


[| VXxA-dS.= $ Ag(P)dse, (7) 


where Ag(P) is the value of the vector potential at the 
point, P, now assumed to lie on the periphery of the 
receiving loop of radius, p=rs. Substituting eq (7) 
in eq (4) gives 


| 
Bay = $ Ag(P)dsz2. (8) 


Since the receiving and transmitting loops are coaxial, 
it can be seen from the circular symmetry involved, 
that Ag(P) will be constant around the periphery of 
the receiving loop, so that eq (8) becomes simply 


r2 
Pay — S. étq ). (9) 


Substituting eq (6) in eq (9), letting B, 
So= 77r3, gives 


w= MH, and 


Ir, T e JBR 


T7T2 Jo R- 


H.a= cos @ dd, 


where now 


R=[@+r3+rR—2rnre cos b|'?, (11) 
as can be determined from figure 2. H,, is the average 
normal component of magnetic field strength over the 
area of the circular receiving loop. 


4. Evaluation of the Complex Intearal 

The problem, now, is to evaluate the Hy in 
eq (10). Maxwell [7] was the first to accomp]i$i: t:is tor 
the static case (8=0), using elliptic inte me which 
can be expressed in several series forms [8]. Numerous 
other workers have since evaluated this integral in 
various ways for the static case. Literally, over one- 
hundred different series expressions appear in the 
literature [9] for the mutual inductance between 
coaxial circular filaments, to which the present prob- 
lem is directly related. 

To the above will now be added still another series 
solution, this time also including the effect of retarded 
time. Equation (10) will be evaluated using the rela- 
tionship [10] 

eniBe wexem - 
—e 
where h2\(8R) is the zero-order spherical Hankel 
function of the second kind. Such a function is often 
used to represent an outward-traveling (+) wave, 


FIGURE 2. Geometry used for determining the vector potential and 
near-zone magnetic field of a small circular transmitting loop. 
The loc 


»p lies in a plane normal to the xz plane with its center at the origin of coordinates. 


whereas a similar function of the first kind, h‘))(BR), 
would have been used to represent an inward-travel- 
ing (—) wave. The function, AG (BR). will now be 
expanded into an infinite series of higher order spheri- 
cal Hankel functions. These will be substituted back 
into eq (10) and integrated term by term. 

The expansion can be accomplished using the fol- 
lowing relationship based on Erdélyi [11}. 

h2)(B =h2(Vu+v) = =S — | ee 

m! 


|’ he(Va). 
ye Var 


(13) 
where m is an integer. This series converges when 
|u|>|v|. Let 

Vutv= B(R§—2rrz cos h)"? 
where Rj= @ + + 75. 
Let u= B°R3, 
and v=— B?(2rr2 cos ). 


Substituting eq (16) and eq (17) in eq (13) gives 


AS) (BR) = s = E ay h2)(BRo) cos” o. 


m=0 


Substituting eq (18) back into eq (10) gives 


Aut - 


Be 


m=0 


H.=- 


=e aay" h®2)(BRo) cos™*! odd. 
(19) 
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The infinite series in eq (19) will converge [11] when 
Ri > 2rnr2 if 6=0 or if 6=7. Furthermore it can be 
shown [12] using, for example, the Weierstrass M 
Test, that eq (19) will converge uniformly in the inter- 
val 0 <@<7. Equation (19) can then be integrated 
term by term, yielding 


_; Bln < =| nal (2) c _ 
J ar > m! B y. hi (BRo) ¥ cos em 
(ZU) 


- m=0 


The integral in eq (20) can be evaluated from the follew- 
ing relations [13]: 


(" cos’ bdddb= 0. (if p is odd) 
stad (21) 


od [ero « «>. p—! ee 
[ cos’ ddd = 7. oe 7. (if p is even) 
0 &* A eS aa oe 4 (22) 


When these results are substituted into eg (20), with an 
appropriate change in the index of summation, eq (20) 
becomes 





re Se (mt ijiz*4-6 . . . 2m+21 Re 


URE "Gen SEES, Ss Sh SR | Ee Mena (Be). 


The spherical Hankel functions in eg (23) can now be 
changed into their exponential polynomial forms as 
given by eq (All) to eq (A16) in the Appendix. When 
this is done, eq (23) can be written 


IS, i 15 /rire\? Set: Oe pean 
A = IaRs | | + jBRo| Tv 8. (Re) [ 7 jBRo- 15 BRE aa 1b pRi| 
ars\4 . i ee | See reer 

=) | I + jBRo— 945 BRS ae — BPRY q a BRS TJ 945 BPRS 


945 : 


»—IBKo | amperes/meter. 


Under the conditions that BR» < 1.0, and rr/R?2 
< 1/16, terms in BRo higher than the first power can 
be ignored without introducing an error in eq (24) of 
more than 0.2 percent. Therefore, if the term (1+jBRo) 
is factored out of eq (24), the magnitude of eq (24) can 
be written as 


IS; f 5 le - 316 le 4 . 
a nll: Sieead she cee “3 -% a. 22 2)\1/2 
\Hy|= cat 8 ( z) 64 (je) ae | (1 +6°R>)*". 


This is the desired expression for the normal com- 
ponent of the magnetic field averaged over the area, | 
S., of the receiving loop. 

The effect of the neglected terms in BRo on the 
accuracy of eq (25) decreases rapidly as BR» is pro- 
gressively decreased below 1.0, other factors remaining 
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the same. At the lower frequencies where correc- 
tions due to BRo can be neglected entirely (BRo < 1), 
the mathematical accuracy of eq (25) is limited only 
by the number of terms used in the infinite series. 

For the conditions previously given, the first cor- 
rection term in the infinite series of eq (25) will con- 
tribute less than 1 percent, so that for this case eq (25) 
can be further simplified to 


IS 
[Hay = onR 


(1+ B2RB)'2 (26) 


, amperes/meter. 


It should be noted that eq (26) is essentially identical 
to the expression for the axial component of the 
magnetic field of a classical, infinitesimal magnetic 
dipole. The difference is that the distance factor, Ro, 
in eq (26) includes correction terms for the finite radii, 
r; and r2,.of the transmitting and receiving loops, 
respectively, as given by eq (15). 

As previously stated, the response of a field-strength 
meter employing a small-loop antenna is directly 
proportional to the average normal component of the 
magnetic field strength, H,.. incident on the loop. 
However, for many years it has been customary to 
express the calibration of such an instrument in terms 
of the equivalent electric field strength, E, that would 
be associated with H for the case of a uniform plane 
wave (where E =1207H). This relationship, when 
substituted in eq (26) with S, = 7r?, gives 


a Omri 


eav| = 


+ B*R?)' 


volts/meter. (27) 
R3 

Needless to say, this manner of expressing the re- 
sponse will be valid on/y when the field-strength meter 
is used to measure a uniform plane wave. The meas- 
urement will be of questionable value for a// cases in 
which the local |£/H| ratio departs from ihe free-space 
value [14]. 


Standard-field formulas appearing in the literature 


[15] to [17] often fail to correct for the finite radii of 


the transmitting or receiving loop antennas, or for 
frequency. It is of interest to evaluate the errors 
resulting from the failure to apply these corrections, 
all of which are included in the formulas derived in 
this paper. 

Failure to correct for the radius, ro, of the receiving 
loop will result in an error of from 1.0 percent (for 
d=1.25 m, r2=0.1 m), to 20 percent (when d= 0.5 m, 
r2=0.2 m). Failure to simultaneously correct for the 
radius, r;, of the transmitting loop will approximately 
double these errors. 


The value of field strength given by eq (25) to 
eq (27) is essentially independent of frequency at 
the lower frequencies (where BRo <1). Failure to 
apply the frequency correction indicated, at the higher 
frequencies, will introduce an error of, e.g., 1.0 per- 
cent at 5 MHz, which increases to approximately 27 
percent at 30 MHz (when d=1.25m, m=r2=0.1m). 
If the transmitting loop circumference exceeds 
approximately A/16, the effect of a nonuniform current 


FIGURE 3. Geometry used to analyze the current distribution on a 
circular transmitting loop antenna. 
The loop is tree? 
circuited at th 
is measured fr 


ed as a circular balanced transmission line, fed at points 1, 2, and short- 
‘iving end (point 3). Distance along the loop periphery, /= 6r,, meters, 
vint 3. 


distribution will also have to be taken into account, 
as discussed in the next section. 


5. Loop-Antenna Current Distribution 


In the preceding derivation, the loop current distri- 
bution was assumed, for simplicity, to be constant in 
amplitude and phase around the loop. Actually, a 
standing wave of current exists on the loop, and 
the assumption of a constant current will not be valid 
unless the loop is small electrically (i.e., unless 
271, Xd). 

The complete solution to this problem for an 
arbitrary current distribution [18] is quite complex. 
However, the effect of a nonuniform current distribu- 
tion at frequencies from roughly 5 to 30 MHz can be 
determined from the following approximate analysis. 
The single-turn loop can be considered as a circular, 
balanced transmission line fed at points 1, 2, as 
indicated in figure 3, and short-circuited at its receiving 
end (point 3). 

If the radius of the wire with which the loop is wound 
is small compared to the loop radius, m, it can be 
assumed that the current distribution around the loop 
is, to a first approximation, of the hyperbolic-cosine 
form [19], 


l=I cosh yl (28) 


where 


/)=reference current at point 3 (fig. 3) 

Y = propagation constant=a+ JB 

a = attenuation constant, nepers/meter 

B=wavelength constant=27/A, radians per meter 

A= free-space wavelength, meters 

/=peripheral distance, in meters, 
point 3. 


measured from 
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The loop will be considered sufficiently loss free 
that @ can be assumed negligible compared to B. 
This is a reasonable assumption for the type of loop 
used, so that eq (28) can be written 

I = Io cosh (jBl) = Io cos Bl. (29) 
As can be seen from eq (29) the current will have a 
constant phase around the loop for this case. 

For the size loop of interest here, the magnitude of 
the magnetic field strength given by eq (24) to eq (26) 
will be, to a first approximation, proportional to the 
average value of the current around the transmitting 
loop. The average current is given by 


hh 
= { 


ie 


cos Bldl =v — . (30) 
P 1 


where /,;=half the loop circumference, meters and 
/=0Or, meters from figure 3. 


TABLE 1. Percent difference between I, and lo 


2rr; = 
- {1- 7, ) * 100 


Current 
difference. 
percent 


Loop 
circumference, 
wavelengths 





Table 1, which is based on eq (30), shows the 
approximate percentage difference between /,, and 
I) for several sizes of low-loss loops (a <), with 
their circumferences expressed in fractions of a 
wavelength, 27r,/A. As can be seen, the correction 
due to the nonuniform current distribution will be 
less than 1.0 percent at the highest operating fre- 
quency, provided the loop circumference does not 
exceed A/16. This correction decreases rapidly 
with decreasing frequency, becoming less than 0.2 
percent at frequencies below 15 MHz for a loop 
having a radius r;=0.1 m. 

These corrections, based on an assumed cosinus- 
vidal current distribution around the loop, are in quite 
close agreement with similar corrections based on a 
recent theoretical analysis by T. T. Wu [20]. The dif- 
ference in the correction values obtained by the two 
methods in less than 0.02 percent for the conditions 
involved here. It is interesting to note that for this case 
the loop current distribution obtained by Wu reduces 
.essentially to the cosine form given by eq (29) above. 
This is shown in figure 4 for a loop having a radius 
r,;=0.1 m operating at a frequency f=30 MHz. 
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FiGURE 4. Relative current distribution around a small circular 
transmitting loop antenna having a radius 1, 0.1 m, fora fre- 
quency f= 30 MHz. 

The current distribution. based on an analysis by T. T. Wu [20]. is 


the text, where 


compared with an 


assumed cosine distribution determined from eg (29) in @r, from figure 3. 


\ wire radius a =0.001 m was used in the Wu analysis 


6. An Alternate Method for Evaluating the 
Integral 


The infinite series of eq (25) can be easily verified 
for the static case (8 = 0) using a much simpler method. 
This is done by rewriting eq (10) in the form 


Ir; 


: F(d) cos ddd. 
7Rore 0 


Ha, 


2rire cos db a2 


where = (32) 
R? 


F(d) = (1 o 


The binomial expansion of eq (32) is then substituted 
back into eq (31) and the result integrated term by 
term, yielding the infinite series of eq (25) directly, 
provided R32 > 2rir2 as before. The first few terms of 
the binomial expansion of eq (32) are as follows: 


nr ao 
—- ) cos” 


F(d)=]1 +{ =e) cos +5 | 


YO 


ran 
P 3: 
cos? @¢+— 
8 


(33) 


63 /rire\° 5 
eee cos? b+. 


Be 

Substituting eq (33) back into eq (31) and perform- 
ing the required integration, making use of eq (21) 
and eq (22), yields for the static case 
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IS, 15 rye 2 Sie rye 4 
av =5.P3 1+>( 2) +22 ( 2) + 
‘ 27R3 8 \ R 64 \ RF 
which is identical to the infinite series in eq (25). 


7. Summary and Conclusions 


An improved high-frequency formula has been de- 
rived for accurately computing the near-zone magnetic 
field of a small (2mr, <A), circular, transmitting 
loop antenna. Such a field is used as a reference 
standard at NBS for the calibration of field-strength 
meters [1], [14] employing small receiving loop an- 
tennas in the frequency range 30 Hz to 30 MHz. 

The new formula is in the form of an infinite series 
and includes correction terms for frequency (due to 
the finite time of propagation), as well as corrections 
for the finite radii of both the transmitting and receiv- 


ing loops. The derivation is based on a current of 


constant amplitude and phase around the transmitting 
loop. However, the effect of a nonuniform current dis- 
tribution at frequencies above 5 or 10 MHz is also 
analyzed. Formulas used by other workers in the field 
[15] to [17] often fail to include these refinements 
which can result in errors of up to 20 percent.and more 
in computing standard-field values. 

The accuracy of the new formula is basically limited 
only by the number of terms used in the series. At the 
highest frequency normally employed, 30° MHz, the 
error in one of the approximate forms given does not 
exceed 0.2 percent for the conditions stated and de- 
creases rapidly with decreasing frequency. The series 
expression is rapidly converging and simple to use by 
workers in the field without recourse to a table of 
functions or a Computer. 
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mathematics group) for their assistance and helpful 
suggestions relating to this problem. The author also 
acknowledges with thanks the assistance of Ezra B. 
Larsen who made a detailed numerical analysis of 
T. T. Wu’s work which was used in the preparation 
of figure 4. 


8. Appendix 


For the convenience of the reader, explicit expres- 
sions for the various Bessel functions used herein, but 
not commonly found in compilations, are listed below. 

The spherical Hankel functions of the second kind, 
used in this paper, are defined in terms of Bessel func- 
tions of the first and second kinds, J,(x) and N,(x) 
respectively, by the relation [21] 


| : 


Kx) = “a (Jn + s/2(x) —JNn+ 1/2(x)), (Al) 


Bessel functions of odd half order have particularly 
simple forms in terms of trignometric functions [22]. 
A few of these relationships are as follows: 
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The recurrence relations for calculating these 
tions for higher orders are [23] 
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If these two relations are substituted into the definition 
of the spherical Hankel functions, given by eq (Al) of 
this appendix, the recurrence relation for spherical 
Hankel functions is obtained as [24] 
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Using the above relations several spherical Hankel 
functions of higher order can be determined. 


» e a 
h(x) =— 


x 


(+7) 


h(x) = 


AY) (x) =- 


h(x) = 





jr 45 
h(x) = (++ 2 -jS- (A15) 


0s 108) 
x ; os 


x: 


esr 15105 .420 945 .945 
acca Nap ad Sala ee ie ee 
x \ : x x™ x* x ee 


(A16) 


It should be noted that the expression for each of 
the above functions is exact, each being a polynomial 
of degree (n+1). These Hankel functions can also 
be derived somewhat more directly from the following 
- . -. = for 
relation based on Erdélyi [25]: 
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Thermal expansion of rutile from 100 to 700 °K, R. K. Kirby, 
J. Res. NBS 71 A, No. 5, 363-369 (Sept.—Oct. 1967). 

The thermal expansion of a single crystal of rutile (TiO) was deter-- 
mined in directions parallel to its crystallographic axes. A unique 
macroscopic technique was used wherein a cube-shaped specimen 
was cut from the anisotropic crystal, each of its six faces polished flat 
and nearly parallel to the face opposite it, and the expansion meas- 
ured along each of its three body axes with an Abbe-Pulfrich inter- 
ferometer. The relationship between the expansion along the body 
axes of the specimen and that in the crystallographic directions 
was determined from a Laue x-ray pattern. The unusual behavior 
of the thermal expansion that was observed can be explained by 
assuming an acoustical and optical contribution where the corre- 
sponding Griineisen constants are 2.80 and 0.75 respectively. 


Key Words: Anisotropic crystal, Griineisen’s constant, interferom- 
eter, rutile, thermal expansion. 


Criterion for the stability of numerical integration methods 
for the solution of systems of differential equations, A. I. 
Karim. J. Res. NBS 71B, Nos. 2 & 3, 91-103 (Apr—Sept. 1967). 
The problem of studying the growth of the ‘error is most important 
for the numerical solution of differential equations. In this paper 
the Wilfs criterion is generalized to be applied for systems of 
differential equations. A general theorem is investigated and regions 
of stability have to be determined. The use of an electronic com- 
puter is more essential for such regions to be characterized. These 
regions of stability have the property that, the error introduced at 
any stage tends to decay. The regions of stability for particular 
numerical methods are explicitly determined. 


Key Words: Numerical integration, systems of differential equations, 
stability, propagation of error. 


Discrete complex functions with prescribed boundary values 
and residues, E Peterson, J. Res. 71B, Nos. 2 & 3, 105-110 
(Apr.Sept. 1967). 

R. Isaacs, J. Ferrand, R. J. Duffin, and several others have developed 
a function theory for “discrete analytic functions” defined on 
“discrete regions” in the “discrete complex plane.” In this paper 
we bring to light some combinatorial-topological properties of 
“simple discrete regions,” and we study some basic properties of 
discrete analytic functions that are defined on simple discrete 
regions. These combinatorial-topological properties and_ basic 
properties are then used to establish an existerice and uniqueness 
theorem for discrete complex functions with prescribed “boundary 
values” and “residues” on an arbitrary simple discrete region. 
Key Words: Analytical functions, Dirichlet 
problem, discrete analytical functions. 


complex analysis, 


Numerical solution of second-order linear difference equa- 
tions, F. W. J. Olver, J. Res. NBS 71B, Nos. 2 & 3, 111-129 (Apr 
Sept. 1967). ; 

A new algorithm is given for computing the solution of any second- 
order linear difference equation which is applicable when simple 
of recurrence steps, (ii) applying to inhomogeneous difference equa- 
tions, (iii) enabling more powerful error analyses to be constructed. 
The method is illustrated by numerical computations, including 
error analyses, of Anger-Weber, Struve, and Bessel functions, and 
the solution of a differential equation in Chebyshev series. 


Key Words: Chebyshev series, difference equations, error analysis, 
Miller algorithm, recurrence methods, special functions. 


Indefinite integrals involving Bessel functions, B. A. Peavy, 
J’ Res. NBS 71B, Nos. 2 & 3, 131-141 (Apr.—Sept. 1967). 


Expressions are derived for the indefinite integrals, 
[y (r)\Colar)dr 


| rf (r)C2(ar)dr 


[ rf (r)Co(ar)Co(Br)dr 


[y (r)Co(ar)Zo(Ar)dr 


where Co(ar) are zero order Bessel functions, Zo(Ar) are zero order 
modified Bessel functions and f(r) is a polynomial in r. In general, 
the expressions given for the integrals are given in terms of pre- 
scribed functions of the Bessel functions, and the coefficients of 
these functions are determined from a finite series, the terms of 
which are found from recurrence relationships that involve only the 
polynomial f(r). Coefficients of the terms of the finite series are 
given in tabular form for up to an eleventh degree polynomial. 


Key Words: Bessel functions, indefinite integrals. 


Bibliography and index on vacuum and low pressure meas- 
urement. January 1960to December 1965, W. G. Brombacher, 
Mono. 35 Supplement 1 (May 31, 1967), 65 cents. 

The bibliography, a supplement to NBS Monograph 35, of the same 
title, covers the years 1960 to 1965, inclusive, and includes 109 
references for prior years not listed in Monograph 35. It contains 
1787 references, 38 of which are to books. The references, besides 
those directly concerned with pressure measurement, include those 
on vacuum technology which may bear on the technique of vacuum 
measurement, such as on adsorption and desorption, diffusion of 
gases into solids, gas conductance and hardware such as pumps, 
seals and traps. Author and subject indices are provided. 


Key Words: Bibliography, low pressure measurement, techniques 
of vacuum measurement, vacuum index, vacuum measurement. 


Automatic typographic-quality typesetting techniques: a 
state-of-the-art review, M. E. Stevens and J. L. Little, Mono. 99 
(April 7, 1967), 70 cents. 

This report describes the current state-of-the-art in automation of 
graphic arts composition, starting from either of two sources: 
(1) keyboard entry of manuscript material, or (2) mechanized input 
in the form of available perforated tapes or magnetic tapes. The 
gamut is covered from one extreme in which a skilled keyboard 
operator performs all of the compositor functions required to operate 
a typesetting machine, to the other extreme in which the input merely 
provides text whether or not including designation of desired font 
changes, followed by a high degree of automation through all opera- 
tions leading to type set for printing. Intermediate automation aids 
for the compositor functions, including characteristics of special- 
purpose digital computers and functions performed by typography 
programs for general-purpose digital computers, are reviewed. 
Characteristics of automatically operated typesetting mechanisms, 
including hot metal casting machines and photocomposers, slow, 
medium and high speed, are outlined. Applications of new tech- 
niques for typographic-quality automated composition that are of 
interest in scientific and technical information centers, libraries and 
other documentation operations include sequential card camera 
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listings, computer-generated KWIC indexes, photocomposition of 
technical journals, automatic composition of books containing both 
computer-produced tabular data and natural language texts, and the 
incorporation of mechanized processes throughout the publication 
cycle from the author’s original manuscript preparation to the final 
printing. A bibliography of 363 references is included. 


Key Words: Automated type composition, computerized typesetting, 
graphic arts, photocomposition, printing and publication. 


NBS standard frequency and time services, Radio stations 
WWV, WWVH, WWVB, WWVL, Misc. Publ. 236, 1967 Edition 
(1967), 15 cents. 

Detailed descriptions are given of eight technical services provided 
by the National Bureau «{ Standards radio stations WWV, WWVH, 
WWVB, and WWVL. These services are: 1. Standard radio fre- 
quencies; 2. Standard audio frequencies; 3. Standard musical pitch; 
4. Standard time intervals: 5. Time signals; 6. UT2 corrections: 
7. Radio propagation forecasts: and 8. Geophysical alerts. In order 
to provide users with the best possible services, occasional changes 
in the broadcasting schedules are required. This publication shows 
the schedules in effect on June 1, 1967. Annual revisions will be 
made. Advance notices of changes occurring between revisions 
will be sent to regular users of these services upon request. Current 
data relating to standard frequencies and time signals are also 
available monthly in the Time and Frequency Services Bulletin. 


Key Words: Broadcast of standard frequencies, high frequency, 
low frequency, standard frequencies, 
quency. 


time signals, very low fre- 


Units of weight and measure. International (metric) and U.S. 
Customary, L. J. Chisholm, Misc. Publ. 286 (May 1967), $1.50 
(Supersedes Miscellaneous Publication 233). 

The primary purpose of this publication is to make available the 
most often needed weights and measures conversion tables —con- 
versions between the U.S. Customary System and International 
(Metric) System. A secondary purpose is to present a brief historical 
outline of the International (Metric) System —following it from its 
country of origin, France, through its progress in the United States. 


Key Words: Conversion tables, International System (SI), Metric 
System, U.S. Customary System, weights and measures, weights 
and measures abbreviations, weights and measures systems, weights 
and measures units. 


Discussion of errors in gain measurements of standard 
electromagnetic horns, R. W. Beatty, Tech. Note 35] (March 
1967), 30 cents. 

In setting up a calibration service for measuring the gain of standard 
electromagnetic horns, one needs a reference horn in which one has 


developed a high degree of confidence. Although it is possible to 
calculate the gain of horns of certain design, confidence can be 
increased by carefully measuring the gain. This note examines a 
method for measuring the gain of two identical horns, listing the 
assumptions made in making such a measurement. The theory of 
2-port waveguide junctions is applied to the analysis of the measure- 
ment technique. The method is shown to be essentially an attenua- 
tion measurement which has additional sources of error. Although 
these errors are not analyzed and evaluated in this note, the problem 
is perhaps more clearly stated than it was previously. The mismatch 
error in comparing two horns as receiving antennas is analyzed. 
Data is given on the aperture efficiency of standard horns which 
indicates that improvements in the design of such horns are feasible. 
It is concluded that, at present, an uncertainty limit of the order of 
tenths of decibels seems realistic, but hundredths of decibels seems 
unattainable until further refinements are made both in the standard 
horns themselves and in the measurement techniques. 


Key Words: Antennas, calibration, comparison, effective aperture, 
electromagnetic horns, gain, measurement, microwave, mismatch 
errors, Rayleigh distance, standard gain horns. 


The viscosity and thermal conductivity coefficients of dilute 
neon, krypton, and xenon, H. J. M. Hanley and G. E. Childs, 
Tech. Note 352 (March 23, 1967), 25 cents. 


The coefficients of viscosity and thermal conductivity for dilute 
neon, krypton, and xenon were examined by a method already proved 
successful for dilute argon, oxygen, and nitrogen. This method 
selects a suitable potential function, and its parameters, which is 
then used to correlate theory with experimental data, given the 
kinetic theory expressions for the transport coefficients. The method 
has recently been expanded and generalized and the results of this 
general study are applied in this note. The potential functions ex- 
amined were members of the m-6, Kihara, Exp: 6, and Morse families. 
It was found that the Kihara was most suitable for neon, and the 
m-6, with m=17 and m=24, was most suitable for krypton and 
and xenon, respectively. Viscosity and thermal conductivities were 
calculated from these functions and tables are given between 100 


and 100 °K. 


Key Words: Dilute gases, neon, krypton, xenon, transport coefh- 
cients, correlations, m-6, Kihara, Exp: 6, Morse, potential functions. 


Connector for saturated standard cells, J. J. Barth, Tech. Note 
353 (April 21, 1967), 10 cents. 

This paper describes a connector used for making electrical connec- 
tions to saturated standard cells in oil baths. The connector has at 
least three advantages over other methods in common use. These 
advantages are namely: (1) it does not generate thermoelectromotive 
forces; (2) it allows cell racks to be placed close together in the bath; 
and (3) it is economical to fabricate since it can be made from a 
common inside caliper. 


Key Words: Connector, electromotive force, oil baths, standard 
cells, thermoelectric. 

Thermal radiation property measurement techniques, S. T. 
Dunn, J. C. Geist, D. G. Moore, H. E. Clark, and J. C. Richmond, 
Tech. Note 415 (April 27, 1967), 50 cents. 

This is an annual summary report of work completed on NASA Con- 
tract R-09-022-032. The work comprised (1) completion of the 
development and calibration of a rotating cylinder procedure for 
measuring normal spectral emittance of non-conducting materials 
at temperatures in the range of 1200 to 1600 °K, (2) analysis and 
calibration of an ellipsoidal mirror reflectometer, and (3) a study of 
the relation between surface roughness and geometric distribution 
of flux reflected from a surface. 


Key Words: Averaging spheres, diffuse reflectance, emissivity, 
emittance, flux averaging devices, infrared reflectance, spectral 
reflectance, spectral emittance, specular reflectance. 


Adhesive bonding of various materials to hard tooth tissues. 
IV. Forces developing in direct-filling materials during 
hardening, R. L. Bowen J. Am. Dental Assoc. 74, No. 3, 440-445 
(Feb. 1967). 

Tensile forces will develop when dental filling materials shrink as 
they harden within a cavity, if there is bonding to the cavity walls. 
A test method was devised by which these forces could be measured. 
This was done by observing the loading required to maintain constant 
distance between two opposing walls (having known area) of an 
artificial cavity. If several broad assumptions are made, it can be 
estimated that a tensile stress of at least 49 kg/cm? (700 psi) at the 
walls of a cavity may develop during the hardening of a direct filling 
resin. The corresponding stress for a silicate cement would be about 
35 kg/cm? (500 psi). Other direct filling materials were also investi- 
gated. The magnitude of the stress depend on the exact test method, 
the material, and other factors. 

The strength of adhesive bonding between a direct filling material 
and the cavity walls of a tooth must exceed the tensile stresses that 
develop during the filling material, if the bonding is to remain intact. 


Key Words: Adhesion, bonding, dental, filling materials, hardening, 
shrinkage, strain, stress, tensile. 


Effects of packing pressures on the properties of spherical 
alloy amalgams, G. T. Eden and R. M. Waterstrat, J. Am. Dental 
Assoc. 74, No. 5, 1024-1029 (Apr. 1967). 

The properties of dental amalgams prepared from spherical-alloy 
particles are compared with the properties of amalgams prepared 
from irregularly shaped conventional alloy particles. The spherical- 
alloy amalgams were found to be as strong or stronger than the 
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conventional amalgams and in addition the spherical-alloy amalgams 
possess a superior property of “adapting” their shape to that of 
the cavity walls particularly when they are condensed under low 
packing pressure. The spherical alloy amalgams also are capable of 
retaining a greater fraction of their strength at low packing pressures 
where the conventional amalgams lose their strengths rapidly. 

The use of spherical-alloy amalgams offers the dentist superior 
“adaptability” and better control over such variables as _rate-of- 
hardening, setting expansion strength, carvability, mercury content, 
etc., through a more effective control of alloy particle size and shape. 


Key Words: Amalgams, spherical-alloy, dentistry, packing 
spherical-alloy amalgam, mercury-silver alloys, 
particle size, adaptability, dental amalgam. 


pressure, 
silver-tin alloys, 


Bonding porcelain teeth to acrylic resin denture bases, 
G. C. Paffenbarger, W. T. Sweeney, and R. L. Bowen, J. Am. Dental 
Assoc. 74, 1018-1023 (Apr. 1967). 

The chemical bonding of porcelain teeth to both cold- and heat- 
curing denture bases of acrylic resins was accomplished by the use 
of a silane coupling agent — gamma-methacryloxy-propyltrimethoxy- 
silane. In tension the porcelain-resin interface did not rupture— 
invariably the teeth did. Untreated teeth showed no measureable 
bonding to the resin. Because of a lack of adhesion between porce- 
lain teeth and acrylic resins the design of most current porcelain 
teeth provides for mechanical locking using metal pins or undercut 
wells and canals (diatoric form). These designs weaken the teeth 
unduly and can be eliminated by the use of a: solid tooth which has 
been treated with an appropriate silane. The bonding strengthens 
the denture since it makes the effective cross-sectional area greater 
and prevents seepage between the tooth and the plastic and so aids 
oral hygiene and esthetics. 


Key Words: Teeth, resin, silane, denture, plastic, porcelain, bonding, 
adhesion, coupling agent. 


Flow properties of aqueous suspensions containing kaolins 
of varying degrees of crystallinity, W. Ormsby and J. H. 
Marcus, J. Am. Ceramic Soc. 50, No. 4, 190-195 (Apr. 1967). 

A study was made of the influence which varying degrees of crystal 
perfection of several Georgia kaolins had on the flow properties of 
aqueous systems containing these clays. The fractionated clay 
samples were characterized with respect to clay mineral composi- 
tion, crystallinity, surface area, and cation exchange capacity. A 
recording Ferranti cone-plate viscometer equipped with truncated 
cones was used to examine the clay-water mixtures under a variety 
of experimental conditions. Emphasis was placed on high shear-iate 
measurements of systems of intermediate water content. 

A comparison of results for bulk samples and fine particle-size 
fractions of different clays showed a general qualitative correlation 
between viscosity (or consistency) and degree of crystallinity. In 
these samples, which were usually pseudo-plastic, the viscosity 
tended to increase with decrease in crystallinity. Coarse particle- 
size fractions of the various clays exhibited essentially Newtonian 
behavior while samples of intermediate particle size and relatively 
good crystallinity exhibited dilatancy. The extent 6f dilatancy was 
usually a function of the degree of crystallinity. 

The viscosity of the fine particle-size fractions of a single, well- 
crystallized clay, which exhibited an apparent increase in crystal- 
linity with decrease in particle size, increased as the particle size 
decreased. This correlation between viscosity and particle size was 
also noted for the fine fractions. of a poorly crystallized clay in 
which the crystallinity was essentially _inde »pendent of the particle 
size. 

In general, results indicated that dilatant and associated properties 
were controlled by the degree of crystal perfection whereas the vis- 
cous and plastic properties, at a given water content, were generally 
dependent upon the particle size. 


Key Words: Flow properties, kaolins, crystallinity, cation exchange 
capacity, surface area, cone-plate viscometer. 

NBS hunts for better ways to control corrosion, S. Lichten- 
stein, Mater. Eng. 65, No. 5, 21-23 (May 1967). 

NBS contributions to a mosaic of evidence on the nature of cor- 
rosion processes are described, including thin-film and passivity 
studies under ultra-clean conditions, underground studies, polar- 


ization and cathodic protection, galvanic effects, and stress and 
atmospheric corrosion in marine environments. 


Key Words: Corrosion, passivity, ultra-high: vacua, ellipsometry, 
low-energy electron diffraction, field-electron emission microscopy, 
NBS corrosion section. 


The relationship between laboratory tests on solubility of 
zine oxide-eugenol type cement and their behavior in the 
mouth, G. M. Brauer, NV. Y. J. Dentistry 37, No. 4, 146 (Apr. 1967). 
The solubility and disintegration in water and dilute acids of zinc 
oxide-eugenol (ZOE) and especially modified cements of improved 
physical properties containing o-ethoxybenzoic acid (EBA) is 
considerably lower than that of other dental cements. Contrary 
to these laboratory findings EBA fillings fail rapidly under oral 
conditions. Thus, a low solubility value may be useful in comparing 
the relative solubility of products of similar composition, but is 
not necessarily an indication of the success of such restorations in 
the mouth. Development of a reliable test that accurately predicts 
length of clinical service will require more basic knowledge of the 
material and of the many parameters encountered in the oral 
cavity. 


Key Words: Zine oxide-eugenol cement, solubility, o-ethoxybenzoic 
acid, filling material. 


A new service by NBS. An aid in making the Bureau’s tech- 
nical resources available to private standards organizations, 
manufacturers, and government agencies, G. S. Gordon, 
Proc. 15th Annual Meeting Standards Engineering Society, Atlantic 
City, N.J., Sept. 14, 1966, pp. 73-75 (1966). 

The Office of Engineering Standards Liaison & Analysis was estab- 
lished as a focal point for assistance to standards-making bodies. 
Its functions will be to catalyze, strengthen existing institutions, 
identify new areas of technical activity in the standards field, 
assist the internal technical engineering stds. activities at NBS 
in reaching maximum effectiveness, and build bridges of communi- 
cation between NBS and both government and non-government 
groups in national and international standardization. The LaQue 
Committee in 1965 recognized the need for developing a national 
stds. effort by tying together the efforts of several hundred stds.- 
making organizations and developing the concept of the USA 
Standards Institute, and recommended expansion of the Bureau’s 


cooperation with industry and USASI in making technical experts 


available to serve on committees and as delegates to international 
ISO, IEC, and COPANT meetings. 


Key Words: OESLA established at NBS as focal point for assistance 
to standards-making bodies. 


Absolute calibration of vibration generators with time- 
sharing computer as integral part of system, B. F. Payne, 
Shock Vibration Bull, 36, Pt. 6 183-194 (Feb. 1967). 

Improved shakers with a simplified ceramic moving element and a 
reference accelerometer now permit reciprocity calibration over a 
wider frequency range, from 10 Hz to 5000 Hz. This paper describes 
present NBS shaker calibration procedure, using a time sharing 
computer, in which magnitude and phase lag measurements are 
made. Calibrations to 10 kHz are planned. 

A teletypewriter gives access to a central commercial computer 
which can accept commands and data from punched paper tape or 
from the keyboard. The data may be recorded from digital measur- 
ing instruments by a special coupling system. Data reduction is 
rapid and permits errors to be corrected quickly. The combined use 
of the computer and the digital data recording system cuts the 
calibration time to about one-fourth the time required previously. 
Monitoring the quality of motion of the shaker through the use of 
the computer makes possible improved accuracy of the calibration 
process. 

As an example, a calibration of one of the improved shakers is 
presented. 


Key Words: Time shared computers, shakers, calibration, acceler- 
ometers, vibration standards. 


Cathodic protection of copper in a severely corrosive soil, 
W. J. Schwerdtfeger, IEEE Trans. Ind. Gen. Appl. 1GA-3, No. 1, 
66-69 (Jan.—Feb. 1967). 
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This paper describes a laboratory procedure used in developing 
electrical requirements pertinent to the cathodic protection of 
copper for a given situation. The work was prompted by severe 
corrosion of copper tubing used underground at a housing develop- 
ment. The soil environment was reported as having a low resistivity 
(about 200 Q-cm) and a high concentration of sulfate ion. 
Specimens of the same type of copper tubing, one arranged to 
corrode freely (control) and the other under cathodic protection 
(cathode), were exposed in the laboratory to a sample of the soil for 
a period of 130 days. 

Polarization data obtained on the control were used to evaluate the 
cathodic protection simultaneously applied to the cathode. After 
130 days, the specimens were removed from the soil and cleaned. 
The cathode was well protected, the metal loss being less than 3 per- 
cent of that on the control which was almost perforated in a few 
places 


Key Words: Copper, soil corrosion, corrosion current, corrosion 
rate, potential, polarization, cathodic protection. 


Evaluation of a stark voltmeter, Y. Beers and T. W. Russell, 
IEEE Trans. Instr. Meas. IM-15, No. 4, 380-388 (Dec. 1966). 
An absorption Stark voltmeter using a Fabry-Perot absorption cell 
with a rotational transition of CH;CN at 36 GHz has been built 
and evaluated. The rotational constant and dipole moment of this 
molecule have been redetermined and found to have the values 
9198.56+0.10 GHz and 3.9185+0.0020 Debye, respectively. The 
short term precision of the device as used as a standard cell was 
better than one part in 10°, but there were unexplained drifts of 
2 parts in 10* per day. The accuracy of the device as a means of 
measuring voltage in terms of atomic and molecular constants is 
5 parts in 10*. The principal errors are due to the measurement of 
the plate spacing (which could be reduced by using an optical 
interferometer to measure spacing), pulling the interferometer 
resonance by the external circuit, and the effects of field non- 
uniformities due to the coupling holes. 


Key Words: CHsCN (Acetonitrile), Stark effect, Fabry-Perot inter- 
ferometer, dipole moment. 


Improved transfer standard for vibration pickups, E. Jones, 
D. Lee, and S. Edeiman, J. Acoust. Soc. Am. 41, No. 2, 354-357 
(Feb. 1967). 

This paper describes a vibration transfer standard designed to 
provide comparison calibrations of pickups with minimum degrada- 
‘tion. Features of the design are the use of a ceramic housing of a 
stiff, light, lossy material, provision of a means for evaluating the 
quality of the motion by the use of three integral accelerometers 
oriented parallel to the nominal direction of motion, and design of 
the geometry to minimize differences in motion between the pickup 
being calibrated and the standard. The useful frequency range is 
10 to 10,000 Hz. 


Key Words: Comparison calibration, vibration pickups, acceler- 
ometers, transfer standard, ceramic, axial motion, transverse motion. 


Introduction, J. M. Richardson, (Proc. 1966 Conf. Precision Elec- 
tromagnetic Measurements, Boulder, Colorado, June 21-23, 1966), 
IEEE Trans. Instr. Meas. IM-15, No. 4, 138 (Dec. 1966). 

It is appropriate, at each biennial Conference on Precision Elec- 
tromagnetic Measurements, to examine the importance of this 
series in the scheme of things. Recently we have come to recog- 
nize that the field of physical measurements has functional ele- 
ments and complicated interactions characteristic of a vast system. 
We may call this complex the system of physical measurement. 
There is a comparatively small group of people who work at the 
limits of attainable accuracy in order to keep the system valid and 
current. There is also a huge number of people who use the system 
as a means to the ends of science, engineering, manufacturing, 
trade, and everyday life. This conference exists so that the keepers 
of the system can communicate their work to each other and also 
to users of the system. No other conference addresses itself ex- 
clusively to the dynamically changing techniques of precision 
electromagnetic measurement and to how they can be put to use. 


Stress-corrosion cracking of cold-reduced austenitic stain- 
less steels, H. L. Logan and M. J. McBee, Mater. Res. Std. 7, 
No. 4, 137-145 (Apr. 1967). 


The effects of specimen orientation with respect to the direction 
of rolling or prior rolling on the susceptibility to stress-corrosion 
cracking in a boiling MgCl. solution were determined for both 
annealed. and cold-reduced types 301, 304, 310, and 321 stainless 
steels. Threshold stresses for the most susceptible orientations 
were determined in a boiling 0.5 N NaCl, 0.1 N NaNO: solution to 
be 60 to 110 percent of the yield strengthsof the materials. None 
of the steels cracked in a marine atmosphere during 15 months’ 
exposure stressed at 90 percent of their yield strengths. The effects 
on the threshold stresses of varying amounts of cold work differed 
from steel to steel and are reported. The effects of internal stress 
patterns and welding were also investigated for some of these 
steels and are described. 


Key Words: Stress corrosion, corrosion, cracking, steels, stainless 
steels, cold working, atmospheric corrosion tests, chemical attack. 


An intercomparison of hydrogen and cesium frequency 
standards, R. Vessot, H. Peters, J. Vanier, R. Beehler, D. Halford, 
R. Harrach, D. Allen, D. Glaze, C. Snider, J. Barnes, L. Cutler, and 
L. Bodily, JEEE Trans. Instr. Meas. IM-15, No. 4, 165-176 (Dec. 
1966). 

Intercomparisons of average frequency and of frequency stability 
were made among one Hewlett-Packard 5060A cesium beam, two 
Varian Associates H—10 atomic hydrogen masers, and the National 
Bureau of Standards NBS III cesium beam designated as the 
United States Frequency Standard. Each of the standards displayed 
a white noise frequency fluctuation behavior with a transition into an 
approximate flicker of frequency fluctuation behavior for longer 
time intervals. The rms fractional frequency fluctuation between 
adjacent samples, a(t, N=2), was 6 X 10~''z~'/? down to a flicker 
level of about 3X 10~-" for the hp 5060A cesium beam (107 <7 S 10's). 
1X 10-''7-"? down to a flicker level of less than 1 X 10-'" for NBS 
III cesium beam (10? <7 < 104s), and 5 X 10-'°7~"/? down to a flicker 
level of about 1 X 10~'4 for the H—-10 hydrogen masers (1 <7 =< 104s). 
The accuracy capabilities of NBS III and H—-10 #4 are now 1.1 X10-" 
and 0.47 X 10-'*, resnectively (lo estimate). 

A discrepancy of only 1.1 parts in 10'* was observed between the 
average frequencies of the hp 5060A cesium beam and the NBS III 
cesium beam, with the former being higher in frequency. 

In terms of the frequency of the Cs'** hyperfine transition (F = 4, 
mp =0)< (F =3, mp=0), defined as 9 192 631 770.0000 Hertz, the 
measured frequency of the H' hyperfine transition (F=1, mp =0) 
<—F=0, mp=0 (was vy=1 420 405 751.7864+0.0017 Hertz. This 
is believed to be the most accurate and precise measurement of any 
physical quantity. 


Interval scales, ratio scales, and additive scales for the sizes 
of differences perceived between members of a geodesic 
series of colors, D. B: Judd, J. Opt. Soc. Am. 57, No. 3, 380-386 
(Mar. 1967). 

From larger-smaller judgments of color differences, compared 
visually two at a time, the perceived sizes may be evaluated on an 
interval scale. Given numbers B so evaluated, and such that B is 
linearly connected to some power p of the physical measure D (such 
as distance on any chromaticity diagram) of the differences, the addi- 
tive constant Ky, such that the numbers B+ Ky, are expressed on a 
ratio scale, may be found from the judgments of the ratio of sizes of 
pairs oi differences. To evaluate p, it is sufficient to observe the 
three differences, 12, 23, and 13 among three colors, 1, 2, 3, forming 
a geodesic series, and chosen so that By. is not much different from 
B.3. The scale formed by the numbers (B+ K,,)'/” is additive if 
the D-scale is additive. If the largest of the color differences judged 
exceeds the smallest by a factor not greater than 3, a close approxi- 
mation to the (B+ K»,)'/-scale may be found without evaluat- 
ing Ky, by ratio judgments. This approximation is based on the 
empirical discovery that scales based on the additivity condition: 
(Biz + Koa)!” + (Bos + Koa)” =(Bi3+ Koa)”, though it implies that 
choice of p between 1 and 1/3. It is sufficient therefore, to derive the 
additive scale by setting p=1, and computing Kya as Bi3— Bi 


— Bas. 


Key Words: Additive scales, color geodesics, color scales, geodesics, 
interval scales, ratio scales, scales. 


Reflectance spectrophotometry, D. B. Judd, Enclycopedia of 
Industrial Chemical Analysis, 3, 376-392 (John Wiley & Sons, Inc. 
New York, N.Y., 1966). 
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Five commercially available spectrophotometers adapted to the 
measurement of light-diffusing specimens are described. The use of 
these instruments for qualitative and quantitative analysis of such 
specimens by means of the Kubelka-Munk model is presented to- 
gether with the methods for formulating colorants required to 
produce both nonmetameric and metameric color matches. 


The national measurement systems of various countries, 
H. W. Lance, JEEE Trans. Instr. Meas. IM-15, No. 4, 397-403 
(Dec. 1966). 

This paper is a summary of an informal panel discussion which 
formed one session of the 1966 Conference on Precision Electromag- 
netic Measurements. In contrast to the other sessions, this one was 
concerned primarily with broad management problems rather than 
with detailed technical matters. The members of the panel, from nine 
different countries, had been asked to consider the measurement 
system of their country, their national standards laboratory (or its 
equivalent), and the scientific and technical program of the labora- 
tory. Although time limitations prevented a thorough discussion, 
many worthwhile facts and opinions emerged. 


The variability of color measurement, |. Nineroff, Color Eng. 5, 
No. 2, 24-29 (Mar.—Apr. 1966). 

Because all the components in the equations for computing chro- 
maticity coordinates in spectrophotometric colorimetry are meas- 
ured, these components are subject to measurement uncertainty. 
The problem of determining the chromaticity uncertainty can be 
treated most simply by the application of the theory of propagation 
of errors in a computed result. The ellipse is chosen as the most 
convenient closed plane curve by which to represent the area of 
uncertainty in the chromaticity diagram within which the true 
chromaticity point can be expected to lie a specified fraction of the 
time. The theory has been applied to determining chromaticity 
uncertainty resulting separately from uncertainties in spectral 
irradiance of a source, spectral emittance of objects and spectral 
response of normal human observers. In actual practice the effects 
of these uncertainties should be computed coincidently. On the 
basis of experience it is suggested that a complete standard observer 
system, consisting of means, variances, and covariances of the 
spectral tristimulus values, in xyz and uvw coordinates, be recom- 
mended by Commission Internationale de l’Eclairage (CIE) for use 
in color measurement. 


Key Words: Colorimetry, propagation of errors, uncertainty in 
chromaticity coordinates, uncertainty in spectral tristimulus values, 
uncertainty in spectrophotometric data, uniform chromaticity 
spacing, variability. 

Thermal converters for audio-frequency voltage measure- 
ments of high accuracy, F..L. Hermach and E. S. Williams, 
IEEE Trans. Instr. Meas. IM-15, No. 4, 360-368 (Dec. 1966). 
The ac-de differences of a reference group of thermoelements have 
been evaluated at audio frequencies to a few ppm at currents 
from 5 to 20 mA. A technique for comparing the ac-de differences 
of two thermoelements with an uncertainty of about 2 ppm has been 
developed. Two 5 mA thermoelements are used with a plug-in set 
of resistors of computable reactances to form thermal voltage con- 
verters for voltage measurements. With this same technique adja- 
cent ranges of these converters can be compared to step up from 
0.5 to 500 V to better than 10 ppm. \ 


Key Words: Thermoelements, thermal converters, voltage measure- 
ment, audio-frequency, standards, transfer standards, comparator, 
voltage comparator. 


WWV moves to Colorado. Part I, Y. Beers, QST 51, No. 1, 
11-14 (Jan. 1967); Part Il, QST 51, No. 2, 30-35 (Feb. 1967). 

This is a semi-popular account of the move of WWV from Greenbelt 
to Fort Collins, written at the request of the Editor of the leading 
amateur radio magazine “QST”. It describes the role of NBS in 
providing standards of time and frequency, the history of the 
broadcasts, the usage of the broadcasts, and the new station. 


An optical power limiter, B. L. Danielson, Appl. Opt. 6, No. 1, 
158-159 (Jan. 1967). 

A power limiting detector is described which is suitable for use in 
monitoring laser or other high intensity radiation in the 2900 to 
3400 A region. This device employs the saturable single crystal 


phosphor CaF.(Ce, Mn) in connection with a photomultiplier tube 
which is sensitive to the sensitized Stokes-shifted Mn? * fluorescence. 
The maximum fluorescent intensity of this crystal can be controlled 
by appropriate choice of impurity concentrations. Limiting the light 
flux in this way minimizes fatigue effects and possible damage due 
to inadvertent overloading of the photomultiplier tube. 


Key Words: Optical power limiter, light flux, photomultiplier tube, 
fluorescence, saturation. 


Compensation and band tailing effects in high power room 
temperature GaAs lasers, N. N. Winogradoff and H. K. Kessler, 
Solid-State Communications 5, 155-158 (Pergamon Press Ltd., 
London , England, 1967). 

The introduction of compensating donors into the p-type side of 
vapor grown GaAs lasers showed that: (1) At room temperature, the 
density of states in the conduction band of uncompensated p-type 
GaAs was too high to permit the production of an inverted popula- 
tion at practical current levels. (2) The presence of a suitable degree 
of compensation produced a low density of states tail which per- 
mitted lasing action to occur at low current densities. (3) High light 
power output and optimum lasing characteristics were produced by 
a critical degree of compensation. (4) This critical compensation 
could be .achieved by suitable thermal processing of solution grown 
and diffused lasers, but is more easily obtained by vapor phase 
epitaxy. 


Key Words: Gallium arsenide, compensation, tail banding, room 
temperature laser. 


Correlation of spherical aberration of axially symmetrical 
magnetic lenses, L. Marton, Rev. Sci. Instr. 38, No. 1, 130-131 
(Jan. 1967). 

It is pointed out that in a magnetic lens of rotational symmetry cor- 
rection of the spherical aberration is possible by abandoning parax- 
ial image formation and introducing a diamagnetic corrective ele- 
ment on the axis. 


Key Words: Electron lens, spherical aberration, diamagnetic cor- 
rective element. 


Correlation of large longitudinal deformations with different 
strain histories, L. J. Zapas and T. Craft, Rubber Chem, Technol. 
40, No. 2 , 506-516 (Mar. 1967). 

It is shown that the BKZ incompressible elastic fluid theory is in 
excellent agreement with experimental results obtained in simple 
extension. From single step stress-relaxation data, the stress-strain 
response for a number of other simple extension histories are calcu- 
lated from the theory and are compared with experiments. 


Key Words: BKZ theory, constant rate of strain, creep, elastic 
fluid, nonlinear behavior, polyisobutylene, polyvinylchloride, recov- 
ery, stress relaxation. 


Dislocation node determinations of the stacking fault energy 
in silver-tin alloys, A. W. Ruff, Jr. and L. K. Ives, Acta Met. 15, 
189-198 (Feb. 1967). 

Direct measurements by transmission electron microscopy on ex- 
tended dislocation nodes in alloys of tin in silver have led to values 
for the intrinsic stacking fault energy. The values decreased smoothly 
from 23 erg/cm? for pure silver to 4.5 erg/cm? for 7.8 at.% tin, near 
tne fec phase limit. The results are compared with previous determi- 
nations in other silver-base alloys, and sources of systematic error 
are discussed. A comparison of different theoretical analyses is 
included. 


Key Words: Dislocations, silver-tin alloy, stacking fault energy 
electron microscopy. 


Electrical feedthrough for pressures to 10 kbar, P. L. M. 
Heydemann, Rev. Sci. Instr. Laboratory Notes 38, No. 4, 558- 
559 (Apr. 1967). 

For ultrasonic experiments at high pressures and for measure- 
ments with manganin resistance gages an electrical feed-through 
was developed, which is both easily fabricated and convenient to 
use. 

The electrical feed-through (EFT) described here is designed for 
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low current applications at frequencies ranging from d.c. to over 
50 MHz. EFTs with one, two and four leads have been made. The 
insulation resistance between leads and between a lead and the 
shield is typically above 1000 MQ. The EFT can therefore be used 
in connection with very sensitive bridge circuits. The EFT is made 
from parts that are easily available. It is readily attached to and 
removed from various kinds of high pressure equipment by means 
of a standard type fitting. 


Key Words: High pressure, feed-through, electrical connector, 
ultrasonics, resistance gage. 


Ellipsemetric errors due to multiple reactions in mica 
quarter-wave plates, H. T. Yolken, R. M. Waxler, and J. Kruger, 
J. Opt. Soc. Am. 57, No. 2, 282-284 (Feb. 1967). 

Experiments performed to determine the source of the large dif- 
ferences in polarizer readings between the ellipsometric zones 
found that the origin of these polarizer reading differences lay in 
multiple reflections occurring in the quarter-wave plate. Coating 
the quarter-wave plate with antireflective layers markedly reduced 
the zone differences. 


Key Words: Ellipsometer, multiple reflections, quarter-wave plate. 


Electron monochromator design, C. E. Kuyatt and J. A. Simp- 
son, Rev. Sci. Instr. 38, No. 1, 103-111 (Jan. 1967). 

A study has been made of all the known factors which limit the 
performance of high resolution (0.07 to 0.01 eV FWHM) monochro- 
mators. These limiting factors have been incorporated into design 
equations for the optimum (maximum current output) monochro- 
mator. The conclusions are tested by performance measurements 
on a prototype instrument. The results require the introduction 
into the design equation of a new limiting factor, an anomalous 
energy spread in dense electron beams, which is empirically 
determined. 


Key Words: Anomalous energy spreads, electron beams, electron 
optical limitations, high-resolution electron monochromator, space 
charge, spherical deflectors. 


Estimate of extrinsic stacking-fault energies from disloca- 
tion configurations, W. F. Hartl, R. deWit, and R. E. Howard, 
J. Appl. Phys. 38, No. 1, 1-4 (Jan. 1967). 

In this paper we use straight line models of observed extrinsic nodes 
to estimate extrinsic stacking fault energies. Using isotropic elas- 
ticity theory and neglecting interactions we fd easily applicable 
expressions for the energy in terms of the inscribed radii of the 
nodes. We apply the results to be observed extended nodes in a 
silver—8 at.% tin alloy and find a ratio of éxtrinsic to intrinsic 
stacking fault energy of 3.0. 


Key Words: Stacking fault energy, extrinsic faults, dislocations, 
silver-tin alloy. 


Measurement of thermal neutron flux, R. S. Carter, C. O. 
Muehlhause, and V. W. Myers, JEEE Trans. Nucl. Sci. NS-14, 
No. 1, 414-415 (Feb. 1967). 

A new technique has been developed for the measurement of 
thermal neutron flux. Small glass beads (~ 1 mm radius) are made 
black to thermal neutrons by the inclusion of B'®. In addition, 
a I/v activator such as cobalt, dysprosium, or indium is incorporated 
in the glass to measure the fluence nv¥t). Following activation over 
a measured time interval, the absolute source strength of a Co 
bead is determined by y-y coincidence. The flux is directly pro- 
portional to the gamma ray source strength. The Dy and In beads 
are calibrated relative to a Co bead. 


Key Words: Thermal neutron flux, black glass beads, absolute 
gamma ray source strength. 


Second breakdown and crystallographic defects in transis- 
tors, H. A. Schafit, G. H. Schwuttke, and R. L. Ruggles, Jr., JEEE 
Trans. Electron Devices ED-13, No. 11, 738-742 (Nov. 1966). 
A study was conducted to determine the effect, on second breakdown 
in transistors, of some of the crystallographic defects that can 
develop during handling and fabricating procedures. X-ray dif- 
fraction microscopy was used to detect these defects. The suscep- 


tibility to second breakdown of about 1500 epitaxial planar silicon 
transistors diffused on two wafers was measured and the site of 
the current constriction of second breakdown was registered on 
most of these transistors. These data were then compared with the 
x-ray topographs. Any effect of the observed gross dislocations was 
masked by other factors, such as surface induced effects, that were 
not discernible in the x-ray topographs. Evidence was found, how- 
ever, to indicate that a faulted emitter structure that has been 
correlated with the emitter diffusion step may be of importance. 
The heating in the transistor at the second breakdown current 
constriction site to temperatures above 600 °C and recrystallized 
laser-induced melt regions smaller than about 30 ym in diameter 
produced no crystallographic changes nor any frozen in strain 
fields detectable in the x-ray topographs. 


Key Words: Transistors, second breakdown, dislocations, x-ray 
diffraction microscopy, semiconductors, failure mechanisms. 


Simultaneous diffraction with the three-circle diffrac- 
tometer, M. Zocchi and A. Santoro, Acta Cryst. 22, Pt. 3, 331- 
334 (Mar. 1967). 

The geometrical conditions of simultaneous diffraction with a 
three-circle diffractometer have been derived for all the crystal- 
systems and for any crystal orientation. 

It is shown that A-independent simultaneous diffraction takes 
place among reciprocal lattice points located on the same vertical 
net, that multiplicity is generally odd, and that it obeys special 
conditions in all cases in the cubic system and, for particular 
orientations and classes of reflection, in the hexagonal and tetragonal 
systems. 

Finally, procedures are suggested for avoiding taking intensity 
measurements under conditions of simultaneous diffraction. 
neutrons, simultaneous, three- 


Key Words: X-ray, 


circle, diffractometer. 


diffraction, 


Specific heat of natural rubber and other elastomers above 
the glass transition temperature, L. A. Wood and N. Bekkedahl, 
J. Polymer Sci. Letters 5, 169-175 (1967). 

A survey is given of literature values of the specific heat C, as a 
function of temperature (above the glass transition temperature) for 
the elastomers, including natural rubber, polyisobutylene, poly- 
butadiene, several butadiene-styrene copolymers, and a butadiene- 
acrylonitrile copolymer. The values of dC,,/dt (where t is the tempera- 
ture) found for natural rubber were abnormally low in the range 
—68 to +2 °C and abnormally high in the range 17° to 47 °C. The 
most probable explanation is that crystallization of the unvulcanized 
rubber was predominating in the first range, and melting of these 
crystals was predominating in the second range. 


Key Words: Specific heat, elastomers, rubber, polyisobutylene, 
polybutadiene, butadiene-styrene copolymers, butadiene-acryloni- 
trile copolymer , polyethylene, crystallization of rubber. 


The stacking-fault energy in alpha-silver-tin alloys, A. W. 
Ruff, Jr. and L. K. Ives, Can. J. Phys. 45, 787-795 (1967). 

Direct measurements by transmission electron microscopy on ex- 
tended dislocation nodes in alloys of tin in silver have led to values 
for the intrinsic stacking fault energy. The values decreased smoothly 
from 23 erg/cm* for pure silver to 4.5 erg/cm® for 7.8 at.% tin. The 
results are compared with previous determinations in other silver- 
base alloys. 


Key Words: Dislocations, stacking faults, silver-tin alloys, electron 
microscopy, stacking fault energy, dislocation nodes. 


Electrolysis of nonaqueous systems, A. Brenner, Book, Ad- 
vances in Electrochemistry and Electrochemical Engineering, Ed. 
C. W. Tobias, ch. V, pp. 205-248 (John Wiley & Sons, Inc., New 
York, N.Y., 1967). 

This review covers the electrodeposition of metals from organic 
solutions from ammonia and also the electrolysis of essentially pure 
organic and inorganic liquids. 


Galvanostalametry: A 
A. Brenner and J. L. Sligh, J. Electrochem. Soc. 114, No. 5, 461-465 
(May 1967). 


technique for chemical analysis, 


332 





A method utilizing a novel technique to indicate the end-point of a 
diffusion-controlled electro-reducible or oxidizable reaction has 
shown promising analytical capability. The method involves the drop- 
ping of a column of electrolyte, suspended under tensile stress, 
initiated by the electrolytic discharge of hydrogen or oxygen at an 
indicator electrode. It has been found that in the constant-current 
electrolysis of a suitable supporting electrolyte containing an electro- 
reactive ion that the transition time required for the column to drop is 
a function of the concentration of the latter. The use of this method 
as an analytical tool and its relationship to chronopotentiometry are 
discussed. 


Key Words: Chronopotentiometry, 
water, electrode reaction, 
galvanostalametry. 


analysis, 
limiting 


negative pressure of 
current density, diffusion, 


Measurements on the molecular nitrogen pulsed laser, 
T. Kasuya and D. R. Lide, Jr., Appl. Opt. 6, No. 1, 69-80 (Jan. 
1967). 

Accurate wavelength measurements of the pulsed laser emission 
from the nitrogen molecule have been made. About one hundred 
new emission lines which belong to the first and second positive 
systems were obtained and identified as to rotational quantum 
number. The relative intensities of the laser lines in the rotational 
band are analyzed. Some characteristics of the oscillations, including 
the recovery time after the laser pulse, are discussed. 


Key Words: 


ultraviolet. 


Electric discharge, laser, infrared, nitrogen, spectrum, 


Generalized integrating-sphere theory, D. G. Goebel, 
Opt. 6, No. 1, 125-128 (Jan. 1967). 

A general equation is developed for the efficiency of an integrating 
sphere with a nonuniform coating. The only assumptions used are 
that the interior is a perfect sphere and that all areas reflect per- 
fectly diffusely. Three special cases of the general equation are 
examined for the basic applications of integrating ‘spheres as 
mixing mechanisms in hemispherical-reflectance measurements 
and in absolute-reflectance techniques. 


Appl. 


Key Words: Integrating sphere, reflectance, absolute reflectance, 
sphere efficiency. 


High order correlations in a turbulent field, F. N. Frenkiel 
and P. S. Klebanoff, Phys. Fluids 10, No. 3, 507-520 (Mar. 1967). 
Higher-order time-correlations and the associated skewnesses and 
flatnesses were measured in a turbulent field downstream of a grid 
using high-speed computing techniques. The results were obtained 
using samples of 160,020 digitized data recorded at time intervals 
of 1/12,800 seconds during time periods approximately 12.5 seconds. 
Comparison is made between the measured correlations and the 
higher-order correlation curves corresponding to a Gaussian prob- 
ability density distribution of turbulent velocities. The departures 
from Gaussianity are shown, and non-Gaussian probability distri- 
butions are proposed which correspond considerably better to 
experimental reality. Several relations between correlation co- 
efficients of different orders are obtained for the non-Gaussian 
probability distributions and confirmed by comparison with the 
‘measured correlations, skewnesses, and flatnesses. 


Key Words: Turbulence, grid, correlations, higher-order, gaussian, 
non-gaussian, digital computer, hot-wire anemometer. 


Space-time correlations in a turbulent field, F. N. Frenkiel 
and P. S. Klebanoff (Proc. Symp. Dynamics of Fluids and Plasma, 
University of Maryland, Oct. 7-9, 1965), Book, Dynamics of Fluids 
and Plasmas, pp. 257-274 Vicuilianeie Press Inc., New York, N.Y., 
1966). 

Space-time correlations of the component of turbulent velocities 
along the direction of mean velocity were measured in a turbulent 
field downstream of a grid using high-speed computing techniques. 
These results provide new information in that the space-time evalua- 
tion is carried out to much higher orders than had previously been 
measured. Space-time correlations of even order up to the eighth, 
and of odd order up to the fifth are presented. The applicability of 
Taylor’s space-time approximation, i.e., the concept of a frozen 
pattern of turbulence moving with the mean velocity, to these higher 
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order correlations is examined. The higher even-order correlations 
are compared with the results obtained from the second-order cor- 
relations assuming a Gaussian distribution for the turbulent veloci- 
ties, and the different behavior of the various odd-order correlations 
is demonstrated. 


Key Words: Turbulence, wind tunnel, computer, hot-wire anemom- 
eter, correlations: even-order, odd-order, space-time, Taylor’s 
approximation. 


Terms, definitions, and symbols in reflectrometry, D. B. Judd, 
J. Opt. Soc. Am. 57, No. 4, 445-452 (Apr. 1967). 

Angular conditions ‘of incidence are described as hemispherical, 
conical, or directional; and the same adjectives are used to describe 
the angular conditions of collection. This classification of angular 
conditions leads to nine kinds of reflectance, and symbols for them 
are proposed in which 27; g; and 0, Mo refer to hemispherical, 
conical, and directional incidence; and 27; g’; and 0,, ®, refer to 
the corresponding kinds of collection. Use of the perfectly reflecting 
mirror and of the perfectly reflecting diffuser as reference standards 
in reflectometry is discussed. Three of the nine reflectance ratios, 
specimen to perfect diffuser, in which the collection is directional 
have already been named radiance [luminance] factor. It is pro- 
posed to differentiate them by angular condition of incidence. It is 
also proposed to name the other six ratios reflectance factor qualified 
by the same adjectives identifying the type of incidence and collec- 
tion as are used for reflectance. The interrelationships of these 18 
concepts are shown both by formulas for computing one from an- 
other and by diagrams indicating the process (integration, summa- 
tion, averaging, equality, reflectance of perfect diffuser and reci- 
procity) by which values of one concept may be computed from 
those of another. 


Key Words: Angular reflectometry conditions, collection angles, 
incidence angles, irradiation angles, reflectance terminology, re- 
flectometry, symbols for reflectometry, terminology for reflectance. 


Instrumentation and measurements, J. E. Gray and H. E. 
Bussey, Book, 1965 Digest of Literature on Dielectrics 29, ch. I, 
pp. 1-25 (National Academy of Sciences, National Research Council, 
Wash., D.C. 1967). 

A classified bibliography of abstracts from the literature of instru- 
mentation and measurement applicable to the study of dielectric 
materials has been prepared for publication as the first chapter of 
the Digest of Literature on Dielectrics for 1965. 


dielectric standards, 
constant, dielectric 


Key Words: Instrumentation, measurement, 
electrical insulation, permittivity, dielectric 
breakdown. 


Optical heterodyne measurement of neon laser’s millimeter 
wave difference frequency, J. L. Hall and W. W. Morey, Appl. 
Phys. Letters 10, No. 5, 152-155 (Mar. 1, 1967). 

We report detection and measurement of the millimeter wave dif- 
ference frequency between two near laser lines at 1.152 uw. The two 
spectral transitions, separated by 2.26 A, oscillate in pure neon in 
a single laser device, producing about 150 microwatts total power 
in several longitudinal modes. The many resulting 51.3-kmc beat 
frequencies have been studied by optical heterodyne techniques. 
A type of diode has been found which has enough microwave 
sensitivity and enough optical sensitivity to combine the optical 
detector and the microwave heterodyne functions in a single ele- 
ment. The preliminary value for the difference in frequency between 
the 2S.— 2P, and 2S,— 2P; transitions in neon at 150 millitorr is 
(51,360+ 150) M,.. This value is not in agreement with the value 
calculated from Vol. I of Sitterly’s Atomic Energy Levels (Ref. 12). 


Key Words: Neon, laser, optical heterodyne techniques, wave 
difference frequency. 


Progress in electron physics during the last 20 years, 
L. Marton, Third Czechoslovak Conf. Electronics and Vacuum 
Physics, Prague, Czechoslovakia, Sept. 23-28, 1965, pp. 17-25 
(Czechoslovak, 1967). 

Review of progress in electron physics with a short introduction 
covering pre-war period. Emphasis is laid on the developments of 
the post war years with the conclusion that although the greatest 


333 





accomplishments of electron physics were in the pre-war years, 
the last 20 years have seen very notable progress in the under- 
standing of the nature and of the interactions of the electrons. 


Key Words: Electron physics, history. 
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